Changeset 6153 for LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyetat0_mod.F90

Timestamp:

Apr 3, 2026, 11:21:47 AM (28 hours ago)

Author:

Sebastien Nguyen

Message:

implementation of isotopes in snow according to C.Agosta N.Dutrievoz and S.Sibin

File:

• LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyetat0_mod.F90 (modified) (3 diffs)

LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyetat0_mod.F90

@@ -122 +122 @@
   REAL :: missing_val
 #ifdef ISO
-  REAL xtsnow(niso,klon, nbsrf)
+  REAL xtsnow(niso,klon,nbsrf)
+  REAL xtsurfsnow(niso, klon,nbsrf)
   REAL xtrun_off_lic_0(niso,klon)
   REAL Rland_ice(niso,klon)
@@ -649 +650 @@
      ! initialise les isotopes       
   WRITE(*,*) 'phyetat0 1069' 
-  CALL phyisoetat0 (snow,run_off_lic_0, &
-  &           xtsnow,xtrun_off_lic_0, &
-  &           Rland_ice,Rsol)
+  CALL phyisoetat0(snow,run_off_lic_0, xtsnow, xtsurfsnow, xtrun_off_lic_0, Rland_ice, Rsol)
 #ifdef ISOVERIF 
   WRITE(*,*) 'phyetat0 1074'
   IF (iso_eau.GT.0) THEN
-    CALL iso_verif_egalite_vect2D(  &
-    &           xtsnow,snow, &
-    &           'phyetat0 1101a',niso,klon,nbsrf)
+    CALL iso_verif_egalite_vect2D(xtsnow, snow, 'phyetat0 1101a',niso,klon,nbsrf)
     DO i=1,klon  
-      CALL iso_verif_egalite(Rland_ice(iso_eau,i),1.0, &
-      &         'phyetat0 1101b')
+      CALL iso_verif_egalite(Rland_ice(iso_eau,i),1.0, 'phyetat0 1101b')
     ENDDO
   ENDIF
@@ -806 +802 @@
   !GG
 #ifdef ISO
-  CALL pbl_surface_init_iso(xtsnow,Rland_ice,Rsol)
+  CALL pbl_surface_init_iso(xtsnow, xtsurfsnow, Rland_ice, Rsol)
 #endif
   endif