Changeset 5627 for LMDZ6/trunk/libf/phylmd/dyn1d/scm.f90

Timestamp:

Apr 25, 2025, 2:21:02 PM (2 months ago)

Author:

amaison

Message:

Representation of heterogeneous continental subsurfaces with parameter or flux aggregation in the simplified surface model (bucket) for 1D case studies.

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/scm.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/scm.f90 ¶

@@ -8 +8 @@
        clwcon, detr_therm, &
        qsol, fevap, z0m, z0h, agesno, &
+       frac_tersrf, z0m_tersrf, ratio_z0m_z0h_tersrf, &
+       albedo_tersrf, beta_tersrf, inertie_tersrf, &
+       alpha_soil_tersrf, period_tersrf, hcond_tersrf, &
+       tsurfi_tersrf, tsoili_tersrf, tsoil_depth, &
        du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
        falb_dir, falb_dif, &
@@ -178 +182 @@
       real :: fder(1),snsrf(1,nbsrf),qsurfsrf(1,nbsrf)
       real :: tsoil(1,nsoilmx,nbsrf)
+      ! AM
+      REAL, ALLOCATABLE, DIMENSION(:,:) :: tsoil_ter_srf2 ! resized initial soil temperature on vertical levels (K) 
+      REAL, ALLOCATABLE, DIMENSION(:,:) :: tsoil_depths2  ! resized soil depth at which inititial temperature is given (m)
 
 !---------------------------------------------------------------------
@@ -222 +229 @@
 !                           <> 0, tendencies of forcing are not added
       INTEGER :: flag_inhib_forcing = 0
-
+      CHARACTER(len=80) :: abort_message
+      CHARACTER(len=20) :: modname = 'scm'
 
       print*,'VOUS ENTREZ DANS LE 1D FORMAT STANDARD'
@@ -385 +393 @@
 !     call init_phys_lmdz(1,1,llm,1,(/1/)) ! job now done via iniphysiq
 !     but we still need to initialize dimphy module (klon,klev,etc.)  here.
-      call init_dimphy1D(1,llm)
+      call init_dimphy1D(1,llm,nb_ter_srf,nb_tsoil_depths)
       call suphel
       call init_infotrac
@@ -560 +568 @@
         agesno  = xagesno
         tsoil(:,:,:)=tsurf
+
+        iflag_hetero_surf = 0
+        CALL getin('iflag_hetero_surf',iflag_hetero_surf)
+
+        IF (iflag_hetero_surf .GT. 0) THEN
+          PRINT*, 'scm iflag_hetero_surf', iflag_hetero_surf
+          IF ((nbtersrf .LT. 2) .OR. (nbtersrf .GT. max_nbtersrf)) THEN
+            abort_message='The number of continental sub-surfaces (nb_ter_srf) must be between 2 and 5'
+            CALL abort_physic(modname,abort_message,1)
+          ENDIF
+          ! resized initial soil temperature on vertical levels and soil depth at which inititial temperature is given
+          ALLOCATE(tsoil_ter_srf2(nbtsoildepths,nbtersrf))
+          ALLOCATE(tsoil_depths2(nbtsoildepths,nbtersrf))
+          tsoil_ter_srf2(:,:) = 0.
+          tsoil_depths2(:,:) = 0.
+          DO i=1, nbtersrf
+            DO l=1, nbtsoildepths
+              k = nbtsoildepths*(i-1)+l
+              tsoil_ter_srf2(l,i) = tsoil_ter_srf(k)
+              tsoil_depths2(l,i) = tsoil_depths(k)
+            ENDDO
+          ENDDO
+          !
+          DO i=1, nbtersrf
+            frac_tersrf(:,i) = frac_ter_srf(i)                   ! fraction of land surface heterogeneity (-)
+            z0m_tersrf(:,i) = rugos_ter_srf(i)                   ! roughness length for momentum of land sub-surfaces (m)
+            ratio_z0m_z0h_tersrf(:,i) = ratio_z0m_z0h_ter_srf(i) ! ratio of heat to momentum roughness length of land sub-surfaces (-)
+            albedo_tersrf(:,i) = albedo_ter_srf(i)               ! albedo of land sub-surfaces (-)
+            beta_tersrf(:,i) = beta_ter_srf(i)                   ! evapotranspiration coef of land sub-surfaces (-)
+            inertie_tersrf(:,i) = inertie_ter_srf(i)             ! soil thermal inertia of land sub-surfaces (J/m2/K/s1/2)
+            hcond_tersrf(:,i) = hcond_ter_srf(i)                 ! soil heat conductivity (W/(m.K))
+            tsurfi_tersrf(:,i) = tsurf_ter_srf(i)                ! initial surface temperature (K)
+            DO l=1, nbtsoildepths
+              tsoili_tersrf(:,l,i) = tsoil_ter_srf2(l,i)         ! initial soil temperature on vertical levels (K)
+              tsoil_depth(:,l,i) = tsoil_depths2(l,i)
+            ENDDO
+          ENDDO
+          alpha_soil_tersrf = alpha_soil_ter_srf               ! ratio between the thicknesses of 2 successive layers (-)
+          period_tersrf = period_ter_srf                       ! temperature oscillation amplitude period
+          !
+          DEALLOCATE(tsoil_ter_srf2)
+          DEALLOCATE(tsoil_depths2)
+        ENDIF
+
 !-----------------------------------------------------------------------
         call pbl_surface_init(fder, snsrf, qsurfsrf, tsoil)