Changeset 11442 for branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90

Timestamp:

2019-08-16T12:32:43+02:00 (5 years ago)

Author:

mattmartin

Message:

Introduction of stochastic physics in NEMO, based on Andrea Storto's code.
For details, see ticket ​https://code.metoffice.gov.uk/trac/utils/ticket/251.

File:

• branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90

@@ -30 +30 @@
    USE wrk_nemo        ! Memory Allocation
    USE timing          ! Timing
+   USE stopack
 
    IMPLICIT NONE
@@ -38 +39 @@
 
    INTEGER ::   nldf = -2   ! type of lateral diffusion used defined from ln_dynldf_... namlist logicals)
+
+#if defined key_dynldf_c3d
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ahm10, ahm20, ahm30, ahm40
+#endif
+#if defined key_dynldf_c2d
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:  ) :: ahm10, ahm20, ahm30, ahm40
+#endif
 
    !! * Substitutions
@@ -67 +75 @@
          ztrdv(:,:,:) = va(:,:,:) 
       ENDIF
+
+#if defined key_dynldf_c3d
+         IF( kt == nit000 .AND. ln_stopack .AND. &
+            & ( (nn_spp_ahm1 + nn_spp_ahm2) > 0 ) ) THEN
+             ALLOCATE ( ahm10(jpi,jpj,jpk), ahm20(jpi,jpj,jpk) )
+             ALLOCATE ( ahm30(jpi,jpj,jpk), ahm40(jpi,jpj,jpk) )
+             ahm10 = ahm1
+             ahm20 = ahm2
+             ahm30 = ahm3
+             ahm40 = ahm4
+         ENDIF
+#endif
+#if defined key_dynldf_c2d
+         IF( kt == nit000 .AND. ln_stopack .AND. &
+            & ( (nn_spp_ahm1 + nn_spp_ahm2) > 0 ) ) THEN
+             ALLOCATE ( ahm10(jpi,jpj), ahm20(jpi,jpj) )
+             ALLOCATE ( ahm30(jpi,jpj), ahm40(jpi,jpj) )
+             ahm10 = ahm1
+             ahm20 = ahm2
+             ahm30 = ahm3
+             ahm40 = ahm4
+         ENDIF
+#endif
+
+#if defined key_traldf_c3d || defined key_traldf_c2d
+         IF( ln_stopack ) THEN
+            IF( nn_spp_ahm1 > 0 ) THEN
+               IF( ln_dynldf_lap ) THEN
+                  ahm1 = ahm10
+                  CALL spp_ahm(kt, ahm1, nn_spp_ahm1, rn_ahm1_sd, jk_spp_ahm1)
+               ENDIF
+               IF( ln_dynldf_bilap ) THEN
+                  ahm3 = ahm30
+                  CALL spp_ahm(kt, ahm3, nn_spp_ahm1, rn_ahm1_sd, jk_spp_ahm3)
+               ENDIF
+            ENDIF
+            IF( nn_spp_ahm2 > 0 ) THEN
+               IF( ln_dynldf_lap ) THEN
+               ahm2 = ahm20
+               CALL spp_ahm(kt, ahm2, nn_spp_ahm2, rn_ahm2_sd, jk_spp_ahm2)
+               ENDIF
+               IF( ln_dynldf_bilap ) THEN
+                  ahm4 = ahm40
+                  CALL spp_ahm(kt, ahm4, nn_spp_ahm2, rn_ahm2_sd, jk_spp_ahm4)
+               ENDIF
+            ENDIF
+         ENDIF
+#endif
 
       SELECT CASE ( nldf )                       ! compute lateral mixing trend and add it to the general trend