Changeset 5629 for LMDZ6

Timestamp:

Apr 25, 2025, 4:25:24 PM (6 weeks ago)

Author:

Laurent Fairhead

Message:

Adaptation of heterogeneous continental surfaces to isotopes

Location:

LMDZ6/trunk/libf

Files:

• phylmd/surf_land_bucket_hetero_mod.F90 (modified) (5 diffs)
• phylmdiso/physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/surf_land_bucket_hetero_mod.F90

@@ -42 +42 @@
     USE cdrag_mod
 
+#ifdef ISO
+    use infotrac_phy, ONLY: niso
+#endif
+
 !****************************************************************************************
 ! Bucket calculations for surface. 
@@ -90 +94 @@
     REAL, DIMENSION(klon, nbtersrf), INTENT(OUT) :: fluxsens_tersrf, fluxlat_tersrf
 
+#ifdef ISO
+!    REAL, DIMENSION(ntiso,klon), INTENT(OUT) :: xtevap
+!    REAL, DIMENSION(klon),       INTENT(OUT) :: h1
+!    REAL, DIMENSION(niso,klon),  INTENT(OUT) :: xtrunoff_diag
+    REAL, DIMENSION(klon)  :: runoff_diag
+!    REAL, DIMENSION(niso,klon),  INTENT(IN)  :: Rland_ice
+#endif
 ! Local variables
 !****************************************************************************************
@@ -100 +111 @@
     REAL, DIMENSION(klon) :: dummy_riverflow,dummy_coastalflow 
     INTEGER               :: i, j, k
+#ifdef ISO
+    INTEGER               :: ixt
+    REAL, PARAMETER       :: t_coup = 273.15
+    REAL, DIMENSION(klon) :: fq_fonte_diag
+    REAL, DIMENSION(klon) :: fqfonte_diag
+    REAL, DIMENSION(klon) :: snow_evap_diag 
+    REAL, DIMENSION(klon) :: fqcalving_diag 
+    REAL                  :: max_eau_sol_diag  
+    REAL, DIMENSION(klon) :: run_off_lic_diag 
+    REAL :: coeff_rel_diag
+#endif
     !
     REAL, DIMENSION(klon)           :: zlev, geop1, speed, pblh, ri_in, sst
@@ -203 +225 @@
       
       !* Calculate snow height, run_off, age of snow
+#ifdef ISO
+    DO i=1,knon
+      ! initialisation:
+      fqfonte_diag(i)  =0.0
+      fq_fonte_diag(i) =0.0
+      snow_evap_diag(i)=0.0
+    ENDDO !DO i=1,knon 
+#endif    
 
       CALL fonte_neige( knon, is_ter, knindex, dtime, &
            tsurf, precip_rain, precip_snow, &
-           snow, qsol, tsurf_new, evap, icesub)
+           snow, qsol, tsurf_new, evap, icesub &
+#ifdef ISO    
+     & ,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag   &
+     & ,max_eau_sol_diag,runoff_diag,run_off_lic_diag,coeff_rel_diag   &
+#endif
+     &  )
     
       ! calculate the age of snow
@@ -325 +360 @@
         CALL fonte_neige( knon, is_ter, knindex, dtime, &
              tsurf_tersrf(:,j), precip_rain, precip_snow, &
-             snow, qsol, tsurf_new_tersrf(:,j), evap_tersrf(:,j), icesub)
+             snow, qsol, tsurf_new_tersrf(:,j), evap_tersrf(:,j), icesub &
+#ifdef ISO    
+     & ,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag   &
+     & ,max_eau_sol_diag,runoff_diag,run_off_lic_diag,coeff_rel_diag   &
+#endif
+     &  )
 
       ENDDO ! loop on sub-surfaces

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -3342 +3342 @@
                                 !>nrlmd+jyg
              treedrg &            
+!AM 
+            , tsurf_tersrf, tsoil_tersrf, qsurf_tersrf, tsurf_new_tersrf, &
+            cdragm_tersrf, cdragh_tersrf, &
+            swnet_tersrf, lwnet_tersrf, fluxsens_tersrf, fluxlat_tersrf &
 #ifdef ISO
      &   ,xtrain_fall, xtsnow_fall,xt_seri, &