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 6667 for branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/usrdef.F90 – NEMO

Ignore:
Timestamp:
2016-06-06T07:57:00+02:00 (8 years ago)
Author:
gm
Message:

#1692 - branch SIMPLIF_2_usrdef: reduced domain_cfg.nc file: GYRE OK using usrdef or reading file

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/usrdef.F90

    r6624 r6667  
    218218   END SUBROUTINE usr_def_hgr 
    219219 
    220  
    221    SUBROUTINE usr_def_zgr() 
    222        ! subroutine for vertical grid 
    223    END SUBROUTINE usr_def_zgr 
    224  
    225220  
    226221   SUBROUTINE usr_def_ini( pts ) 
    … …  
    248243            DO ji = 1, jpi 
    249244               pts(ji,jj,jk,jp_tem) =  (  (  16. - 12. * TANH( (gdept_0(ji,jj,jk) - 400) / 700 ) )   & 
    250                     &           * (-TANH( (500-gdept_0(ji,jj,jk)) / 150. ) + 1.) / 2.               & 
    251                     &           + ( 15. * ( 1. - TANH( (gdept_0(ji,jj,jk)-50.) / 1500.) )        & 
    252                     &           - 1.4 * TANH((gdept_0(ji,jj,jk)-100.) / 100.)                    & 
    253                     &           + 7.  * (1500. - gdept_0(ji,jj,jk) ) / 1500.)                     & 
     245                    &           * (-TANH( (500-gdept_0(ji,jj,jk)) / 150. ) + 1.) / 2.                & 
     246                    &           + ( 15. * ( 1. - TANH( (gdept_0(ji,jj,jk)-50.) / 1500.) )            & 
     247                    &           - 1.4 * TANH((gdept_0(ji,jj,jk)-100.) / 100.)                        & 
     248                    &           + 7.  * (1500. - gdept_0(ji,jj,jk) ) / 1500.)                        & 
    254249                    &           * (-TANH( (gdept_0(ji,jj,jk) - 500.) / 150.) + 1.) / 2.  ) * tmask(ji,jj,jk) 
    255250 
    256251               pts(ji,jj,jk,jp_sal) =  (  (  36.25 - 1.13 * TANH( (gdept_0(ji,jj,jk) - 305) / 460 ) )  & 
    257                     &         * (-TANH((500. - gdept_0(ji,jj,jk)) / 150.) + 1.) / 2              & 
    258                     &         + ( 35.55 + 1.25 * (5000. - gdept_0(ji,jj,jk)) / 5000.          & 
    259                     &         - 1.62 * TANH( (gdept_0(ji,jj,jk) - 60.  ) / 650. )             & 
    260                     &         + 0.2  * TANH( (gdept_0(ji,jj,jk) - 35.  ) / 100. )             & 
    261                     &         + 0.2  * TANH( (gdept_0(ji,jj,jk) - 1000.) / 5000.) )           & 
    262                     &         * (-TANH((gdept_0(ji,jj,jk) - 500.) / 150.) + 1.) / 2  ) * tmask(ji,jj,jk) 
     252                    &         * (-TANH((500. - gdept_0(ji,jj,jk)) / 150.) + 1.) / 2                  & 
     253                    &         + ( 35.55 + 1.25 * (5000. - gdept_0(ji,jj,jk)) / 5000.                 & 
     254                    &         - 1.62 * TANH( (gdept_0(ji,jj,jk) - 60.  ) / 650. )                    & 
     255                    &         + 0.2  * TANH( (gdept_0(ji,jj,jk) - 35.  ) / 100. )                    & 
     256                    &         + 0.2  * TANH( (gdept_0(ji,jj,jk) - 1000.) / 5000.) )                  & 
     257                    &         * (-TANH( (gdept_0(ji,jj,jk) - 500.) / 150.) + 1.) / 2  ) * tmask(ji,jj,jk) 
    263258            END DO 
    264259         END DO 
Note: See TracChangeset for help on using the changeset viewer.