Index: trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/iniphysiq_mod.F =================================================================== --- trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/iniphysiq_mod.F (revision 4051) +++ trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/iniphysiq_mod.F (revision 4123) @@ -29,5 +29,5 @@ use time_phylmdz_mod, only: dtphys, daysec,day_ini use planete_mod, only: year_day, periastr, apoastr, peri_day,& - obliquit, z0, lmixmin, emin_turb + obliquit, z0!, lmixmin, emin_turb use surfdat_h, only: emissiv,iceradius, & emisice,dtemisice Index: trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_inputs_physiq_mod.F =================================================================== --- trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_inputs_physiq_mod.F (revision 4051) +++ trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_inputs_physiq_mod.F (revision 4123) @@ -100,5 +100,5 @@ use planete_mod, only: year_day, periastr, apoastr, peri_day,& - obliquit, z0, lmixmin, emin_turb + obliquit, z0!, lmixmin, emin_turb use surfdat_h, only: emissiv,iceradius, & emisice,dtemisice Index: trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F =================================================================== --- trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F (revision 4051) +++ trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F (revision 4123) @@ -48,6 +48,6 @@ HFMAX,ZMAX,& USTM,HFX,& - SLPX,SLPY,RESTART,& - DT_COND,RH_CH4) + SLPX,SLPY,RESTART)!,& + !DT_COND,RH_CH4) ! NB: module_lmd_driver_output1.inc : output arguments generated from Registry @@ -115,6 +115,6 @@ RTHPLATEN,RUPLATEN,RVPLATEN, & HR_SW,HR_LW,HR_DYN,DT_RAD,& - CLOUDFRAC,RH,DQICE,DQVAP,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,& - DT_COND, RH_CH4 + CLOUDFRAC,RH,DQICE,DQVAP,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF!,& + !DT_COND, RH_CH4 REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: & P_Q2 @@ -447,14 +447,14 @@ q_prof(:,2) = moist(i,kps:kpe,j,P_QC) / (1.d0 + moist(i,kps:kpe,j,P_QV)) ! conversion from mass mixing ratio in WRF to specific concentration in Physiq -ELSE IF (TRACER_MODE == 61) THEN - ! to be clean we should have an automatized process that makes sure that moist is sent to igcm_h2o_vap and etc. - q_prof(:,1) = SCALAR(i,kps:kpe,j,P_mu_m0as) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,2) = SCALAR(i,kps:kpe,j,P_mu_m3as) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,3) = SCALAR(i,kps:kpe,j,P_mu_m0af) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,4) = SCALAR(i,kps:kpe,j,P_mu_m3af) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,5) = SCALAR(i,kps:kpe,j,P_mu_m0n) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,6) = SCALAR(i,kps:kpe,j,P_mu_m3n) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,7) = SCALAR(i,kps:kpe,j,P_mu_m3CH4) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) - q_prof(:,8) = SCALAR(i,kps:kpe,j,P_CH4) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +!ELSE IF (TRACER_MODE == 61) THEN +! ! to be clean we should have an automatized process that makes sure that moist is sent to igcm_h2o_vap and etc. +! q_prof(:,1) = SCALAR(i,kps:kpe,j,P_mu_m0as) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,2) = SCALAR(i,kps:kpe,j,P_mu_m3as) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,3) = SCALAR(i,kps:kpe,j,P_mu_m0af) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,4) = SCALAR(i,kps:kpe,j,P_mu_m3af) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,5) = SCALAR(i,kps:kpe,j,P_mu_m0n) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,6) = SCALAR(i,kps:kpe,j,P_mu_m3n) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,7) = SCALAR(i,kps:kpe,j,P_mu_m3CH4) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) +! q_prof(:,8) = SCALAR(i,kps:kpe,j,P_CH4) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4)) ELSE q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1) !! the names were set above !! one dummy tracer in WRF @@ -663,6 +663,6 @@ DQICE,DQVAP,REEVAP,SURFRAIN,& ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,& - FLUXSURF_LW,FLXGRD,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,LATENT_HF,& - DT_COND,RH_CH4) + FLUXSURF_LW,FLXGRD,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,LATENT_HF)!,& + !DT_COND,RH_CH4) !!! !print *, '** ',planet_type,'** OUTPUT PHYSICS DONE' @@ -752,21 +752,21 @@ SCALAR(i,kps:kpe,j,P_MARKER) = SCALAR(i,kps:kpe,j,P_MARKER)*exp(-dt/tau_decay) SCALAR(i,1,j,P_MARKER) = 1. !! this tracer is emitted in the surface layer - CASE(61) !emoisan - scalar(i,kps:kpe,j,P_mu_m0as)=scalar(i,kps:kpe,j,P_mu_m0as) & - +zdqfi_omp(subs,kps:kpe,1)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_mu_m3as)=scalar(i,kps:kpe,j,P_mu_m3as) & - +zdqfi_omp(subs,kps:kpe,2)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_mu_m0af)=scalar(i,kps:kpe,j,P_mu_m0af) & - +zdqfi_omp(subs,kps:kpe,3)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_mu_m3af)=scalar(i,kps:kpe,j,P_mu_m3af) & - +zdqfi_omp(subs,kps:kpe,4)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_mu_m0n)=scalar(i,kps:kpe,j,P_mu_m0n) & - +zdqfi_omp(subs,kps:kpe,5)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_mu_m3n)=scalar(i,kps:kpe,j,P_mu_m3n) & - +zdqfi_omp(subs,kps:kpe,6)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_mu_m3CH4)=scalar(i,kps:kpe,j,P_mu_m3CH4) & - +zdqfi_omp(subs,kps:kpe,7)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) - scalar(i,kps:kpe,j,P_CH4)=scalar(i,kps:kpe,j,P_CH4) & - +zdqfi_omp(subs,kps:kpe,8)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! CASE(61) !emoisan +! scalar(i,kps:kpe,j,P_mu_m0as)=scalar(i,kps:kpe,j,P_mu_m0as) & +! +zdqfi_omp(subs,kps:kpe,1)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_mu_m3as)=scalar(i,kps:kpe,j,P_mu_m3as) & +! +zdqfi_omp(subs,kps:kpe,2)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_mu_m0af)=scalar(i,kps:kpe,j,P_mu_m0af) & +! +zdqfi_omp(subs,kps:kpe,3)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_mu_m3af)=scalar(i,kps:kpe,j,P_mu_m3af) & +! +zdqfi_omp(subs,kps:kpe,4)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_mu_m0n)=scalar(i,kps:kpe,j,P_mu_m0n) & +! +zdqfi_omp(subs,kps:kpe,5)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_mu_m3n)=scalar(i,kps:kpe,j,P_mu_m3n) & +! +zdqfi_omp(subs,kps:kpe,6)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_mu_m3CH4)=scalar(i,kps:kpe,j,P_mu_m3CH4) & +! +zdqfi_omp(subs,kps:kpe,7)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) +! scalar(i,kps:kpe,j,P_CH4)=scalar(i,kps:kpe,j,P_CH4) & +! +zdqfi_omp(subs,kps:kpe,8)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4)) CASE DEFAULT !SCALAR(i,kps:kpe,j,2:nq+1)=SCALAR(i,kps:kpe,j,2:nq+1)+zdqfi_omp(subs,kps:kpe,1:nq)*dt !!! here dt is needed