Changeset 1516 for trunk/LMDZ.GENERIC/libf/phystd/soil.F

Timestamp:

Mar 17, 2016, 9:45:36 AM (9 years ago)

Author:

emillour

Message:

Generic model:

• Added tests to ensure that soil model settings are adequate to resolve sub-surface diurnal and annual thermal waves.

EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/soil.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/soil.F

@@ -4 +4 @@
      &          capcal,fluxgrd)
 
-      use comsoil_h, only: layer, mlayer, volcapa
+      use comsoil_h, only: layer, mlayer, volcapa, inertiedat
+      use comcstfi_mod, only: pi, daysec
+      use planete_mod, only: year_day
+      use comgeomfi_h, only: long, lati
 
       implicit none
@@ -49 +52 @@
 ! local variables:
       integer ig,ik
-
+      real :: inertia_min,inertia_max
+      real :: diurnal_skin ! diurnal skin depth (m)
+      real :: annual_skin ! anuual skin depth (m)
 
 ! 0. Initialisations and preprocessing step
@@ -118 +123 @@
       !write(*,*)'soil: ig=',ig,' capcal(ig)=',capcal(ig)
       enddo ! of do ig=1,ngrid
-            
+      
+      ! Additional checks: is the vertical discretization sufficient
+      ! to resolve diurnal and annual waves?
+      do ig=1,ngrid
+        ! extreme inertia for this column
+        inertia_min=minval(inertiedat(ig,:))
+        inertia_max=maxval(inertiedat(ig,:))
+        ! diurnal and annual skin depth
+        diurnal_skin=(inertia_min/volcapa)*sqrt(daysec/pi)
+        annual_skin=(inertia_max/volcapa)*sqrt(year_day*daysec/pi)
+        if (0.5*diurnal_skin<layer(1)) then
+        ! one should have the fist layer be at least half of diurnal skin depth
+          write(*,*) "soil Error: grid point ig=",ig
+          write(*,*) "            longitude=",long(ig)*(180./pi)
+          write(*,*) "             latitude=",lati(ig)*(180./pi)
+          write(*,*) "  first soil layer depth ",layer(1)
+          write(*,*) "  not small enough for a diurnal skin depth of ",
+     &                diurnal_skin
+          write(*,*) " change soil layer distribution (comsoil_h.F90)"
+          stop
+        endif
+        if (2.*annual_skin>layer(nsoil)) then
+        ! one should have the full soil be at least twice the diurnal skin depth
+          write(*,*) "soil Error: grid point ig=",ig
+          write(*,*) "            longitude=",long(ig)*(180./pi)
+          write(*,*) "             latitude=",lati(ig)*(180./pi)
+          write(*,*) "  total soil layer depth ",layer(nsoil)
+          write(*,*) "  not large enough for an annual skin depth of ",
+     &                annual_skin
+          write(*,*) " change soil layer distribution (comsoil_h.F90)"
+          stop
+        endif
+      enddo ! of do ig=1,ngrid
+      
       else ! of if (firstcall)