Index: LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/lmdz_lscp_old.F90 =================================================================== --- LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/lmdz_lscp_old.F90 (revision 6151) +++ LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/lmdz_lscp_old.F90 (revision 6153) @@ -1822,5 +1822,5 @@ ! 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,6 +2209,6 @@ 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)) Index: LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyaqua_mod.F90 =================================================================== --- LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyaqua_mod.F90 (revision 6151) +++ LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyaqua_mod.F90 (revision 6153) @@ -52,6 +52,4 @@ IMPLICIT NONE - - INTEGER, INTENT (IN) :: nlon, year_len, iflag_phys ! IM ajout latfi, lonfi @@ -69,5 +67,7 @@ #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,5 +356,5 @@ !GG #ifdef ISO - CALL pbl_surface_init_iso(xtsnsrf,Rland_ice,Rsol) + CALL pbl_surface_init_iso(xtsnsrf, xtsurfsnsrf, Rland_ice, Rsol) #endif Index: LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyetat0_mod.F90 =================================================================== --- LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyetat0_mod.F90 (revision 6151) +++ LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyetat0_mod.F90 (revision 6153) @@ -122,5 +122,6 @@ 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,16 +650,11 @@ ! 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,5 +802,5 @@ !GG #ifdef ISO - CALL pbl_surface_init_iso(xtsnow,Rland_ice,Rsol) + CALL pbl_surface_init_iso(xtsnow, xtsurfsnow, Rland_ice, Rsol) #endif endif Index: LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyredem.F90 =================================================================== --- LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyredem.F90 (revision 6151) +++ LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/phyredem.F90 (revision 6153) @@ -91,4 +91,5 @@ #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,5 +115,5 @@ CALL pbl_surface_final(fder, snow, qsurf, tsoil & #ifdef ISO - , xtsnow, Rland_ice, Rsol & + , xtsnow, xtsurfsnow, Rland_ice, Rsol & #endif ) @@ -486,5 +487,5 @@ WRITE(*,*) 'phyredem 486' CALL phyisoredem (pass, & - & xtsnow, & + & xtsnow, xtsurfsnow, & & xtrun_off_lic_0, Rland_ice, Rsol, & & run_off_lic_0) @@ -562,8 +563,8 @@ 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,5 +593,6 @@ ! 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,4 +690,7 @@ 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,:) Index: LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/physiq_mod.F90 =================================================================== --- LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/physiq_mod.F90 (revision 6151) +++ LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/physiq_mod.F90 (revision 6153) @@ -1491,5 +1491,5 @@ 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,5 +1811,5 @@ ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! - !! Un petit travail \`a faire ici. + !! Un petit travail a faire ici. ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! @@ -3690,5 +3690,5 @@ !!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