New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 462 – NEMO

Changeset 462


Ignore:
Timestamp:
2006-05-10T19:20:08+02:00 (18 years ago)
Author:
opalod
Message:

nemo_v1_update_053:RB: light update of trends part (supppress key_partial_steps)

Location:
trunk/NEMO/OPA_SRC/TRD
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/TRD/trdmod.F90

    r247 r462  
    223223                  DO jj = 2, jpjm1 
    224224                     DO ji = fs_2, fs_jpim1   ! vector opt. 
    225 # if defined key_s_coord || defined key_partial_steps 
     225# if defined key_zco 
     226                        zbtr = zbtr2(ji,jj) 
     227                        zfui = 0.5 * e2u(ji  ,jj) * zun(ji,  jj,jk) 
     228                        zfui1= 0.5 * e2u(ji-1,jj) * zun(ji-1,jj,jk) 
     229# else          
    226230                        zbtr = zbtr2(ji,jj) / fse3t(ji,jj,jk) 
    227  
    228231                        zfui = 0.5 * e2u(ji  ,jj) * fse3u(ji,  jj,jk) * zun(ji,  jj,jk) 
    229232                        zfui1= 0.5 * e2u(ji-1,jj) * fse3u(ji-1,jj,jk) * zun(ji-1,jj,jk) 
    230 # else          
    231                         zbtr = zbtr2(ji,jj) 
    232  
    233                         zfui = 0.5 * e2u(ji  ,jj) * zun(ji,  jj,jk) 
    234                         zfui1= 0.5 * e2u(ji-1,jj) * zun(ji-1,jj,jk) 
    235233# endif 
    236234                        ! save i- advective trend  
     
    250248                  DO jj = 2, jpjm1 
    251249                     DO ji = fs_2, fs_jpim1   ! vector opt. 
    252 # if defined key_s_coord || defined key_partial_steps 
     250# if defined key_zco 
     251                        zbtr = zbtr2(ji,jj) 
     252                        zfvj = 0.5 * e1v(ji,jj  ) * zvn(ji,jj  ,jk) 
     253                        zfvj1= 0.5 * e1v(ji,jj-1) * zvn(ji,jj-1,jk) 
     254# else          
    253255                        zbtr = zbtr2(ji,jj) / fse3t(ji,jj,jk) 
    254  
    255256                        zfvj = 0.5 * e1v(ji,jj  ) * fse3v(ji,jj  ,jk) * zvn(ji,jj  ,jk) 
    256257                        zfvj1= 0.5 * e1v(ji,jj-1) * fse3v(ji,jj-1,jk) * zvn(ji,jj-1,jk) 
    257 # else          
    258                         zbtr = zbtr2(ji,jj) 
    259  
    260                         zfvj = 0.5 * e1v(ji,jj  ) * zvn(ji,jj  ,jk) 
    261                         zfvj1= 0.5 * e1v(ji,jj-1) * zvn(ji,jj-1,jk) 
    262258# endif 
    263259                        ! save j- advective trend  
     
    322318                  DO jj = 2, jpjm1       
    323319                     DO ji = fs_2, fs_jpim1   ! vector opt. 
    324 # if defined key_s_coord || defined key_partial_steps 
     320# if defined key_zco 
     321                        zbtr = zbtr2(ji,jj) 
     322                        zfui = e2u(ji  ,jj) * zun(ji,  jj,jk)   & 
     323                           & - e2u(ji-1,jj) * zun(ji-1,jj,jk) 
     324# else       
    325325                        zbtr = zbtr2(ji,jj) / fse3t(ji,jj,jk) 
    326326                        zfui =  e2u(ji  ,jj) * fse3u(ji,  jj,jk) * zun(ji,  jj,jk)   & 
    327327                           & -  e2u(ji-1,jj) * fse3u(ji-1,jj,jk) * zun(ji-1,jj,jk) 
    328 # else       
    329                         zbtr = zbtr2(ji,jj) 
    330                         zfui = e2u(ji  ,jj) * zun(ji,  jj,jk)   & 
    331                            & - e2u(ji-1,jj) * zun(ji-1,jj,jk) 
    332328# endif 
    333329                        ! save i- advective trend  
     
    347343                  DO jj = 2, jpjm1       
    348344                     DO ji = fs_2, fs_jpim1   ! vector opt. 
    349 # if defined key_s_coord || defined key_partial_steps 
     345# if defined key_zco 
     346                        zbtr = zbtr2(ji,jj) 
     347                        zfvj = e1v(ji,jj  ) * zvn(ji,jj  ,jk)   & 
     348                           & - e1v(ji,jj-1) * zvn(ji,jj-1,jk) 
     349# else       
    350350                        zbtr = zbtr2(ji,jj) / fse3t(ji,jj,jk) 
    351351                        zfvj =  e1v(ji,jj  ) * fse3v(ji,jj  ,jk) * zvn(ji,jj  ,jk)   & 
    352352                           & -  e1v(ji,jj-1) * fse3v(ji,jj-1,jk) * zvn(ji,jj-1,jk) 
    353 # else       
    354                         zbtr = zbtr2(ji,jj) 
    355                         zfvj = e1v(ji,jj  ) * zvn(ji,jj  ,jk)   & 
    356                            & - e1v(ji,jj-1) * zvn(ji,jj-1,jk) 
    357353# endif 
    358354                        ! save j- advective trend  
  • trunk/NEMO/OPA_SRC/TRD/trdvor.F90

    r352 r462  
    2323   USE daymod          ! calandar 
    2424   USE dianam          ! build the name of file (routine) 
    25    USE ldfslp          ! iso-neutral slopes 
    2625   USE zdfmxl          ! mixed layer depth 
    2726   USE ioipsl          ! NetCDF library 
Note: See TracChangeset for help on using the changeset viewer.