Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Tue Nov 21, 2017 12:44 am




Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Author Message
 Post subject: NLM_LBC in ROMS 3.7
PostPosted: Fri Dec 09, 2016 11:09 am 
Offline

Joined: Tue Jan 14, 2014 3:02 pm
Posts: 16
Location: University of Copenhagen
I've been running a modified version of the upwelling case with (among other modifications), a uniform 1 m thick sediment bed in ROMS 3.7
When trying to change the sediment bed to have a non-uniform thickness in the entire domain, I added the following code to ana_sediment:
Code:
> #  if defined unised
>       DO j=JstrT,JendT
>    DO i=IstrT,IendT
>      DO k=1,Nbed
>             bed(i,j,k,iaged)=time(ng)
>             bed(i,j,k,ithck)=0.10_r8
>             bed(i,j,k,iporo)=0.90_r8
>             bed_frac(i,j,k,1)=1.0_r8
>           END DO
> !  Set exposed sediment layer properties.
>
>           bottom(i,j,irlen)=0.10_r8
>           bottom(i,j,irhgt)=0.01_r8
>           bottom(i,j,izdef)=Zob(ng)
> #  elif defined SHELFSED
>       DO j=JstrT,JendT
>    DO i=IstrT,IendT
>      val4=REAL(Lm(ng)+1-i,r8)
> !  Set bed layer properties.
>      IF(val4.le.20.0_r8) THEN
>             DO k=1,Nbed
>               bed(i,j,k,iaged)=time(ng)
>               bed(i,j,k,ithck)=0.10_r8
>               bed(i,j,k,iporo)=0.90_r8
>               bed_frac(i,j,k,1)=1.0_r8
>             END DO
> !  Set exposed sediment layer properties.
>
>             bottom(i,j,irlen)=0.10_r8
>             bottom(i,j,irhgt)=0.01_r8
>             bottom(i,j,izdef)=Zob(ng)
>           ELSE IF(val4.le.25.0_r8) THEN 
>             DO k=1,Nbed
>               bed(i,j,k,iaged)=time(ng)
>               bed(i,j,k,ithck)=0.10_r8*cos((val4-21)*pi/8)
>               bed(i,j,k,iporo)=0.90_r8
>               bed_frac(i,j,k,1)=1.0_r8
>             END DO
> !  Set exposed sediment layer properties.
>
>             bottom(i,j,irlen)=0.10_r8
>             bottom(i,j,irhgt)=0.01_r8
>             bottom(i,j,izdef)=Zob(ng)
>           ELSE
>             DO k=1,Nbed
>               bed(i,j,k,iaged)=time(ng)
>               bed(i,j,k,ithck)=0.0_r8
>               bed(i,j,k,iporo)=0.9_r8
>               bed_frac(i,j,k,1)=1.0_r8
>             END DO
> !  Set exposed sediment layer properties.
>
>             bottom(i,j,irlen)=0.0_r8
>             bottom(i,j,irhgt)=0.0_r8
>             bottom(i,j,izdef)=Zob(ng)
>           END IF
>         END DO
>       END DO
> #  endif


My model runs in three phases with different forcing, so I'm restarting from the RST-file a lot. This has never been a problem until now.
Since making this change, I've been getting the LBC_GETATT error when trying to restart from a previous solution.
I have tried changing the ana_sediment.h file back, I've tried overwriting it with a different file that I know to work, but this does not solve the problem.

Does anyone have experience with this issue? I guess I could just activate the #define NO_LBC_ATT option, but it seems like this should not be necessary when running 3.7


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group