Changeset 4457


Ignore:
Timestamp:
2014-02-04T13:16:53+01:00 (7 years ago)
Author:
trackstand2
Message:

Add mbkmax and z-first opt to dynhpg (hpg_sco)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

    r3211 r4457  
    405405      ENDIF 
    406406 
    407       ! Surface value 
     407#if ! defined key_z_first 
     408     ! Surface value 
    408409      DO jj = 2, jpjm1 
    409410         DO ji = fs_2, fs_jpim1   ! vector opt.    
     
    423424         END DO   
    424425      END DO    
    425              
     426#endif 
     427 
    426428      ! interior value (2=<jk=<jpkm1) 
    427429#if defined key_z_first 
    428430      DO jj = 2, jpjm1      
    429431         DO ji = 2, jpim1 
    430             DO jk = 2, jpkm1                                   
     432 
     433            ! Surface value 
     434 
     435            ! hydrostatic pressure gradient along s-surfaces 
     436            zhpi(ji,jj,1) = zcoef0 / e1u(ji,jj) * ( fse3w(ji+1,jj  ,1) * ( znad + rhd(ji+1,jj  ,1) )   & 
     437               &                                  - fse3w(ji  ,jj  ,1) * ( znad + rhd(ji  ,jj  ,1) ) ) 
     438            zhpj(ji,jj,1) = zcoef0 / e2v(ji,jj) * ( fse3w(ji  ,jj+1,1) * ( znad + rhd(ji  ,jj+1,1) )   & 
     439               &                                  - fse3w(ji  ,jj  ,1) * ( znad + rhd(ji  ,jj  ,1) ) ) 
     440            ! s-coordinate pressure gradient correction 
     441            zuap = -zcoef0 * ( rhd   (ji+1,jj,1) + rhd   (ji,jj,1) + 2._wp * znad )   & 
     442               &           * ( fsde3w(ji+1,jj,1) - fsde3w(ji,jj,1) ) / e1u(ji,jj) 
     443            zvap = -zcoef0 * ( rhd   (ji,jj+1,1) + rhd   (ji,jj,1) + 2._wp * znad )   & 
     444               &           * ( fsde3w(ji,jj+1,1) - fsde3w(ji,jj,1) ) / e2v(ji,jj) 
     445            ! add to the general momentum trend 
     446            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) + zuap 
     447            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) + zvap 
     448 
     449            DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1                                   
    431450#else 
    432451      DO jk = 2, jpkm1                                   
Note: See TracChangeset for help on using the changeset viewer.