Changeset 13239

Timestamp:

2020-07-03T12:30:58+02:00 (3 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0.1_GM_rossby_radius_cutoff : Yet another shape option, and new slope diagnostic.

Location:

NEMO/branches/UKMO/NEMO_4.0.1_GM_rossby_radius_cutoff

Files:

• cfgs/SHARED/field_def_nemo-oce.xml (modified) (1 diff)
• src/OCE/LDF/ldftra.F90 (modified) (5 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_GM_rossby_radius_cutoff/cfgs/SHARED/field_def_nemo-oce.xml

@@ -71 +71 @@
         <field id="rhop"         long_name="potential density (sigma0)"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D" />
         <field id="N_2d"         long_name="Depth-mean N"                                                              unit="m/s" />
+        <field id="modslp"       long_name="sqrt( slpi^2 + slpj^2 )"                                                   unit="1"      grid_ref="grid_T_3D" />
         <field id="RossRad"      long_name="Rossby radius"                                                             unit="m" />
         <field id="RossRadlim"   long_name="Rossby radius (limited)"                                                   unit="m" />

NEMO/branches/UKMO/NEMO_4.0.1_GM_rossby_radius_cutoff/src/OCE/LDF/ldftra.F90

@@ -647 +647 @@
       REAL(wp) ::   zfw, ze3w, zn2, z1_f20, zaht, zaht_min, zzaei, z2_3    ! local scalars
       REAL(wp), DIMENSION(jpi,jpj) ::   zn, zah, zhw, zRo, zRo_lim, zTclinic_recip, zaeiw, zratio   ! 2D workspace
+      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zmodslp ! 3D workspace
       !!----------------------------------------------------------------------
       !
       zn (:,:) = 0._wp        ! Local initialization
+      zmodslp(:,:,:) = 0._wp
       zhw(:,:) = 5._wp
       zah(:,:) = 0._wp
@@ -689 +691 @@
                   ! T^-1 = sqrt(m_jpk(N^2*(r1^2+r2^2)*e3w))
                   ze3w = e3w_n(ji,jj,jk) * tmask(ji,jj,jk)
-                  zah(ji,jj) = zah(ji,jj) + zn2 * ( wslpi(ji,jj,jk) * wslpi(ji,jj,jk)   &
-                     &                            + wslpj(ji,jj,jk) * wslpj(ji,jj,jk) ) * ze3w
+                  zmodslp(ji,jj,jk) =  wslpi(ji,jj,jk) * wslpi(ji,jj,jk)   &
+                     &               + wslpj(ji,jj,jk) * wslpj(ji,jj,jk)
+                  zah(ji,jj) = zah(ji,jj) + zn2 * zmodslp(ji,jj,jk) * ze3w
                   zhw(ji,jj) = zhw(ji,jj) + ze3w
                END DO
@@ -709 +712 @@
       END DO
       IF( iom_use('N_2d') ) CALL iom_put('N_2d',zn(:,:)/zhw(:,:))
+      IF( iom_use('modslp') ) CALL iom_put('modslp',SQRT(zmodslp(:,:,:)) )
       CALL iom_put('RossRad',zRo)
       CALL iom_put('RossRadlim',zRo_lim)
@@ -715 +719 @@
       z1_f20 = 1._wp / (  2._wp * omega * sin( rad * 20._wp )  )
       z2_3 = 2._wp/3._wp
-      DO jj = 2, jpjm1
-         DO ji = fs_2, fs_jpim1   ! vector opt.
-            zzaei = MIN( 1._wp, ABS( ff_t(ji,jj) * z1_f20 ) ) * zaeiw(ji,jj)     ! tropical decrease
-            zaeiw(ji,jj) = MIN( zzaei, paei0 )
-         END DO
-      END DO
 
       SELECT CASE(nn_ldfeiv_shape)
+         CASE(0) !! Standard shape applied - decrease in tropics and cap. 
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  zzaei = MIN( 1._wp, ABS( ff_t(ji,jj) * z1_f20 ) ) * zaeiw(ji,jj)     ! tropical decrease
+                  zaeiw(ji,jj) = MIN( zzaei, paei0 )
+               END DO
+            END DO
+
          CASE(1) !! Abrupt cut-off on Rossby radius:
 ! JD : modifications here to introduce scaling by local rossby radius of deformation vs local grid scale
@@ -761 +767 @@
                END DO
             END DO
+
+         CASE(5) !! Rossby radius ramp type 1 applied to Treguier et al coefficient rather than cap:
+                 !! Note the ramp is RR/GS=[2.0,1.0] (not [2.0,0.5] as for cases 2,3) and we ramp up 
+                 !! to 5% of the Treguier et al coefficient, aiming for peak values of around 100m2/s
+                 !! at high latitudes rather than 2000m2/s which is what you get in eORCA025 with an 
+                 !! uncapped coefficient.
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  zratio(ji,jj) = zRo(ji,jj)/MIN(e1t(ji,jj),e2t(ji,jj))
+                  zaeiw(ji,jj) = MAX( 0._wp, MIN( 1._wp, 2._wp - zratio(ji,jj) ) ) * 0.05 * zaeiw(ji,jj)
+                  zaeiw(ji,jj) = MIN( zaeiw(ji,jj), paei0 )
+               END DO
+            END DO
+            CALL iom_put('RR_GS',zratio)
+
+         CASE DEFAULT
+               CALL ctl_stop('ldf_eiv: Unrecognised option for nn_ldfeiv_shape.')         
 
       END SELECT