Changeset 15676 for NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE/USR/usrdef_istate.F90

Timestamp:

2022-01-27T16:28:43+01:00 (2 years ago)

Author:

dbruciaferri

Message:

Adding Shchepetkin & McWilliams? 2003 initial profile

File:

• NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE/USR/usrdef_istate.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_hpge_ovf/src/OCE/USR/usrdef_istate.F90

@@ -170 +170 @@
            ! b0 = -1.1 
            ! b1 = 0.01 
-           ! b2 = 34.85 
+           ! b2 = 34.85
+         CASE(6) ! Analytic T in the tropical Atlantic ocean
+            rn_sal_sf = 35
+            a0 = 0.2 
+            a1 = 6.0 
+            a2 = 20. 
+            b0 = 2500.
+            b1 = 250.
       END SELECT
 
@@ -181 +188 @@
             END DO
          END DO
-    
-      ELSE  ! Use an analytic function
+
+      ELSE IF (nn_tsd_type == 5) THEN ! Siddorn & Furner 2013 analytical function
+
          DO jk = 1, jpk
             DO jj = 1, jpj
@@ -192 +200 @@
             END DO
          END DO
+
+      ELSE IF (nn_tsd_type == 6) THEN ! Shchepetkin and McWilliams 2003 analytical function
+
+         DO jk = 1, jpk
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  pts(ji,jj,jk,jp_sal) = rn_sal_sf !* ptmask(ji,jj,jk)
+                  pts(ji,jj,jk,jp_tem) = a0 + a1*EXP(-pdept(ji,jj,jk)/b0) &
+                                            + a2*EXP(-pdept(ji,jj,jk)/b1) !* ptmask(ji,jj,jk) 
+               END DO
+            END DO
+         END DO
+
       END IF