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 15676 – NEMO

Changeset 15676


Ignore:
Timestamp:
2022-01-27T16:28:43+01:00 (2 years ago)
Author:
dbruciaferri
Message:

Adding Shchepetkin & McWilliams? 2003 initial profile

Location:
NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE/DOM/dtatsd.F90

    r15366 r15676  
    3939                                       !: 4 = Jan off-shelf 
    4040                                       !: 5 = Siddorn & Furner 2013 
     41                                       !: 6 = Shchepetkin & McWilliams 2003 
    4142   REAL(wp), PUBLIC ::   rn_sal_sf, rn_mld_sf, rn_maxdep_sf, rn_c0_sf, rn_c1_sf !: Siddorn & Furner 2013  
    4243                                                                                !: idealised profile  
     
    104105         WRITE(numout,*) '      Type of intial profile nn_tsd_type = ', nn_tsd_type 
    105106         IF (nn_tsd_type == 5) THEN 
    106             WRITE(numout,*) '      Siddorn and Furner 2013 idalised initial profile' 
     107            WRITE(numout,*) '      Siddorn and Furner 2013 idealised initial profile' 
    107108            WRITE(numout,*) '            rn_sal_sf    = ', rn_sal_sf 
    108109            WRITE(numout,*) '            rn_mld_sf    = ', rn_mld_sf 
     
    110111            WRITE(numout,*) '            rn_c0_sf     = ', rn_c0_sf 
    111112            WRITE(numout,*) '            rn_c1_sf     = ', rn_c1_sf 
     113         ENDIF 
     114         IF (nn_tsd_type == 6) THEN 
     115            WRITE(numout,*) '      Shchepetkin and McWilliams 2003 idealised initial profile' 
    112116         ENDIF 
    113117      ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE/USR/usrdef_istate.F90

    r15366 r15676  
    170170           ! b0 = -1.1  
    171171           ! b1 = 0.01  
    172            ! b2 = 34.85  
     172           ! b2 = 34.85 
     173         CASE(6) ! Analytic T in the tropical Atlantic ocean 
     174            rn_sal_sf = 35 
     175            a0 = 0.2  
     176            a1 = 6.0  
     177            a2 = 20.  
     178            b0 = 2500. 
     179            b1 = 250. 
    173180      END SELECT 
    174181 
     
    181188            END DO 
    182189         END DO 
    183      
    184       ELSE  ! Use an analytic function 
     190 
     191      ELSE IF (nn_tsd_type == 5) THEN ! Siddorn & Furner 2013 analytical function 
     192 
    185193         DO jk = 1, jpk 
    186194            DO jj = 1, jpj 
     
    192200            END DO 
    193201         END DO 
     202 
     203      ELSE IF (nn_tsd_type == 6) THEN ! Shchepetkin and McWilliams 2003 analytical function 
     204 
     205         DO jk = 1, jpk 
     206            DO jj = 1, jpj 
     207               DO ji = 1, jpi 
     208                  pts(ji,jj,jk,jp_sal) = rn_sal_sf !* ptmask(ji,jj,jk) 
     209                  pts(ji,jj,jk,jp_tem) = a0 + a1*EXP(-pdept(ji,jj,jk)/b0) & 
     210                                            + a2*EXP(-pdept(ji,jj,jk)/b1) !* ptmask(ji,jj,jk)  
     211               END DO 
     212            END DO 
     213         END DO 
     214 
    194215      END IF 
    195216 
Note: See TracChangeset for help on using the changeset viewer.