Changeset 3364 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jun 10, 2024, 4:38:06 PM (7 months ago)

Author:

bhatnags

Message:

Generic-PCM: Added dt_gm (Gent-McWilliams?) in the calculation for fluxgrdocean (GND) and in ocean_slab_get_vars.
SB

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• ocean_slab_mod.F90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (9 diffs)

trunk/LMDZ.GENERIC/libf/phystd/ocean_slab_mod.F90

@@ -1090 +1090 @@
 !
   SUBROUTINE ocean_slab_get_vars(ngrid, tslab_loc, tice_loc, seaice_loc, &
-                                 flux_g_loc, dt_hdiff_loc,dt_ekman_loc)
+                                 flux_g_loc, dt_hdiff_loc,dt_ekman_loc,dt_gm_loc)
 
 ! "Get some variables from module ocean_slab_mod"
@@ -1101 +1101 @@
     REAL, DIMENSION(ngrid,nslay), INTENT(OUT) :: dt_hdiff_loc ! [in W/m2]
     REAL, DIMENSION(ngrid,nslay), INTENT(OUT) :: dt_ekman_loc ! [in W/m2]
+    REAL, DIMENSION(ngrid,nslay), INTENT(OUT) :: dt_gm_loc ! [in W/m2]
     INTEGER :: i
 
@@ -1111 +1112 @@
     dt_hdiff_loc(:,:)=0.
     dt_ekman_loc(:,:)=0.
+    dt_gm_loc(:,:)=0.
 
 !!      dt_hdiff_loc(:,i) = dt_hdiff(:,i)*slabh(i)*1000.*4228. !Convert en W/m2
@@ -1125 +1127 @@
         ENDDO
     ENDIF
+    IF (slab_gm) THEN
+        DO i=1,nslay
+           dt_gm_loc(:,i) = dt_gm(:,i)
+        ENDDO
+    ENDIF
 
 

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -327 +327 @@
       real zdtrain_generic(ngrid,nlayer)                      ! Rain_generic routine.
       real dtmoist(ngrid,nlayer)                              ! Moistadj routine.
-      real dt_ekman(ngrid,nslay), dt_hdiff(ngrid,nslay) ! Slab_ocean routine.
+      real dt_ekman(ngrid,nslay), dt_hdiff(ngrid,nslay), dt_gm(ngrid,nslay) ! Slab_ocean routine.
       real zdtsw1(ngrid,nlayer), zdtlw1(ngrid,nlayer)         ! Callcorrk routine.
       real zdtchim(ngrid,nlayer)                              ! Calchim routine.
@@ -669 +669 @@
            ! ensure that "slab" variables match "physiq" ones
            call ocean_slab_get_vars(ngrid, tslab, tsea_ice, sea_ice, flux_g, &
-                      dt_hdiff, dt_ekman)
+                      dt_hdiff, dt_ekman, dt_gm)
            
          endif ! end of 'ok_slab_ocean'.
@@ -689 +689 @@
                   fluxgrd(ig)=0.
                   ! Dividing by cell area to have flux in W/m2
-                  fluxgrdocean(ig)=flux_g(ig)+(1-pctsrf_sic(ig))*(dt_hdiff(ig,1)+dt_ekman(ig,1))/cell_area(ig)
+                  fluxgrdocean(ig)=flux_g(ig)+(1-pctsrf_sic(ig))*(dt_hdiff(ig,1)+dt_ekman(ig,1)+dt_gm(ig,1))/cell_area(ig)
                   do iq=1,nsoilmx
                     tsoil(ig,iq)=tsurf(ig)
@@ -1992 +1992 @@
 
             call ocean_slab_get_vars(ngrid, tslab, tsea_ice, sea_ice, flux_g, &
-                      dt_hdiff, dt_ekman)
+                      dt_hdiff, dt_ekman, dt_gm)
                       
 !! sea_ice defines the sea ice thickness in kg/m2
@@ -2100 +2100 @@
 !!               fluxgrdocean(ig)=pctsrf_sic(ig)*flux_g(ig)+(1-pctsrf_sic(ig))*(dt_hdiff(ig,1)+dt_ekman(ig,1))
                ! Dividing by cell area to have flux in W/m2
