Index: trunk/LMDZ.MARS/libf/phymars/callkeys.h =================================================================== --- trunk/LMDZ.MARS/libf/phymars/callkeys.h (revision 228) +++ trunk/LMDZ.MARS/libf/phymars/callkeys.h (revision 234) @@ -12,5 +12,5 @@ & ,callg2d,linear,rayleigh,tracer,active,doubleq,submicron & & ,lifting,callddevil,scavenging,sedimentation,activice,water & - & ,caps,photochem,calltherm,outptherm + & ,caps,photochem,calltherm,outptherm,callslope COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche & @@ -24,5 +24,5 @@ & ,callnirco2,callnlte,callthermos,callconduct, & & calleuv,callmolvis,callmoldiff,thermochem,thermoswater & - & ,calltherm,outptherm + & ,calltherm,outptherm,callslope Index: trunk/LMDZ.MARS/libf/phymars/inifis.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/inifis.F (revision 228) +++ trunk/LMDZ.MARS/libf/phymars/inifis.F (revision 234) @@ -4,4 +4,7 @@ $ ,plat,plon,parea $ ,prad,pg,pr,pcpp +#ifdef MESOSCALE +#include "meso_inc/meso_inc_inifisinvar.F" +#endif $ ) ! @@ -21,4 +24,5 @@ ! stored in the q(:,:,:,:) array ! E.M. (june 2009) use getin routine to load parameters +! adapted to the mesoscale use - Aymeric Spiga - 01/2007-07/2011 ! ! @@ -56,5 +60,9 @@ #include "yomaer.h" #include "datafile.h" - +#include "slope.h" +#ifdef MESOSCALE +#include "comsoil.h" !!MESOSCALE -- needed to fill volcapa +#include "meso_inc/meso_inc_inifisvar.F" +#endif REAL prad,pg,pr,pcpp,pdaysec @@ -85,4 +93,7 @@ daysec=pdaysec dtphys=ptimestep +#ifdef MESOSCALE +#include "meso_inc/meso_inc_inifisini.F" +#endif ! -------------------------------------------------------- @@ -150,5 +161,9 @@ write(*,*) "Save statistics in file stats.nc ?" +#ifdef MESOSCALE + callstats=.false. ! default value +#else callstats=.true. ! default value +#endif call getin("callstats",callstats) write(*,*) " callstats = ",callstats @@ -190,4 +205,14 @@ call getin("callrad",callrad) write(*,*) " callrad = ",callrad + + write(*,*) "call slope insolation scheme ?", + & "(matters only if callrad=T)" +#ifdef MESOSCALE + callslope=.true. ! default value +#else + callslope=.false. ! default value (not supported yet) +#endif + call getin("callslope",callslope) + write(*,*) " callslope = ",callslope write(*,*) "call NLTE radiative schemes ?", Index: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisini.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisini.F (revision 228) +++ trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisini.F (revision 234) @@ -8,5 +8,4 @@ c ! in 'comcstfi.h' - daysec=wdaysec omeg=womeg mugaz=wmugaz @@ -89,4 +88,4 @@ c It must be set now, because it is used afterwards c***************************************************** - dtphys=wdt*float(wappel_phys) + dtphys=wdt*float(ptimestep) print*,'Physical timestep (s) ',dtphys Index: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisinvar.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisinvar.F (revision 228) +++ trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisinvar.F (revision 234) @@ -1,2 +1,3 @@ + $ ,nq,wdt $ ,womeg,wmugaz $ ,wyear_day,wperiheli,waphelie,wperi_day,wobliquit Index: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisvar.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisvar.F (revision 228) +++ trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisvar.F (revision 234) @@ -1,5 +1,6 @@ - INTEGER nq,wday_ini + INTEGER nq + REAL wdt - REAL womeg,wmugaz,wdaysec + REAL womeg,wmugaz REAL wyear_day,wperiheli,waphelie,wperi_day,wobliquit REAL wz0,wemin_turb,wlmixmin @@ -12,4 +13,2 @@ REAL wtheta(ngrid),wpsi(ngrid) REAL wvolcapa - REAL wdt - INTEGER wappel_phys Index: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_slope.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_slope.F (revision 228) +++ (revision ) @@ -1,54 +1,0 @@ -cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -ccccc -ccccc PARAM SLOPE : Insolation (direct + scattered) -ccccc - DO ig=1,ngrid - sl_the = theta_sl(ig) - IF (sl_the .ne. 0.) THEN - ztim1=fluxsurf_sw(ig,1)+fluxsurf_sw(ig,2) - DO l=1,2 - sl_lct = ptime*24. + 180.*long(ig)/pi/15. - sl_ra = pi*(1.0-sl_lct/12.) - sl_lat = 180.*lati(ig)/pi - sl_tau = tau(ig,1) - sl_alb = albedo(ig,l) - sl_psi = psi_sl(ig) - sl_fl0 = fluxsurf_sw(ig,l) - sl_di0 = 0. - if (mu0(ig) .gt. 0.) then - sl_di0 = mu0(ig)*(exp(-sl_tau/mu0(ig))) - sl_di0 = sl_di0*1370./dist_sol/dist_sol - sl_di0 = sl_di0/ztim1 - sl_di0 = fluxsurf_sw(ig,l)*sl_di0 - endif - ! sait-on jamais (a cause des arrondis) - if (sl_fl0 .lt. sl_di0) sl_di0=sl_fl0 - !!!!!!!!!!!!!!!!!!!!!!!!!! - CALL param_slope( mu0(ig), declin, sl_ra, sl_lat, - & sl_tau, sl_alb, - & sl_the, sl_psi, sl_di0, sl_fl0, sl_flu) - !!!!!!!!!!!!!!!!!!!!!!!!!! - fluxsurf_sw(ig,l) = sl_flu - !! sl_ls = 180.*zls/pi - !! sl_lct = ptime*24. + 180.*long(ig)/pi/15. - !! sl_lat = 180.*lati(ig)/pi - !! sl_tau = tau(ig,1) - !! sl_alb = albedo(ig,l) - !! sl_the = theta_sl(ig) - !! sl_psi = psi_sl(ig) - !! sl_fl0 = fluxsurf_sw(ig,l) - !! CALL param_slope_full(sl_ls, sl_lct, sl_lat, - !! & sl_tau, sl_alb, - !! & sl_the, sl_psi, sl_fl0, sl_flu) - ENDDO - !!! compute correction on IR flux as well - sky= (1.+cos(pi*theta_sl(ig)/180.))/2. - fluxsurf_lw(ig)= fluxsurf_lw(ig)*sky - ENDIF - ENDDO -ccccc -ccccc PARAM SLOPE -ccccc -cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc Index: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F (revision 228) +++ trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F (revision 234) @@ -10,8 +10,4 @@ REAL output_tab2d(ngridmx,n2d) REAL output_tab3d(ngridmx,nlayer,n3d) - REAL sl_ls, sl_lct, sl_lat, sl_tau, sl_alb, sl_the, sl_psi - REAL sl_fl0, sl_flu - REAL sl_ra, sl_di0 - REAL sky REAL hfx(ngridmx) !! pour LES avec isfflx!=0 REAL ust(ngridmx) !! pour LES avec isfflx!=0 Index: trunk/LMDZ.MARS/libf/phymars/meso_inifis.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_inifis.F (revision 228) +++ (revision ) @@ -1,694 +1,0 @@ - SUBROUTINE meso_inifis( - $ ngrid,nlayer -#ifdef MESOSCALE - $ ,nq,wdt,wday_ini,wdaysec,wappel_phys -#else - $ ,day_ini,pdaysec,ptimestep -#endif - $ ,plat,plon,parea - $ ,prad,pg,pr,pcpp -#ifdef MESOSCALE -#include "meso_inc/meso_inc_inifisinvar.F" -#endif - $ ) -! -!======================================================================= -! -! purpose: -! ------- -! -! Initialisation for the physical parametrisations of the LMD -! martian atmospheric general circulation modele. -! -! author: Frederic Hourdin 15 / 10 /93 -! ------- -! modified: Sebastien Lebonnois 11/06/2003 (new callphys.def) -! Ehouarn Millour (oct. 2008) tracers are now identified -! by their names and may not be contiguously -! stored in the q(:,:,:,:) array -! E.M. (june 2009) use getin routine to load parameters -! adapted to the WRF use - Aymeric Spiga - Jan 2009 - Jan 2007 -! -! -! arguments: -! ---------- -! -! input: -! ------ -! -! ngrid Size of the horizontal grid. -! All internal loops are performed on that grid. -! nlayer Number of vertical layers. -! pdayref Day of reference for the simulation -! pday Number of days counted from the North. Spring -! equinoxe. -! -!======================================================================= -! -!----------------------------------------------------------------------- -! declarations: -! ------------- -! to use 'getin' - USE ioipsl_getincom - IMPLICIT NONE -#include "dimensions.h" -#include "dimphys.h" -#include "planete.h" -#include "comcstfi.h" -#include "comsaison.h" -#include "comdiurn.h" -#include "comgeomfi.h" -#include "callkeys.h" -#include "surfdat.h" -#include "dimradmars.h" -#include "yomaer.h" -#include "datafile.h" -#ifdef MESOSCALE -#include "slope.h" -#include "comsoil.h" !!MESOSCALE -- needed to fill volcapa -#include "meso_inc/meso_inc_inifisvar.F" -#endif - REAL prad,pg,pr,pcpp,pdaysec -#ifndef MESOSCALE - REAL ptimestep - INTEGER day_ini -#endif - INTEGER ngrid,nlayer - REAL plat(ngrid),plon(ngrid),parea(ngridmx) - INTEGER ig,ierr - -! EXTERNAL iniorbit,orbite - EXTERNAL SSUM - REAL SSUM - - CHARACTER ch1*12 - CHARACTER ch80*80 - -! logical chem, h2o - -! chem = .false. -! h2o = .false. - - rad=prad - cpp=pcpp - g=pg - r=pr - rcp=r/cpp -#ifndef MESOSCALE - daysec=pdaysec - dtphys=ptimestep -#else -#include "meso_inc/meso_inc_inifisini.F" -#endif - -! -------------------------------------------------------- -! The usual Tests -! -------------- - IF (nlayer.NE.nlayermx) THEN - PRINT*,'STOP in inifis' - PRINT*,'Probleme de dimensions :' - PRINT*,'nlayer = ',nlayer - PRINT*,'nlayermx = ',nlayermx - STOP - ENDIF - - IF (ngrid.NE.ngridmx) THEN - PRINT*,'STOP in inifis' - PRINT*,'Probleme de dimensions :' - PRINT*,'ngrid = ',ngrid - PRINT*,'ngridmx = ',ngridmx - STOP - ENDIF - -! -------------------------------------------------------------- -! Reading the "callphys.def" file controlling some key options -! -------------------------------------------------------------- - - ! check that 'callphys.def' file is around - OPEN(99,file='callphys.def',status='old',form='formatted' - & ,iostat=ierr) - CLOSE(99) - - IF(ierr.EQ.0) THEN - PRINT* - PRINT* - PRINT*,'--------------------------------------------' - PRINT*,' inifis: Parameters for the physics (callphys.def)' - PRINT*,'--------------------------------------------' - - write(*,*) "Directory where external input files are:" - datafile="/u/forget/WWW/datagcm/datafile" - call getin("datadir",datafile) ! default path - write(*,*) " datafile = ",trim(datafile) - - write(*,*) "Run with or without tracer transport ?" - tracer=.false. ! default value - call getin("tracer",tracer) - write(*,*) " tracer = ",tracer - - write(*,*) "Diurnal cycle ?" - write(*,*) "(if diurnal=False, diurnal averaged solar heating)" - diurnal=.true. ! default value - call getin("diurnal",diurnal) - write(*,*) " diurnal = ",diurnal - - write(*,*) "Seasonal cycle ?" - write(*,*) "(if season=False, Ls stays constant, to value ", - & "set in 'start'" - season=.true. ! default value - call getin("season",season) - write(*,*) " season = ",season - - write(*,*) "Write some extra output to the screen ?" - lwrite=.false. ! default value - call getin("lwrite",lwrite) - write(*,*) " lwrite = ",lwrite - - write(*,*) "Save statistics in file stats.nc ?" -#ifdef MESOSCALE - callstats=.false. ! default value -#else - callstats=.true. ! default value -#endif - call getin("callstats",callstats) - write(*,*) " callstats = ",callstats - - write(*,*) "Save EOF profiles in file 'profiles' for ", - & "Climate Database?" - calleofdump=.false. ! default value - call getin("calleofdump",calleofdump) - write(*,*) " calleofdump = ",calleofdump - - write(*,*) "Dust scenario: 1=constant dust (read from startfi", - & " or set as tauvis); 2=Viking scenario; =3 MGS scenario,", - & "=4 Mars Year 24 from TES assimilation, ", - & "=24,25 or 26 :Mars Year 24,25 or 26 from TES assimilation" - iaervar=3 ! default value - call getin("iaervar",iaervar) - write(*,*) " iaervar = ",iaervar - - write(*,*) "Reference (visible) dust opacity at 700 Pa ", - & "(matters only if iaervar=1)" - ! NB: default value of tauvis is set/read in startfi.nc file - call getin("tauvis",tauvis) - write(*,*) " tauvis = ",tauvis - - write(*,*) "Dust vertical distribution:" - write(*,*) "(=1 top set by topdustref parameter;", - & " =2 Viking scenario; =3 MGS scenario)" - iddist=3 ! default value - call getin("iddist",iddist) - write(*,*) " iddist = ",iddist - - write(*,*) "Dust top altitude (km). (Matters only if iddist=1)" - topdustref= 90.0 ! default value - call getin("topdustref",topdustref) - write(*,*) " topdustref = ",topdustref - - write(*,*) "call radiative transfer ?" - callrad=.true. ! default value - call getin("callrad",callrad) - write(*,*) " callrad = ",callrad - - write(*,*) "call NLTE radiative schemes ?", - & "(matters only if callrad=T)" - callnlte=.false. ! default value - call getin("callnlte",callnlte) - write(*,*) " callnlte = ",callnlte - - write(*,*) "call CO2 NIR absorption ?", - & "(matters only if callrad=T)" - callnirco2=.false. ! default value - call getin("callnirco2",callnirco2) - write(*,*) " callnirco2 = ",callnirco2 - - write(*,*) "call turbulent vertical diffusion ?" - calldifv=.true. ! default value - call getin("calldifv",calldifv) - write(*,*) " calldifv = ",calldifv - - write(*,*) "call thermals ?" - calltherm=.false. ! default value - call getin("calltherm",calltherm) - write(*,*) " calltherm = ",calltherm - - write(*,*) "output thermal diagnostics ?" - outptherm=.false. ! default value - call getin("outptherm",outptherm) - write(*,*) " outptherm = ",outptherm - - write(*,*) "call convective adjustment ?" - calladj=.true. ! default value - call getin("calladj",calladj) - write(*,*) " calladj = ",calladj - - if (calltherm .and. (.not. calladj)) then - print*,'Convadj has to be activated when using thermals' - stop - endif - - write(*,*) "call CO2 condensation ?" - callcond=.true. ! default value - call getin("callcond",callcond) - write(*,*) " callcond = ",callcond - - write(*,*)"call thermal conduction in the soil ?" - callsoil=.true. ! default value - call getin("callsoil",callsoil) - write(*,*) " callsoil = ",callsoil - - - write(*,*)"call Lott's gravity wave/subgrid topography ", - & "scheme ?" - calllott=.true. ! default value - call getin("calllott",calllott) - write(*,*)" calllott = ",calllott - - - write(*,*)"rad.transfer is computed every iradia", - & " physical timestep" - iradia=1 ! default value - call getin("iradia",iradia) - write(*,*)" iradia = ",iradia - - - write(*,*)"Output of the exchange coefficient mattrix ?", - & "(for diagnostics only)" - callg2d=.false. ! default value - call getin("callg2d",callg2d) - write(*,*)" callg2d = ",callg2d - - write(*,*)"Rayleigh scattering : (should be .false. for now)" - rayleigh=.false. - call getin("rayleigh",rayleigh) - write(*,*)" rayleigh = ",rayleigh - - -! TRACERS: - - write(*,*)"Transported dust ? (if >0, use 'dustbin' dust bins)" - dustbin=0 ! default value - call getin("dustbin",dustbin) - write(*,*)" dustbin = ",dustbin - - write(*,*)"Radiatively active dust ? (matters if dustbin>0)" - active=.false. ! default value - call getin("active",active) - write(*,*)" active = ",active - -! Test of incompatibility: -! if active is used, then dustbin should be > 0 - - if (active.and.(dustbin.lt.1)) then - print*,'if active is used, then dustbin should > 0' - stop - endif - - write(*,*)"use mass and number mixing ratios to predict", - & " dust size ?" - doubleq=.false. ! default value - call getin("doubleq",doubleq) - write(*,*)" doubleq = ",doubleq - - submicron=.false. ! default value - call getin("submicron",submicron) - write(*,*)" submicron = ",submicron - -! Test of incompatibility: -! if doubleq is used, then dustbin should be 2 - - if (doubleq.and.(dustbin.ne.2)) then - print*,'if doubleq is used, then dustbin should be 2' - stop - endif - if (doubleq.and.submicron.and.(nqmx.LT.3)) then - print*,'If doubleq is used with a submicron tracer,' - print*,' then the number of tracers has to be' - print*,' larger than 3.' - stop - endif - - write(*,*)"dust lifted by GCM surface winds ?" - lifting=.false. ! default value - call getin("lifting",lifting) - write(*,*)" lifting = ",lifting - -! Test of incompatibility: -! if lifting is used, then dustbin should be > 0 - - if (lifting.and.(dustbin.lt.1)) then - print*,'if lifting is used, then dustbin should > 0' - stop - endif - - write(*,*)" dust lifted by dust devils ?" - callddevil=.false. !default value - call getin("callddevil",callddevil) - write(*,*)" callddevil = ",callddevil - - -! Test of incompatibility: -! if dustdevil is used, then dustbin should be > 0 - - if (callddevil.and.(dustbin.lt.1)) then - print*,'if dustdevil is used, then dustbin should > 0' - stop - endif - - write(*,*)"Dust scavenging by CO2 snowfall ?" - scavenging=.false. ! default value - call getin("scavenging",scavenging) - write(*,*)" scavenging = ",scavenging - - -! Test of incompatibility: -! if scavenging is used, then dustbin should be > 0 - - if (scavenging.and.(dustbin.lt.1)) then - print*,'if scavenging is used, then dustbin should > 0' - stop - endif - - write(*,*) "Gravitationnal sedimentation ?" - sedimentation=.true. ! default value - call getin("sedimentation",sedimentation) - write(*,*) " sedimentation = ",sedimentation - - write(*,*) "Radiatively active transported atmospheric ", - & "water ice ?" - activice=.false. ! default value - call getin("activice",activice) - write(*,*) " activice = ",activice - - write(*,*) "Compute water cycle ?" - water=.false. ! default value - call getin("water",water) - write(*,*) " water = ",water - -! Test of incompatibility: - - if (activice.and..not.water) then - print*,'if activice is used, water should be used too' - stop - endif - - if (water.and..not.tracer) then - print*,'if water is used, tracer should be used too' - stop - endif - -! Test of incompatibility: - - write(*,*) "Permanent water caps at poles ?", - & " .true. is RECOMMENDED" - write(*,*) "(with .true., North cap is a source of water ", - & "and South pole is a cold trap)" - caps=.true. ! default value - call getin("caps",caps) - write(*,*) " caps = ",caps - - write(*,*) "photochemistry: include chemical species" - photochem=.false. ! default value - call getin("photochem",photochem) - write(*,*) " photochem = ",photochem - - -! THERMOSPHERE - - write(*,*) "call thermosphere ?" - callthermos=.false. ! default value - call getin("callthermos",callthermos) - write(*,*) " callthermos = ",callthermos - - - write(*,*) " water included without cycle ", - & "(only if water=.false.)" - thermoswater=.false. ! default value - call getin("thermoswater",thermoswater) - write(*,*) " thermoswater = ",thermoswater - - write(*,*) "call thermal conduction ?", - & " (only if callthermos=.true.)" - callconduct=.false. ! default value - call getin("callconduct",callconduct) - write(*,*) " callconduct = ",callconduct - - write(*,*) "call EUV heating ?", - & " (only if callthermos=.true.)" - calleuv=.false. ! default value - call getin("calleuv",calleuv) - write(*,*) " calleuv = ",calleuv - - write(*,*) "call molecular viscosity ?", - & " (only if callthermos=.true.)" - callmolvis=.false. ! default value - call getin("callmolvis",callmolvis) - write(*,*) " callmolvis = ",callmolvis - - write(*,*) "call molecular diffusion ?", - & " (only if callthermos=.true.)" - callmoldiff=.false. ! default value - call getin("callmoldiff",callmoldiff) - write(*,*) " callmoldiff = ",callmoldiff - - - write(*,*) "call thermospheric photochemistry ?", - & " (only if callthermos=.true.)" - thermochem=.false. ! default value - call getin("thermochem",thermochem) - write(*,*) " thermochem = ",thermochem - - write(*,*) "date for solar flux calculation:", - & " (1985 < date < 2002)" - write(*,*) "(Solar min=1996.4 ave=1993.4 max=1990.6)" - solarcondate=1993.4 ! default value - call getin("solarcondate",solarcondate) - write(*,*) " solarcondate = ",solarcondate - - - if (.not.callthermos) then - if (thermoswater) then - print*,'if thermoswater is set, callthermos must be true' - stop - endif - if (callconduct) then - print*,'if callconduct is set, callthermos must be true' - stop - endif - if (calleuv) then - print*,'if calleuv is set, callthermos must be true' - stop - endif - if (callmolvis) then - print*,'if callmolvis is set, callthermos must be true' - stop - endif - if (callmoldiff) then - print*,'if callmoldiff is set, callthermos must be true' - stop - endif - if (thermochem) then - print*,'if thermochem is set, callthermos must be true' - stop - endif - endif - -! Test of incompatibility: -! if photochem is used, then water should be used too - - if (photochem.and..not.water) then - print*,'if photochem is used, water should be used too' - stop - endif - -! if callthermos is used, then thermoswater should be used too -! (if water not used already) - - if (callthermos .and. .not.water) then - if (callthermos .and. .not.thermoswater) then - print*,'if callthermos is used, water or thermoswater - & should be used too' - stop - endif - endif - - PRINT*,'--------------------------------------------' - PRINT* - PRINT* - ELSE - write(*,*) - write(*,*) 'Cannot read file callphys.def. Is it here ?' - stop - ENDIF - -8000 FORMAT(t5,a12,l8) -8001 FORMAT(t5,a12,i8) - - PRINT* - PRINT*,'inifis: daysec',daysec - PRINT* - PRINT*,'inifis: The radiative transfer is computed:' - PRINT*,' each ',iradia,' physical time-step' - PRINT*,' or each ',iradia*dtphys,' seconds' - PRINT* -! -------------------------------------------------------------- -! Managing the Longwave radiative transfer -! -------------------------------------------------------------- - -! In most cases, the run just use the following values : -! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - callemis=.true. -! ilwd=10*int(daysec/dtphys) ! bug before 22/10/01 - ilwd=1 - ilwn=1 !2 - ilwb=1 !2 - linear=.true. - ncouche=3 - alphan=0.4 - semi=0 - -! BUT people working hard on the LW may want to read them in 'radia.def' -! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - - OPEN(99,file='radia.def',status='old',form='formatted' - . ,iostat=ierr) - IF(ierr.EQ.0) THEN - write(*,*) 'inifis: Reading radia.def !!!' - READ(99,fmt='(a)') ch1 - READ(99,*) callemis - WRITE(*,8000) ch1,callemis - - READ(99,fmt='(a)') ch1 - READ(99,*) iradia - WRITE(*,8001) ch1,iradia - - READ(99,fmt='(a)') ch1 - READ(99,*) ilwd - WRITE(*,8001) ch1,ilwd - - READ(99,fmt='(a)') ch1 - READ(99,*) ilwn - WRITE(*,8001) ch1,ilwn - - READ(99,fmt='(a)') ch1 - READ(99,*) linear - WRITE(*,8000) ch1,linear - - READ(99,fmt='(a)') ch1 - READ(99,*) ncouche - WRITE(*,8001) ch1,ncouche - - READ(99,fmt='(a)') ch1 - READ(99,*) alphan - WRITE(*,*) ch1,alphan - - READ(99,fmt='(a)') ch1 - READ(99,*) ilwb - WRITE(*,8001) ch1,ilwb - - - READ(99,fmt='(a)') ch1 - READ(99,'(l1)') callg2d - WRITE(*,8000) ch1,callg2d - - READ(99,fmt='(a)') ch1 - READ(99,*) semi - WRITE(*,*) ch1,semi - end if - CLOSE(99) - -!----------------------------------------------------------------------- -! Some more initialization: -! ------------------------ - - ! in 'comgeomfi.h' - CALL SCOPY(ngrid,plon,1,long,1) - CALL SCOPY(ngrid,plat,1,lati,1) - CALL SCOPY(ngrid,parea,1,area,1) - totarea=SSUM(ngridmx,area,1) - - ! in 'comdiurn.h' - DO ig=1,ngrid - sinlat(ig)=sin(plat(ig)) - coslat(ig)=cos(plat(ig)) - sinlon(ig)=sin(plon(ig)) - coslon(ig)=cos(plon(ig)) - ENDDO - - pi=2.*asin(1.) ! NB: pi is a common in comcstfi.h - -! managing the tracers, and tests: -! ------------------------------- -! Ehouarn: removed; as these tests are now done in initracer.F -! if(tracer) then -! -!! when photochem is used, nqchem_min is the rank -!! of the first chemical species -! -!! Ehouarn: nqchem_min is now meaningless and no longer used -!! nqchem_min = 1 -! if (photochem .or. callthermos) then -! chem = .true. -! end if -! -! if (water .or. thermoswater) h2o = .true. -! -!! TESTS -! -! print*,'inifis: TRACERS:' -! write(*,*) " chem=",chem," h2o=",h2o -!! write(*,*) " doubleq=",doubleq -!! write(*,*) " dustbin=",dustbin -! -! if ((doubleq).and.(h2o).and. -! $ (chem)) then -! print*,' 2 dust tracers (doubleq)' -! print*,' 1 water vapour tracer' -! print*,' 1 water ice tracer' -! print*,nqmx-4,' chemistry tracers' -! endif -! -! if ((doubleq).and.(h2o).and. -! $ .not.(chem)) then -! print*,' 2 dust tracers (doubleq)' -! print*,' 1 water vapour tracer' -! print*,' 1 water ice tracer' -! if (nqmx.LT.4) then -! print*,'nqmx should be at least equal to' -! print*,'4 with these options.' -! stop -! endif -! endif -! -! if (.not.(doubleq).and.(h2o).and. -! $ (chem)) then -! if (dustbin.gt.0) then -! print*,dustbin,' dust bins' -! endif -! print*,nqmx-2-dustbin,' chemistry tracers' -! print*,' 1 water vapour tracer' -! print*,' 1 water ice tracer' -! endif -! -! if (.not.(doubleq).and.(h2o).and. -! $ .not.(chem)) then -! if (dustbin.gt.0) then -! print*,dustbin,' dust bins' -! endif -! print*,' 1 water vapour tracer' -! print*,' 1 water ice tracer' -! if (nqmx.gt.(dustbin+2)) then -! print*,'nqmx should be ',(dustbin+2), -! $ ' with these options...' -! print*,'(or check callphys.def)' -! endif -! if (nqmx.lt.(dustbin+2)) then -! write(*,*) "inifis: nqmx.lt.(dustbin+2)" -! stop -! endif -! endif -! -! endif ! of if (tracer) -! -! RETURN - END Index: trunk/LMDZ.MARS/libf/phymars/meso_physiq.F =================================================================== --- trunk/LMDZ.MARS/libf/phymars/meso_physiq.F (revision 228) +++ (revision ) @@ -1,1680 +1,0 @@ - SUBROUTINE meso_physiq( - $ ngrid,nlayer,nq - $ ,firstcall,lastcall - $ ,pday,ptime,ptimestep - $ ,pplev,pplay,pphi - $ ,pu,pv,pt,pq - $ ,pw - $ ,pdu,pdv,pdt,pdq,pdpsrf,tracerdyn -#ifdef MESOSCALE -#include "meso_inc/meso_inc_invar.F" -#endif - $ ) - - IMPLICIT NONE -c======================================================================= -c -c subject: -c -------- -c -c Organisation of the physical parametrisations of the LMD -c martian atmospheric general circulation model. -c -c The GCM can be run without or with tracer transport -c depending on the value of Logical "tracer" in file "callphys.def" -c Tracers may be water vapor, ice OR chemical species OR dust particles -c -c SEE comments in initracer.F about numbering of tracer species... -c -c It includes: -c -c 1. Initialization: -c 1.1 First call initializations -c 1.2 Initialization for every call to physiq -c 1.2.5 Compute mean mass and cp, R and thermal conduction coeff. -c 2. Compute radiative transfer tendencies -c (longwave and shortwave) for CO2 and aerosols. -c 3. Gravity wave and subgrid scale topography drag : -c 4. Vertical diffusion (turbulent mixing): -c 5. Convective adjustment -c 6. Condensation and sublimation of carbon dioxide. -c 7. TRACERS : -c 7a. water and water ice -c 7b. call for photochemistry when tracers are chemical species -c 7c. other scheme for tracer (dust) transport (lifting, sedimentation) -c 7d. updates (CO2 pressure variations, surface budget) -c 8. Contribution to tendencies due to thermosphere -c 9. Surface and sub-surface temperature calculations -c 10. Write outputs : -c - "startfi", "histfi" (if it's time) -c - Saving statistics (if "callstats = .true.") -c - Dumping eof (if "calleofdump = .true.") -c - Output any needed variables in "diagfi" -c 11. Diagnostic: mass conservation of tracers -c -c author: -c ------- -c Frederic Hourdin 15/10/93 -c Francois Forget 1994 -c Christophe Hourdin 02/1997 -c Subroutine completly rewritten by F.Forget (01/2000) -c Introduction of the photochemical module: S. Lebonnois (11/2002) -c Introduction of the thermosphere module: M. Angelats i Coll (2002) -c Water ice clouds: Franck Montmessin (update 06/2003) -c Radiatively active tracers: J.-B. Madeleine (10/2008-06/2009) -c Nb: See callradite.F for more information. -c Mesoscale version: Aymeric Spiga (2007 - 2011) -- check MESOSCALE flags -c -c arguments: -c ---------- -c -c input: -c ------ -c ecri period (in dynamical timestep) to write output -c ngrid Size of the horizontal grid. -c All internal loops are performed on that grid. -c nlayer Number of vertical layers. -c nq Number of advected fields -c firstcall True at the first call -c lastcall True at the last call -c pday Number of days counted from the North. Spring -c equinoxe. -c ptime Universal time (0