Changeset 15676
- Timestamp:
- 2022-01-27T16:28:43+01:00 (3 years ago)
- 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 39 39 !: 4 = Jan off-shelf 40 40 !: 5 = Siddorn & Furner 2013 41 !: 6 = Shchepetkin & McWilliams 2003 41 42 REAL(wp), PUBLIC :: rn_sal_sf, rn_mld_sf, rn_maxdep_sf, rn_c0_sf, rn_c1_sf !: Siddorn & Furner 2013 42 43 !: idealised profile … … 104 105 WRITE(numout,*) ' Type of intial profile nn_tsd_type = ', nn_tsd_type 105 106 IF (nn_tsd_type == 5) THEN 106 WRITE(numout,*) ' Siddorn and Furner 2013 id alised initial profile'107 WRITE(numout,*) ' Siddorn and Furner 2013 idealised initial profile' 107 108 WRITE(numout,*) ' rn_sal_sf = ', rn_sal_sf 108 109 WRITE(numout,*) ' rn_mld_sf = ', rn_mld_sf … … 110 111 WRITE(numout,*) ' rn_c0_sf = ', rn_c0_sf 111 112 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' 112 116 ENDIF 113 117 ENDIF -
NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE/USR/usrdef_istate.F90
r15366 r15676 170 170 ! b0 = -1.1 171 171 ! 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. 173 180 END SELECT 174 181 … … 181 188 END DO 182 189 END DO 183 184 ELSE ! Use an analytic function 190 191 ELSE IF (nn_tsd_type == 5) THEN ! Siddorn & Furner 2013 analytical function 192 185 193 DO jk = 1, jpk 186 194 DO jj = 1, jpj … … 192 200 END DO 193 201 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 194 215 END IF 195 216
Note: See TracChangeset
for help on using the changeset viewer.