-               fluxgrdocean(ig)=flux_g(ig)+(1-pctsrf_sic(ig))*(dt_hdiff(ig,1)+dt_ekman(ig,1))/cell_area(ig)
+               fluxgrdocean(ig)=flux_g(ig)+(1-pctsrf_sic(ig))*(dt_hdiff(ig,1)+dt_ekman(ig,1)+dt_gm(ig,1))/cell_area(ig)
+!               fluxgrdocean(ig)=(flux_g(ig)-fluxabs_sw(ig))+(1-pctsrf_sic(ig))*(dt_hdiff(ig,1)+dt_ekman(ig,1)+dt_gm(ig,1))/cell_area(ig)
                do iq=1,nsoilmx
                   tsoil(ig,iq)=tsurf(ig)
@@ -2389 +2390 @@
             ! fetch "ocean variables" to ensure they are stored
             call ocean_slab_get_vars(ngrid, tslab, tsea_ice, sea_ice, flux_g, &
-                                     dt_hdiff, dt_ekman)
+                                     dt_hdiff, dt_ekman, dt_gm)
 
             call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq, &
@@ -2474 +2475 @@
 
          if (ok_slab_ocean) then
-            call wstats(ngrid,"dt_hdiff1","dt_hdiff1","K/s",2,dt_hdiff(:,1))
-            call wstats(ngrid,"dt_hdiff2","dt_hdiff2","K/s",2,dt_hdiff(:,2))
-            call wstats(ngrid,"dt_ekman1","dt_ekman1","K/s",2,dt_ekman(:,1))
-            call wstats(ngrid,"dt_ekman2","dt_ekman2","K/s",2,dt_ekman(:,2))
+            call wstats(ngrid,"dt_hdiff1","dt_hdiff1","W m-2",2,dt_hdiff(:,1))
+            call wstats(ngrid,"dt_hdiff2","dt_hdiff2","W m-2",2,dt_hdiff(:,2))
+            call wstats(ngrid,"dt_ekman1","dt_ekman1","W m-2",2,dt_ekman(:,1))
+            call wstats(ngrid,"dt_ekman2","dt_ekman2","W m-2",2,dt_ekman(:,2))
+            call wstats(ngrid,"dt_gm1","dt_gm1","W m-2",2,dt_gm(:,1))
+            call wstats(ngrid,"dt_gm2","dt_gm2","W m-2",2,dt_gm(:,2))
             call wstats(ngrid,"tslab1","tslab1","K",2,tslab(:,1))
             call wstats(ngrid,"tslab2","tslab2","K",2,tslab(:,2))
@@ -2529 +2532 @@
          call writediagfi(ngrid,"tslab1","tslab1","K",2,tslab(:,1))
          call writediagfi(ngrid,"tslab2","tslab2","K",2,tslab(:,2))
-         call writediagfi(ngrid,"dt_hdiff1","dt_hdiff1","K/s",2,dt_hdiff(:,1))
-         call writediagfi(ngrid,"dt_hdiff2","dt_hdiff2","K/s",2,dt_hdiff(:,2))
-         call writediagfi(ngrid,"dt_ekman1","dt_ekman1","K/s",2,dt_ekman(:,1))
-         call writediagfi(ngrid,"dt_ekman2","dt_ekman2","K/s",2,dt_ekman(:,2))
+         call writediagfi(ngrid,"dt_hdiff1","dt_hdiff1","W m-2",2,dt_hdiff(:,1))
+         call writediagfi(ngrid,"dt_hdiff2","dt_hdiff2","W m-2",2,dt_hdiff(:,2))
+         call writediagfi(ngrid,"dt_ekman1","dt_ekman1","W m-2",2,dt_ekman(:,1))
+         call writediagfi(ngrid,"dt_ekman2","dt_ekman2","W m-2",2,dt_ekman(:,2))
+         call writediagfi(ngrid,"dt_gm1","dt_gm1","W m-2",2,dt_gm(:,1))
+         call writediagfi(ngrid,"dt_gm2","dt_gm2","W m-2",2,dt_gm(:,2))
          call writediagfi(ngrid,"rnat","nature of the surface","",2,rnat)
          call writediagfi(ngrid,"sensibFlux","sensible heat flux","w.m^-2",2,sensibFlux)
@@ -2563 +2568 @@
       if(callrad.and.(.not.newtonian))then
 
-         !call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
+         call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
          !call writediagfi(ngrid,"ALB_1st","First Band Surface albedo"," ",2,albedo(:,1))
          call writediagfi(ngrid,"ISR","incoming stellar rad.","W m-2",2,fluxtop_dn)