Changeset 1516 for trunk/LMDZ.GENERIC/libf/phystd/soil_settings.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_settings.F (modified) (4 diffs)

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

@@ -2 +2 @@
 
 !      use netcdf
-      use comsoil_h, only: layer, mlayer, inertiedat, volcapa
+      use comsoil_h, only: layer, mlayer, inertiedat, volcapa,
+     &                     lay1, alpha
       use iostart, only: inquire_field_ndims, get_var, get_field,
      &                   inquire_field, inquire_dimension_length
@@ -70 +71 @@
       real,dimension(:),allocatable :: newval
 
-      real alpha,lay1 ! coefficients for building layers
+      real malpha,mlay1 ! coefficients for building layers
       real xmin,xmax ! to display min and max of a field
 
@@ -133 +134 @@
       ! default mlayer distribution, following a power law:
       !  mlayer(k)=lay1*alpha**(k-1/2)
-        lay1=2.e-4
-        alpha=2
         do iloop=0,nsoil-1
           mlayer(iloop)=lay1*(alpha**(iloop-0.5))
@@ -142 +141 @@
       ! Assuming layer distribution follows mid-layer law:
       ! layer(k)=lay1*alpha**(k-1)
-      lay1=sqrt(mlayer(0)*mlayer(1))
-      alpha=mlayer(1)/mlayer(0)
+      mlay1=sqrt(mlayer(0)*mlayer(1))
+      malpha=mlayer(1)/mlayer(0)
+      ! Check that these values are the same as those prescibed for mlayers
+      if ((abs(mlay1-lay1)/lay1).gt.1.e-6) then
+         write(*,*) "soil_settings error: mlay1=",mlay1
+         write(*,*) " does not match comsoil_h lay1=",lay1
+         stop
+      endif
+      if ((abs(malpha-alpha)/alpha).gt.1.e-6) then
+         write(*,*) "soil_settings error: malpha=",malpha
+         write(*,*) " does not match comsoil_h alpha=",alpha
+         stop
+      endif
       do iloop=1,nsoil
-        layer(iloop)=lay1*(alpha**(iloop-1))
+        layer(iloop)=mlay1*(malpha**(iloop-1))
       enddo
+
 ! 2. Volumetric heat capacity (note: it is declared in comsoil_h)
 ! ---------------------------