Index: trunk/LMDZ.TITAN/libf/phytitan/comm_wrf.F90 =================================================================== --- trunk/LMDZ.TITAN/libf/phytitan/comm_wrf.F90 (revision 3660) +++ trunk/LMDZ.TITAN/libf/phytitan/comm_wrf.F90 (revision 3661) @@ -18,11 +18,14 @@ REAL,SAVE,ALLOCATABLE :: comm_FLUXSURF_LW(:) REAL,SAVE,ALLOCATABLE :: comm_FLXGRD(:) + REAL,SAVE,ALLOCATABLE :: comm_zqfi_omp(:,:,:) + REAL,SAVE,ALLOCATABLE :: comm_zdtlc(:,:) contains - subroutine allocate_comm_wrf(ngrid,nlayer) + subroutine allocate_comm_wrf(ngrid,nlayer,nq) implicit none - integer,intent(in) :: ngrid ! number of atmospheric columns + integer,intent(in) :: ngrid ! number of atmospheric columns integer,intent(in) :: nlayer ! number of atmospheric layers + integer,intent(in) :: nq ! number of tracers allocate(comm_HR_SW(ngrid,nlayer)) allocate(comm_HR_LW(ngrid,nlayer)) @@ -38,4 +41,6 @@ allocate(comm_FLUXSURF_LW(ngrid)) allocate(comm_FLXGRD(ngrid)) + allocate(comm_zqfi_omp(ngrid,nlayer,nq)) + allocate(comm_zdtlc(ngrid,nlayer)) end subroutine allocate_comm_wrf @@ -56,4 +61,6 @@ deallocate(comm_FLUXSURF_LW) deallocate(comm_FLXGRD) + deallocate(comm_zqfi_omp) + deallocate(comm_zdtlc) end subroutine deallocate_comm_wrf Index: trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 =================================================================== --- trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (revision 3660) +++ trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (revision 3661) @@ -46,7 +46,8 @@ #else use comm_wrf, only : comm_HR_SW, comm_HR_LW, & - comm_FLUXTOP_DN,comm_FLUXABS_SW,& - comm_FLUXTOP_LW,comm_FLUXSURF_SW,& - comm_FLUXSURF_LW,comm_FLXGRD + comm_FLUXTOP_DN, comm_FLUXABS_SW,& + comm_FLUXTOP_LW, comm_FLUXSURF_SW,& + comm_FLUXSURF_LW, comm_FLXGRD,& + comm_zqfi_omp, comm_zdtlc #endif #ifdef CPP_XIOS @@ -138,5 +139,5 @@ ! pdv(ngrid,nlayer) \ Temporal derivative of the corresponding ! pdt(ngrid,nlayer) / variables due to physical processes. -! pdq(ngrid,nlayer) / +! pdq(ngrid,nlayer,nq) / ! pdpsrf(ngrid) / ! @@ -331,5 +332,5 @@ character*2 :: str2 -#ifndef MESOSCALE +!#ifndef MESOSCALE ! Local variables for Titan chemistry and microphysics @@ -363,5 +364,23 @@ #endif - logical file_ok +#ifdef MESOSCALE + LOGICAL, SAVE :: moyzon_ch ! used for zonal averages in Titan + REAL, ALLOCATABLE :: zplevbar(:,:) + REAL, ALLOCATABLE :: zplaybar(:,:) + REAL, ALLOCATABLE :: ztfibar(:,:) + REAL, ALLOCATABLE :: zqfibar(:,:,:) + REAL, ALLOCATABLE :: zphibar(:,:) + REAL, ALLOCATABLE :: zphisbar(:) + REAL, ALLOCATABLE :: zzlevbar(:,:) + REAL, ALLOCATABLE :: zzlaybar(:,:) +! REAL, DIMENSION(:,:) :: zplevbar(:,:) +! REAL, DIMENSION(:,:) :: zplaybar(:,:) +! REAL, DIMENSION(:,:) :: ztfibar(:,:) +! REAL, DIMENSION(:,:,:) :: zqfibar(:,:,:) +! REAL, DIMENSION(:,:) :: zphibar(:,:) +! REAL, DIMENSION(:) :: zphisbar(:) +! REAL, DIMENSION(:,:) :: zzlevbar(:,:) +! REAL, DIMENSION(:,:) :: zzlaybar(:,:) +#endif !----------------------------------------------------------------------------- @@ -385,6 +404,8 @@ END INTERFACE -#endif - +!#endif + + logical file_ok + !================================================================================================== @@ -447,5 +468,5 @@ int_dtauv(:,:,:) = 0.D0 -#ifndef MESOSCALE +!#ifndef MESOSCALE IF (callmufi .AND. (.NOT. uncoupl_optic_haze)) THEN haze_opt_file=trim(datadir)//'/optical_tables/HAZE_OPTIC_'//trim(adjustl(tmp1))//'x'//trim(adjustl(tmp2))//'.DAT' @@ -459,12 +480,15 @@ endif ENDIF -#endif +!#endif endif + +#ifdef MESOSCALE + moyzon_ch = .false. !no zonal mean for mesoscale +#endif #ifndef MESOSCALE ! Initialize names and timestep for chemistry ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - if (callchim) then @@ -479,4 +503,5 @@ endif +#endif ! Initialize microphysics. @@ -491,5 +516,4 @@ ENDIF -#endif #ifdef CPP_XIOS @@ -717,5 +741,5 @@ ! JVO 19 : We shall always have correct altitudes in chemistry no matter what's in physics ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ -#ifndef MESOSCALE +!#ifndef MESOSCALE if (moyzon_ch) then ! Zonal averages @@ -748,5 +772,5 @@ else ! if not moyzon -#endif +!#endif DO ig=1,ngrid @@ -763,7 +787,7 @@ ENDDO -#ifndef MESOSCALE +!#ifndef MESOSCALE endif ! moyzon -#endif +!#endif ! ------------------------------------------------------------------------------------- @@ -1065,5 +1089,4 @@ if (tracer) then -#ifndef MESOSCALE ! ------------------- ! V.1 Microphysics @@ -1072,5 +1095,4 @@ ! We must call microphysics before chemistry, for condensation ! if (callmufi) then - zzlev(:,nlayer+1)=zzlay(:,nlayer)+(zzlay(:,nlayer)-zzlev(:,nlayer)) ! JVO 19 : We assume zzlev isn't reused later on (could be done cleaner) @@ -1148,4 +1170,5 @@ endif ! callmufi +#ifndef MESOSCALE ! ----------------- ! V.2. Chemistry @@ -1292,5 +1315,5 @@ endif ! end of 'callchim' -! END MESOSCALE +!! END ifndef MESOSCALE #endif @@ -2010,4 +2033,6 @@ comm_FLXGRD(1:ngrid)=fluxgrd(1:ngrid) sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid) + comm_zqfi_omp(1:ngrid,1:nlayer,1:nq) = zq(1:ngrid,1:nlayer,1:nq) + comm_zdtlc(1:ngrid,1:nlayer) = zdtlc(1:ngrid,1:nlayer) #endif Index: trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/callphysiq_mod.F =================================================================== --- trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/callphysiq_mod.F (revision 3661) +++ trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/callphysiq_mod.F (revision 3661) @@ -0,0 +1,72 @@ +! +! $Id: $ +! +MODULE callphysiq_mod + +IMPLICIT NONE + +CONTAINS + +SUBROUTINE call_physiq(planet_type,klon,llm,nqtot, & + debut_split,lafin_split) + + USE variables_mod + USE physiq_mod, ONLY: physiq + use phys_state_var_mod , only : qsurf + use comm_wrf, only : allocate_comm_wrf + use tracer_h, only: noms + + IMPLICIT NONE + character(len=10),INTENT(IN) :: planet_type ! planet type ('earth','mars',...) + INTEGER,INTENT(IN) :: klon ! (local) number of atmospheric columns + INTEGER,INTENT(IN) :: llm ! number of atmospheric layers + INTEGER,INTENT(IN) :: nqtot ! number of tracers + LOGICAL,INTENT(IN) :: debut_split ! .true. if very first call to physics + LOGICAL,INTENT(IN) :: lafin_split ! .true. if last call to physics + + ! Local variables +! CHARACTER(len=11) :: modname="call_physiq" + +! Sanity check on physics package type + IF (debut_split) THEN + IF (planet_type.ne."titan") THEN + PRINT*,"wrong planet_type for this physics package" + STOP + ENDIF + ENDIF + + !qsurf(:,:)=zqfi_omp(:,1,:) + !JL22 what is this, why in the world would we want to put the value of the last layer at the surface + ! these two quantities do not even have the same units (kg/kg and kg/m^2)!!! + ! seems it is the same for some other physics. It is probably a bug. + !EMo25 zqfi_omp in kg/kghumidair? + + call allocate_comm_wrf(klon,llm,nqtot) +! Call physics package with required inputs/outputs + CALL physiq(klon, & ! ngrid + llm, & ! nlayer + nqtot, & ! nq + noms, & ! nametrac + debut_split, & ! firstcall + lafin_split, & ! lastcall + jD_cur, & ! pday + jH_cur_split, & ! ptime + zdt_split, & ! ptimestep + zplev_omp, & ! pplev + zplay_omp, & ! pplay + zphi_omp, & ! pphi + zufi_omp, & ! pu + zvfi_omp, & ! pv + ztfi_omp, & ! pt + zqfi_omp, & ! pq + flxwfi_omp, & ! flxw + zdufi_omp, & ! pdu + zdvfi_omp, & ! pdv + zdtfi_omp, & ! pdt + zdqfi_omp, & ! pdq + zdpsrf_omp) ! pdpsrf + + +END SUBROUTINE call_physiq + +END MODULE callphysiq_mod Index: trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/iniphysiq_mod.F =================================================================== --- trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/iniphysiq_mod.F (revision 3661) +++ trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/iniphysiq_mod.F (revision 3661) @@ -0,0 +1,173 @@ +MODULE iniphysiq_mod + +CONTAINS + +subroutine iniphysiq(ngrid,nlayer,nq,piphysiq,& + punjours, pdayref, & + prad,pg,pr,pcpp,iflag_phys) + +!use control_mod, only: nday +!use surf_heat_transp_mod, only: ini_surf_heat_transp +!use infotrac, only : nqtot ! number of advected tracers +!USE comvert_mod, ONLY: ap,bp,preff +use inifis_mod, only: inifis +use ioipsl_getin_p_mod, only: getin_p + +!use inigeomphy_mod, only: inigeomphy +!use geometry_mod, only: cell_area, & ! physics grid area (m2) +! longitude, & ! longitudes (rad) +! latitude ! latitudes (rad) +! necessary to get klon_omp +!USE mod_phys_lmdz_para, ONLY: klon_omp ! number of columns (on local omp grid) +USE mod_phys_lmdz_para, ONLY: Init_phys_lmdz_para +USE dimphy, ONLY: init_dimphy +USE phys_state_var_mod +!use planete_mod, only: year_day, periastr, apoastr, peri_day,& +! obliquit, z0, lmixmin, emin_turb +! init_planete_mod +use planete_mod +use time_phylmdz_mod, only: dtphys, daysec,day_ini +use planete_mod, only: year_day, periastr, apoastr, peri_day,& + obliquit, z0, lmixmin, emin_turb +use surfdat_h, only: emissiv,iceradius, & + emisice,dtemisice +use comcstfi_mod, only: omeg,mugaz +use tracer_h, only: nqtot_p +!use comm_wrf, only : allocate_comm_wrf +USE variables_mod, only: zdt_split +use mod_grid_phy_lmdz, only: nbp_lev, klon_glo + +implicit none + +INCLUDE 'mpif.h' + +REAL,intent(in) :: prad +REAL,intent(in) :: pg +REAL,intent(in) :: pr +REAL,intent(in) :: pcpp +REAL,intent(in) :: punjours +!DOUBLE PRECISION,intent(in) :: ptimestep + +!real,intent(in) :: prad ! radius of the planet (m) +!real,intent(in) :: pg ! gravitational acceleration (m/s2) +!real,intent(in) :: pr ! ! reduced gas constant R/mu +!real,intent(in) :: pcpp ! specific heat Cp +!real,intent(in) :: punjours ! length (in s) of a standard day [daysec] +integer,intent(in) :: pdayref ! reference day of for the simulation [day_ini] +integer,intent(in) :: iflag_phys ! type of physics to be called + +integer :: nday=0 ! this is dummy for mesoscale (in dyn3d/control_mod) + +integer,intent(in) :: ngrid ! number of physics columns for this MPI process +integer,intent(in) :: nlayer ! number of atmospheric layers +integer,intent(in) :: nq ! number of tracers +!real,intent(in) :: phour_ini ! start time (fraction of day) of the run +!0= initialize physics distribution, global fields and geometry + ! (i.e. things in phy_common or dynphy_lonlat) + + ! the distinct part for all planetary physics (ie. things in phystd) + !------------------------------------------ + +nbp_lev = nlayer ! emoisan +klon_glo = 1 !emoisan !we run the physics as if we were in 1D +CALL Init_phys_lmdz_para(1,1,1,MPI_COMM_WORLD) + +!call phys_state_var_init +print*,'ngrid',ngrid,'nlayer',nlayer +call init_dimphy(ngrid,nlayer) +call phys_state_var_init(nqtot_p) +! copy over preff , ap() and bp() +!call ini_planete_mod(nlayer,prefff,apf,bpf) + +! for slab ocean, copy over some arrays +ok_slab_ocean=.false. ! default value +!call getin_p("ok_slab_ocean",ok_slab_ocean) +!if (ok_slab_ocean) then +! call ini_surf_heat_transp(ip1jm,ip1jmp1,unsairez,fext,unsaire, & +! cu,cuvsurcv,cv,cvusurcu,aire,apoln,apols, & +! aireu,airev) +!endif + +!dummy=1. !used before instead of zdt_split (EMoi) +lat(:)=0. +long(:)=0. +cellarea(:)=1. +print*,'pg',pg +!ppunjours=punjours +ppunjours=1. +pprad=prad +ppg=pg +ppr=pr +ppcpp=pcpp +call inifis(ngrid,nlayer,nq,pdayref,ppunjours,nday,zdt_split, & + lat,long,cellarea,pprad,ppg,ppr,ppcpp) + + open(17,file='controle.txt',form='formatted',status='old') + rewind(17) + read(17,*) + read(17,*) + read(17,*) day_ini !(tab0+3) + read(17,*) + read(17,*) !tab0+5) + read(17,*) !omeg !(tab0+6) + read(17,*) !(tab0+7) + read(17,*) !mugaz + read(17,*) !(tab0+9) + read(17,*) daysec + read(17,*) dtphys !tab0+11) + read(17,*) + read(17,*) + read(17,*) year_day !(tab0+14) + read(17,*) periastr !tab0+15) + read(17,*) apoastr !tab0+16) + read(17,*) peri_day !tab0+17) + read(17,*) obliquit !tab0+18) + read(17,*) z0 + read(17,*) + read(17,*) + read(17,*) + read(17,*) + read(17,*) emisice(1) + read(17,*) emisice(2) + read(17,*) emissiv + read(17,*) + read(17,*) + read(17,*) + read(17,*) + read(17,*) iceradius(1) + read(17,*) iceradius(2) + read(17,*) dtemisice(1) + read(17,*) dtemisice(2) + close(17) + !print*,'g',g + + !emissiv(:)=EMIS +! cloudfrac(:,:)=0.5 +! totcloudfrac(:)=0.5 +! hice(:)=0. +! rnat(:)=0. +! pctsrf_sic(:)=0. +! tsea_ice(:)=0. + !qsurf(:,:) = 0. + print*,'check' + print*,'iceradius',iceradius,'dtemisice',dtemisice + print*,'apoastr,periastr,year_day,peri_day,obliq',apoastr,periastr,year_day,peri_day,obliquit + print*,'emissiv',emissiv + print*,'mugaz',mugaz + +end subroutine iniphysiq + + +END MODULE iniphysiq_mod Index: trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/update_inputs_physiq_mod.F =================================================================== --- trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/update_inputs_physiq_mod.F (revision 3661) +++ trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/update_inputs_physiq_mod.F (revision 3661) @@ -0,0 +1,574 @@ +MODULE update_inputs_physiq_mod + +CONTAINS + +!SUBROUTINE update_inputs_physiq_time +!SUBROUTINE update_inputs_physiq_tracers +!SUBROUTINE update_inputs_physiq_constants +!SUBROUTINE update_inputs_physiq_geom +!SUBROUTINE update_inputs_physiq_surf +!SUBROUTINE update_inputs_physiq_soil +!SUBROUTINE update_inputs_physiq_turb +!SUBROUTINE update_inputs_physiq_rad +!SUBROUTINE update_inputs_physiq_slope + +!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! +!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! +SUBROUTINE update_inputs_physiq_time(& + JULYR,JULDAY,GMT,& + elaps,& + lct_input,lon_input,ls_input,& + MY) + USE variables_mod, only: JD_cur,JH_cur_split,phour_ini + use callkeys_mod, only : tlocked + INTEGER, INTENT(IN) :: JULDAY, JULYR + REAL, INTENT(IN) :: GMT,elaps,lon_input,ls_input,lct_input + REAL,INTENT(OUT) :: MY + REAL :: sec,nsec + !IF (JULYR .le. 8999) THEN + if (tlocked .eqv. .false.) THEN + JH_cur_split = (GMT + elaps/3600.) !! universal time (0