Changeset 6153 for LMDZ6/branches/LMDZISOdevSN/libf

Timestamp:

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

Author:

Sebastien Nguyen

Message:

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

Location:

LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso

Files:

• lmdz_lscp_old.F90 (modified) (2 diffs)
• phyaqua_mod.F90 (modified) (3 diffs)
• phyetat0_mod.F90 (modified) (3 diffs)
• phyredem.F90 (modified) (6 diffs)
• physiq_mod.F90 (modified) (3 diffs)

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

@@ -1822 +1822 @@
         ! condensation, et le condensat en un tag spécifique
         ! Attention: zqn,zxtn ne servent qu'au calcul de zcond,zxtcond.
-        ! C'est en fait zxt qui se fait retranché
+        ! C'est en fait zxt qui se fait retrancher
         ! correction le 31 mars 2010: c'est à qn qu'il faut retrancher
         ! le condensat, car la condensation LS est un processus sous
@@ -2209 +2209 @@
        IF (rneb(i,k).GT.0.0) THEN
         if (zcond(i).gt.ridicule**2) then
-            ! le 21 dec: on change .gt.0 en .gt.ridicule**2 pour éviter valeur
-            ! ridiculement petites
+           ! le 21 dec: on change .gt.0 en .gt.ridicule**2 pour éviter valeur ridiculement petites
+           !cCAa condensed water in cloud ("zoliq") conserves its isotopic composition
            do ixt=1,ntraciso
              zxtoliq(ixt,i)=zoliq(i)*(zxtcond(ixt,i)/zcond(i))

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

@@ -52 +52 @@
     IMPLICIT NONE
 
-
-
     INTEGER, INTENT (IN) :: nlon, year_len, iflag_phys
     ! IM ajout latfi, lonfi
@@ -69 +67 @@
 #ifdef ISO
     REAL :: xtrun_off_lic_0(niso,nlon)
-    REAL :: xtsolsrf(niso,nlon, nbsrf), xtsnsrf(niso,nlon, nbsrf)
+    REAL :: xtsolsrf(niso,nlon, nbsrf)
+    REAL :: xtsnsrf(niso,nlon, nbsrf)
+    REAL :: xtsurfsnsrf(niso, nlon, nbsrf)
     REAL :: Rland_ice(niso,nlon)
     REAL :: Rsol(niso,nlon)
@@ -356 +356 @@
 !GG
 #ifdef ISO
-    CALL pbl_surface_init_iso(xtsnsrf,Rland_ice,Rsol)
+    CALL pbl_surface_init_iso(xtsnsrf, xtsurfsnsrf, Rland_ice, Rsol)
 #endif
 

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

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

@@ -91 +91 @@
 #ifdef ISO
   REAL xtsnow(niso,klon, nbsrf)
+  REAL xtsurfsnow(niso,klon, nbsrf)
   REAL xtrun_off_lic_0(niso,klon)
   REAL Rland_ice(niso,klon)
@@ -114 +115 @@
      CALL pbl_surface_final(fder, snow, qsurf, tsoil &
 #ifdef ISO
-       , xtsnow, Rland_ice, Rsol &
+       , xtsnow, xtsurfsnow, Rland_ice, Rsol &
 #endif       
        )
@@ -486 +487 @@
     WRITE(*,*) 'phyredem 486'
     CALL phyisoredem (pass, &
-   &           xtsnow, &
+   &           xtsnow, xtsurfsnow, &
    &           xtrun_off_lic_0, Rland_ice, Rsol, &
    &           run_off_lic_0)
@@ -562 +563 @@
 
 SUBROUTINE phyisoredem (pass, &
-     &          xtsnow, &
-     &          xtrun_off_lic_0, Rland_ice, Rsol, &
-     &          run_off_lic_0) 
- USE dimphy
+     &    xtsnow, xtsurfsnow, &
+     &    xtrun_off_lic_0, Rland_ice, Rsol, &
+     &    run_off_lic_0) 
+  USE dimphy
  !USE mod_grid_phy_lmdz
  !USE mod_phys_lmdz_para
@@ -592 +593 @@
   ! inputs
   !REAL xtsol(niso,klon)
-  REAL, INTENT(IN)    :: xtsnow(niso,klon,nbsrf) 
+  REAL, INTENT(IN)    :: xtsnow(niso,klon,nbsrf)
+  REAL, INTENT(IN)    :: xtsurfsnow(niso,klon,nbsrf)
   !REAL xtevap(ntiso,klon,nbsrf)     
   REAL, INTENT(IN)    :: xtrun_off_lic_0(niso,klon)
@@ -688 +690 @@
       iso_tmp_lonsrf(:,:)=xtsnow(ixt,:,:)
       CALL put_field_srf1(pass, "XTSNOW"//TRIM(outiso), "NEIGE",       iso_tmp_lonsrf)
+      
+      iso_tmp_lonsrf(:,:)=xtsurfsnow(ixt,:,:)
+      CALL put_field_srf1(pass, "XTSURFSNOW"//TRIM(outiso), "SURFACE SNOW", iso_tmp_lonsrf)
 
       iso_tmp(:)=xtsol(ixt,:)

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

@@ -1491 +1491 @@
        WRITE(lunout,*) '================================================='
        CALL phys_local_var_init
-       WRITE(lunout,*) 'MAXVAL(runoff_diag(:))=',MAXVAL(runoff_diag)
+!SN a enlever       WRITE(lunout,*) 'MAXVAL(runoff_diag(:))=',MAXVAL(runoff_diag)
        !
        !     appel a la lecture du run.def physique
@@ -1811 +1811 @@
 
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-       !! Un petit travail \`a faire ici.
+       !! Un petit travail a faire ici.
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
@@ -3690 +3690 @@
 !!SN    ENDIF
 !!SN    zx_ajustq = .FALSE.
-!!SN    IF (iflag_con.EQ.2) zx_ajustq=.TRUE.
+!!SN    IF (iflag_con.EQ.2) zx_ajustq=.TRUE. !SN ISOTOPES PAS PREVUS !
 !!SN    IF (zx_ajustq) THEN
 !!SN       DO i = 1, klon