Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM (revision 2018) @@ -106,4 +106,21 @@ state real DT_VDF ikj misc 1 - - "DT_VDF" "PBL CLMAIN TEMPERATURE TENDENCY" "K/s" state real DT_AJS ikj misc 1 - - "DT_AJS" "AJUSTEMENT SEC TEMPERATURE TENDENCY" "K/s" +state real CLOUDFRAC ikj misc 1 - - "CLOUDFRAC" "CLOUD FRACTION" "" +state real TOTCLOUDFRAC ij misc 1 - - "TOTCLOUDFRAC" "COLUMN CLOUD FRACTION" "" +state real GRAIN ikj misc 1 - - "RAIN" "RAIN kg/kg.s-1" "" +state real GSNOW ikj misc 1 - - "SNOW" "SNOW kg/kg.s-1" "" +state real ALBEQ ij misc 1 - - "ALBEQ" "ALBEDO EQUIVALENT" "" +state real FLUXTOP_DN ij misc 1 - - "FLUXTOP_DN" "Incoming SW radiation at the TOA" "W.m-2" +state real FLUXABS_SW ij misc 1 - - "FLUXABS_SW" "Absorbed SW flux" "W.m-2" +state real FLUXTOP_LW ij misc 1 - - "FLUXTOP_LW" "Outgoing LW flux TOA" "W.m-2" +state real FLUXSURF_SW ij misc 1 - - "FLUXSURF_SW" "SW Surface flux" "W.m-2" +state real FLUXSURF_LW ij misc 1 - - "FLUXSURF_LW" "LW Surface flux" "W.m-2" +state real FLXGRD ij misc 1 - - "FLXGRD" "heat flux from ground" "W.m-2" +state real LSCEZ ikj misc 1 - - "LSCEZ" "heat from largescale" "W.m-2" +state real H2OICE_REFF ikj misc 1 - - "H2OICE_REFF" "Radii of H2O ice" "m" +state real LATENT_HF ij misc 1 - - "LATENT_HF" "latent heat flux" "W.m-2" +state real REEVAP ij misc 1 - - "REEVAP" "re-evaporation flux of precipitation (integrated over the atmospheric column)" +state real SURFRAIN ij misc 1 - - "SURFRAIN" "Rainfall" "kg m-2 s-1" + #### #### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys (revision 2018) @@ -102,4 +102,21 @@ state real ZMAX ij misc 1 - - "ZMAX" "MAXIMUM LEVEL THERMALS" "m" state real HFMAX ij misc 1 - - "HFMAX" "MAXIMUM PBL HEAT FLUX" "m.K/s" +#### GENERIC OUTPUTS +state real CLOUDFRAC ikj misc 1 - - "CLOUDFRAC" "CLOUD FRACTION" "" +state real TOTCLOUDFRAC ij misc 1 - - "TOTCLOUDFRAC" "COLUMN CLOUD FRACTION" "" +state real GRAIN ikj misc 1 - - "RAIN" "RAIN kg/kg.s-1" "" +state real GSNOW ikj misc 1 - - "SNOW" "SNOW kg/kg.s-1" "" +state real ALBEQ ij misc 1 - - "ALBEQ" "ALBEDO EQUIVALENT" "" +state real FLUXTOP_DN ij misc 1 - - "FLUXTOP_DN" "Incoming SW radiation at the TOA" "W.m-2" +state real FLUXABS_SW ij misc 1 - - "FLUXABS_SW" "Absorbed SW flux" "W.m-2" +state real FLUXTOP_LW ij misc 1 - - "FLUXTOP_LW" "Outgoing LW flux TOA" "W.m-2" +state real FLUXSURF_SW ij misc 1 - - "FLUXSURF_SW" "SW Surface flux" "W.m-2" +state real FLUXSURF_LW ij misc 1 - - "FLUXSURF_LW" "LW Surface flux" "W.m-2" +state real FLXGRD ij misc 1 - - "FLXGRD" "heat flux from ground" "W.m-2" +state real LSCEZ ikj misc 1 - - "LSCEZ" "heat from largescale" "W.m-2" +state real H2OICE_REFF ikj misc 1 - - "H2OICE_REFF" "Radii of H2O ice" "m" +state real LATENT_HF ij misc 1 - - "LATENT_HF" "Latent heat flux" "W.m-2" +state real REEVAP ij misc 1 - - "REEVAP" "Re-evaporation flux of precipitation (integrated over the atmospheric column)" +state real SURFRAIN ij misc 1 - - "SURFRAIN" "Rainfall" "kg m-2 s-1" #### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_first_rk_step_part1.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_first_rk_step_part1.F (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_first_rk_step_part1.F (revision 2018) @@ -655,4 +655,21 @@ & ,DT_VDF=grid%DT_VDF & & ,DT_AJS=grid%DT_AJS & + & ,CLOUDFRAC=grid%CLOUDFRAC & + & ,TOTCLOUDFRAC=grid%TOTCLOUDFRAC & + & ,GRAIN=grid%GRAIN & + & ,GSNOW=grid%GSNOW & + & ,REEVAP=grid%REEVAP & + & ,SURFRAIN=grid%SURFRAIN & + & ,ALBEQ=grid%ALBEQ & + & ,FLUXTOP_DN=grid%FLUXTOP_DN & + & ,FLUXABS_SW=grid%FLUXABS_SW & + & ,FLUXTOP_LW=grid%FLUXTOP_LW & + & ,FLUXSURF_SW=grid%FLUXSURF_SW & + & ,FLUXSURF_LW=grid%FLUXSURF_LW & + & ,FLXGRD=grid%FLXGRD & + & ,LSCEZ=grid%LSCEZ & + & ,H2OICE_REFF=grid%H2OICE_REFF & + & ,LATENT_HF=grid%LATENT_HF & + & ,SWDOWNZ=grid%SWDOWNZ & & ,SWDOWNZ=grid%SWDOWNZ & & ,TAU_DUST=grid%TAU_DUST & Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM (revision 2018) @@ -117,4 +117,20 @@ state real DT_VDF ikj misc 1 - - "DT_VDF" "PBL CLMAIN TEMPERATURE TENDENCY" "K/s" state real DT_AJS ikj misc 1 - - "DT_AJS" "AJUSTEMENT SEC TEMPERATURE TENDENCY" "K/s" +state real CLOUDFRAC ikj misc 1 - - "CLOUDFRAC" "CLOUD FRACTION" "" +state real TOTCLOUDFRAC ij misc 1 - - "TOTCLOUDFRAC" "COLUMN CLOUD FRACTION" "" +state real GRAIN ikj misc 1 - - "RAIN" "RAIN kg/kg.s-1" "" +state real GSNOW ikj misc 1 - - "SNOW" "SNOW kg/kg.s-1" "" +state real ALBEQ ij misc 1 - - "ALBEQ" "ALBEDO EQUIVALENT" "" +state real FLUXTOP_DN ij misc 1 - - "FLUXTOP_DN" "Incoming SW radiation at the TOA" "W.m-2" +state real FLUXABS_SW ij misc 1 - - "FLUXABS_SW" "Absorbed SW flux" "W.m-2" +state real FLUXTOP_LW ij misc 1 - - "FLUXTOP_LW" "Outgoing LW flux TOA" "W.m-2" +state real FLUXSURF_SW ij misc 1 - - "FLUXSURF_SW" "SW Surface flux" "W.m-2" +state real FLUXSURF_LW ij misc 1 - - "FLUXSURF_LW" "LW Surface flux" "W.m-2" +state real FLXGRD ij misc 1 - - "FLXGRD" "heat flux from ground" "W.m-2" +state real LSCEZ ikj misc 1 - - "LSCEZ" "heat from largescale" "W.m-2" +state real H2OICE_REFF ikj misc 1 - - "H2OICE_REFF" "Radii of H2O ice" "m" +state real LATENT_HF ij misc 1 - - "LATENT_HF" "latent heat flux" "W.m-2" +state real REEVAP ij misc 1 - - "REEVAP" "re-evaporation flux of precipitation (integrated over the atmospheric column)" +state real SURFRAIN ij misc 1 - - "SURFRAIN" "Rainfall" "kg m-2 s-1" #### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys (revision 2018) @@ -99,4 +99,22 @@ state real ZMAX ij misc 1 - h "ZMAX" "MAXIMUM LEVEL THERMALS" "m" state real HFMAX ij misc 1 - h "HFMAX" "MAXIMUM PBL HEAT FLUX" "m.K/s" +#### GENERIC OUTPUTS +state real CLOUDFRAC ikj misc 1 - - "CLOUDFRAC" "CLOUD FRACTION" "" +state real TOTCLOUDFRAC ij misc 1 - - "TOTCLOUDFRAC" "COLUMN CLOUD FRACTION" "" +state real GRAIN ikj misc 1 - - "RAIN" "RAIN kg/kg.s-1" "" +state real GSNOW ikj misc 1 - - "SNOW" "SNOW kg/kg.s-1" "" +state real ALBEQ ij misc 1 - - "ALBEQ" "ALBEDO EQUIVALENT" "" +state real FLUXTOP_DN ij misc 1 - - "FLUXTOP_DN" "Incoming SW radiation at the TOA" "W.m-2" +state real FLUXABS_SW ij misc 1 - - "FLUXABS_SW" "Absorbed SW flux" "W.m-2" +state real FLUXTOP_LW ij misc 1 - - "FLUXTOP_LW" "Outgoing LW flux TOA" "W.m-2" +state real FLUXSURF_SW ij misc 1 - - "FLUXSURF_SW" "SW Surface flux" "W.m-2" +state real FLUXSURF_LW ij misc 1 - - "FLUXSURF_LW" "LW Surface flux" "W.m-2" +state real FLXGRD ij misc 1 - - "FLXGRD" "heat flux from ground" "W.m-2" +state real LSCEZ ikj misc 1 - - "LSCEZ" "heat from largescale" "W.m-2" +state real H2OICE_REFF ikj misc 1 - - "H2OICE_REFF" "Radii of H2O ice" "m" +state real LATENT_HF ij misc 1 - - "LATENT_HF" "Latent heat flux" "W.m-2" +state real REEVAP ij misc 1 - - "REEVAP" "Re-evaporation flux of precipitation (integrated over the atmospheric column)" +state real SURFRAIN ij misc 1 - - "SURFRAIN" "Rainfall" "kg m-2 s-1" + #### #### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F (revision 2018) @@ -897,4 +897,20 @@ & ,DT_VDF=grid%DT_VDF & & ,DT_AJS=grid%DT_AJS & + & ,CLOUDFRAC=grid%CLOUDFRAC & + & ,TOTCLOUDFRAC=grid%TOTCLOUDFRAC & + & ,GRAIN=grid%GRAIN & + & ,GSNOW=grid%GSNOW & + & ,REEVAP=grid%REEVAP & + & ,SURFRAIN=grid%SURFRAIN & + & ,ALBEQ=grid%ALBEQ & + & ,FLUXTOP_DN=grid%FLUXTOP_DN & + & ,FLUXABS_SW=grid%FLUXABS_SW & + & ,FLUXTOP_LW=grid%FLUXTOP_LW & + & ,FLUXSURF_SW=grid%FLUXSURF_SW & + & ,FLUXSURF_LW=grid%FLUXSURF_LW & + & ,FLXGRD=grid%FLXGRD & + & ,LSCEZ=grid%LSCEZ & + & ,H2OICE_REFF=grid%H2OICE_REFF & + & ,LATENT_HF=grid%LATENT_HF & & ,SWDOWNZ=grid%SWDOWNZ & & ,TAU_DUST=grid%TAU_DUST & Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_generic_lmd_new/update_outputs_physiq_mod.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_generic_lmd_new/update_outputs_physiq_mod.F (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_generic_lmd_new/update_outputs_physiq_mod.F (revision 2018) @@ -158,6 +158,10 @@ SWDOWNZ,TAU_DUST,QSURFDUST,& MTOT,ICETOT,TAU_ICE,& - HR_SW,HR_LW,HR_DYN,DT,DTRAD,DTVDF,DTAJS,& - RDUST,VMR_ICE,RICE) + HR_SW,HR_LW,HR_DYN,DT,DTRAD,& + RDUST,VMR_ICE,RICE,& + CLOUDFRAC,TOTCLOUDFRAC,& + RAIN,SNOW,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF) + USE comm_wrf !! to get fields to be written from physiq @@ -167,6 +171,9 @@ SWDOWNZ,TAU_DUST,QSURFDUST,& MTOT,ICETOT,TAU_ICE,& + TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: & - HR_SW,HR_LW,RDUST,VMR_ICE,RICE,HR_DYN,DT,DTRAD,DTVDF,DTAJS + HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DTRAD,RAIN,SNOW,& + LSCEZ,H2OICE_REFF INTEGER :: i,j,subs @@ -180,13 +187,22 @@ !! get diagnostics from physics - SWDOWNZ(i,j) = comm_SWDOWNZ(subs) - MTOT(i,j) = comm_MTOT(subs) HR_SW(i,kps:kpe,j) = comm_HR_SW(subs,kps:kpe) HR_LW(i,kps:kpe,j) = comm_HR_LW(subs,kps:kpe) - HR_DYN(i,kps:kpe,j) = comm_HR_DYN(subs,kps:kpe) - DT(i,kps:kpe,j) = comm_DT(subs,kps:kpe) - DTRAD(i,kps:kpe,j) = comm_DT_RAD(subs,kps:kpe) - DTVDF(i,kps:kpe,j) = comm_DT_VDF(subs,kps:kpe) - DTAJS(i,kps:kpe,j) = comm_DT_AJS(subs,kps:kpe) + CLOUDFRAC(i,kps:kpe,j)= comm_CLOUDFRAC(subs,kps:kpe) + TOTCLOUDFRAC(i,j)= comm_TOTCLOUDFRAC(subs) + RAIN(i,kps:kpe,j) = comm_RAIN(subs,kps:kpe) + SNOW(i,kps:kpe,j) = comm_SNOW(subs,kps:kpe) + ALBEQ(i,j)=comm_ALBEQ(subs) + FLUXTOP_DN(i,j)=comm_FLUXTOP_DN(subs) + FLUXABS_SW(i,j)=comm_FLUXABS_SW(subs) + FLUXTOP_LW(i,j)=comm_FLUXTOP_LW(subs) + FLUXSURF_SW(i,j)=comm_FLUXSURF_SW(subs) + FLUXSURF_LW(i,j)=comm_FLUXSURF_LW(subs) + FLXGRD(i,j)=comm_FLXGRD(subs) + LSCEZ(i,kps:kpe,j)=comm_LSCEZ(subs,kps:kpe) + H2OICE_REFF(i,kps:kpe,j)=comm_H2OICE_REFF(subs,kps:kpe) + LATENT_HF(i,j)=comm_LATENT_HF(subs) + REEVAP(i,j)=comm_REEVAP(subs) + SURFRAIN(i,j)=comm_SURFRAIN(subs) ENDDO Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_mars_lmd_new/update_outputs_physiq_mod.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_mars_lmd_new/update_outputs_physiq_mod.F (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_mars_lmd_new/update_outputs_physiq_mod.F (revision 2018) @@ -167,7 +167,10 @@ SWDOWNZ,TAU_DUST,QSURFDUST,& MTOT,ICETOT,TAU_ICE,& - HR_SW,HR_LW,HR_DYN,DT,DTRAD,DTVDF,DTAJS,& - RDUST,VMR_ICE,RICE) - + HR_SW,HR_LW,HR_DYN,DT,DTRAD,& + RDUST,VMR_ICE,RICE,& + CLOUDFRAC,TOTCLOUDFRAC,& + RAIN,SNOW,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF) + USE comm_wrf !! to get fields to be written from physiq @@ -176,8 +179,12 @@ REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: & SWDOWNZ,TAU_DUST,QSURFDUST,& - MTOT,ICETOT,TAU_ICE + MTOT,ICETOT,TAU_ICE,& + TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: & - HR_SW,HR_LW,HR_DYN,RDUST,VMR_ICE,RICE,DT,DTRAD,DTVDF,DTAJS - INTEGER :: i,j,subs + HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DTRAD,RAIN,SNOW,& + LSCEZ,H2OICE_REFF + INTEGER :: i,j,subs + DO j = jps,jpe Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_outputs_physiq_mod.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_outputs_physiq_mod.F (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_outputs_physiq_mod.F (revision 2018) @@ -71,6 +71,9 @@ SWDOWNZ,TAU_DUST,QSURFDUST,& MTOT,ICETOT,TAU_ICE,& - HR_SW,HR_LW,HR_DYN,DT,DTRAD,DTVDF,DTAJS,& - RDUST,VMR_ICE,RICE) + HR_SW,HR_LW,HR_DYN,DT,DTRAD,& + RDUST,VMR_ICE,RICE,& + CLOUDFRAC,TOTCLOUDFRAC,& + RAIN,SNOW,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF) INTEGER, INTENT(IN) :: ims,ime,jms,jme,kms,kme @@ -78,9 +81,11 @@ REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: & SWDOWNZ,TAU_DUST,QSURFDUST,& - MTOT,ICETOT,TAU_ICE + MTOT,ICETOT,TAU_ICE,& + TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: & - HR_SW,HR_LW,HR_DYN,RDUST,VMR_ICE,RICE,DT,DTRAD,DTVDF,DTAJS - INTEGER :: i,j,subs,k - REAL, DIMENSION(kms:kme) :: phr_sw,phr_lw,phr_dyn + HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DTRAD,RAIN,SNOW,& + LSCEZ,H2OICE_REFF + INTEGER :: i,j,subs open(17,file='prescribed_sw.txt',form='formatted',status='old') Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_outputs_physiq_mod.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_outputs_physiq_mod.F (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_outputs_physiq_mod.F (revision 2018) @@ -160,7 +160,10 @@ SWDOWNZ,TAU_DUST,QSURFDUST,& MTOT,ICETOT,TAU_ICE,& - HR_SW,HR_LW,HR_DYN,DT,DTRAD,DTVDF,DTAJS,& - RDUST,VMR_ICE,RICE) - + HR_SW,HR_LW,HR_DYN,DT,DTRAD,& + RDUST,VMR_ICE,RICE,& + CLOUDFRAC,TOTCLOUDFRAC,& + RAIN,SNOW,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF) + USE comm_wrf !! to get fields to be written from physiq @@ -169,8 +172,12 @@ REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: & SWDOWNZ,TAU_DUST,QSURFDUST,& - MTOT,ICETOT,TAU_ICE + MTOT,ICETOT,TAU_ICE,& + TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: & - HR_SW,HR_LW,HR_DYN,RDUST,VMR_ICE,RICE,DT,DTRAD,DTVDF,DTAJS - INTEGER :: i,j,subs + HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DTRAD,RAIN,SNOW,& + LSCEZ,H2OICE_REFF + INTEGER :: i,j,subs + DO j = jps,jpe Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new (revision 2017) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new (revision 2018) @@ -40,5 +40,8 @@ isfflx, diff_opt, km_opt, & HISTORY_INTERVAL, & - HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,DT_VDF,DT_AJS,SWDOWNZ,& + HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,DT_VDF,DT_AJS,& + CLOUDFRAC,TOTCLOUDFRAC, & + GRAIN,GSNOW,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF,SWDOWNZ,& TAU_DUST,RDUST,QSURFDUST,& MTOT,ICETOT,VMR_ICE,TAU_ICE,RICE,& @@ -104,5 +107,6 @@ MTOT,ICETOT,TAU_ICE,& HFMAX,ZMAX,& - USTM,HFX + USTM,HFX,TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN ! 3D arrays REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ) :: & @@ -110,5 +114,6 @@ REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(OUT ) :: & RTHBLTEN,RUBLTEN,RVBLTEN, & - HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,DT_VDF,DT_AJS,RDUST,VMR_ICE,RICE + HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,DT_VDF,DT_AJS,RDUST,VMR_ICE,RICE,& + CLOUDFRAC,GRAIN,GSNOW,LSCEZ,H2OICE_REFF REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: & M_Q2 @@ -734,6 +739,10 @@ SWDOWNZ,TAU_DUST,QSURFDUST,& MTOT,ICETOT,TAU_ICE,& - HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,DT_VDF,DT_AJS,& - RDUST,VMR_ICE,RICE) + HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,& + RDUST,VMR_ICE,RICE,& + CLOUDFRAC,TOTCLOUDFRAC,& + GRAIN,GSNOW,REEVAP,SURFRAIN,& + ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& + FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF) !!! print *, '** ',planet_type,'** OUTPUT PHYSICS DONE'