Changeset 2021 for trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars
- Timestamp:
- Oct 19, 2018, 6:16:41 PM (6 years ago)
- Location:
- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars
- Files:
-
- 3 edited
-
Registry.EM (modified) (2 diffs)
-
Registry.EM.newphys (modified) (2 diffs)
-
module_initialize_les.F (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM
r2018 r2021 147 147 ####MARS MARS MARS MARS MARS 148 148 ####MARS MARS MARS MARS MARS 149 149 ###VENUS TRACER 150 state real qco2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCO2" "CO2 number mixing ratio" "kg kg-1" 151 state real qco ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCO" "CO number mixing ratio" "kg kg-1" 152 state real qh2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2" "H2 number mixing ratio" "kg kg-1" 153 state real qo1d ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO1D" "O1d number mixing ratio" "kg kg-1" 154 state real qo ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO" "O number mixing ratio" "kg kg-1" 155 state real qo2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO2" "O2 number mixing ratio" "kg kg-1" 156 state real qo2dg ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO2dg" "O2dg number mixing ratio" "kg kg-1" 157 state real qo3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO3" "O3 number mixing ratio" "kg kg-1" 158 state real qh ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH" "H number mixing ratio" "kg kg-1" 159 state real qoh ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QOH" "OH number mixing ratio" "kg kg-1" 160 state real qho2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHO2" "HO2 number mixing ratio" "kg kg-1" 161 state real qh2o2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2O2" "H2O2 number mixing ratio" "kg kg-1" 162 state real qcl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCL" "Cl number mixing ratio" "kg kg-1" 163 state real qclo ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLO" "ClO number mixing ratio" "kg kg-1" 164 state real qcl2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCL2" "Cl2 number mixing ratio" "kg kg-1" 165 state real qhcl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHCL" "HCl number mixing ratio" "kg kg-1" 166 state real qhocl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHOCL" "HOCl number mixing ratio" "kg kg-1" 167 state real qclco ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLCO" "ClCO number mixing ratio" "kg kg-1" 168 state real qclco3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLCO3" "ClCO3 number mixing ratio" "kg kg-1" 169 state real qcocl2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCOCL2" "COCl2 number mixing ratio" "kg kg-1" 170 state real qss ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QS" "S number mixing ratio" "kg kg-1" 171 state real qso ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QSO" "SO number mixing ratio" "kg kg-1" 172 state real qso2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QSO2" "SO2 number mixing ratio" "kg kg-1" 173 state real qso3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QSO3" "SO3 number mixing ratio" "kg kg-1" 174 state real qs2o2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QS2O2" "S2O2 number mixing ratio" "kg kg-1" 175 state real qocs ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QOCS" "OCS number mixing ratio" "kg kg-1" 176 state real qhso3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHSO3" "HSO3 number mixing ratio" "kg kg-1" 177 state real qh2so4 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2SO4" "H2SO4 number mixing ratio" "kg kg-1" 178 state real qs2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QS2" "S2 number mixing ratio" "kg kg-1" 179 state real qclso2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLSO2" "ClSO2 number mixing ratio" "kg kg-1" 180 state real qoscl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QOSCL" "OSCl number mixing ratio" "kg kg-1" 181 state real qh2oliq ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2OLIQ" "H2O liq number mixing ratio" "kg kg-1" 182 state real qh2so4liq ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2SO4LIQ" "H2SO4 liq number mixing ratio" "kg kg-1" 150 183 151 184 … … 1517 1550 package radioac mars==20 - scalar:qtrac1 1518 1551 package radioac2 mars==21 - scalar:upward,downward 1552 package VenusPC mars==34 - scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq 1519 1553 ##### MARS OPTIONS 1520 1554 ##### MARS OPTIONS -
trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys
r2018 r2021 144 144 ####MARS MARS MARS MARS MARS 145 145 ####MARS MARS MARS MARS MARS 146 ###VENUS TRACER 147 state real qco2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCO2" "CO2 number mixing ratio" "kg kg-1" 148 state real qco ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCO" "CO number mixing ratio" "kg kg-1" 149 state real qh2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2" "H2 number mixing ratio" "kg kg-1" 150 state real qo1d ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO1D" "O1d number mixing ratio" "kg kg-1" 151 state real qo ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO" "O number mixing ratio" "kg kg-1" 152 state real qo2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO2" "O2 number mixing ratio" "kg kg-1" 153 state real qo2dg ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO2dg" "O2dg number mixing ratio" "kg kg-1" 154 state real qo3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QO3" "O3 number mixing ratio" "kg kg-1" 155 state real qh ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH" "H number mixing ratio" "kg kg-1" 156 state real qoh ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QOH" "OH number mixing ratio" "kg kg-1" 157 state real qho2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHO2" "HO2 number mixing ratio" "kg kg-1" 158 state real qh2o2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2O2" "H2O2 number mixing ratio" "kg kg-1" 159 state real qcl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCL" "Cl number mixing ratio" "kg kg-1" 160 state real qclo ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLO" "ClO number mixing ratio" "kg kg-1" 161 state real qcl2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCL2" "Cl2 number mixing ratio" "kg kg-1" 162 state real qhcl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHCL" "HCl number mixing ratio" "kg kg-1" 163 state real qhocl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHOCL" "HOCl number mixing ratio" "kg kg-1" 164 state real qclco ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLCO" "ClCO number mixing ratio" "kg kg-1" 165 state real qclco3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLCO3" "ClCO3 number mixing ratio" "kg kg-1" 166 state real qcocl2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCOCL2" "COCl2 number mixing ratio" "kg kg-1" 167 state real qss ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QS" "S number mixing ratio" "kg kg-1" 168 state real qso ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QSO" "SO number mixing ratio" "kg kg-1" 169 state real qso2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QSO2" "SO2 number mixing ratio" "kg kg-1" 170 state real qso3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QSO3" "SO3 number mixing ratio" "kg kg-1" 171 state real qs2o2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QS2O2" "S2O2 number mixing ratio" "kg kg-1" 172 state real qocs ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QOCS" "OCS number mixing ratio" "kg kg-1" 173 state real qhso3 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QHSO3" "HSO3 number mixing ratio" "kg kg-1" 174 state real qh2so4 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2SO4" "H2SO4 number mixing ratio" "kg kg-1" 175 state real qs2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QS2" "S2 number mixing ratio" "kg kg-1" 176 state real qclso2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCLSO2" "ClSO2 number mixing ratio" "kg kg-1" 177 state real qoscl ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QOSCL" "OSCl number mixing ratio" "kg kg-1" 178 state real qh2oliq ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2OLIQ" "H2O liq number mixing ratio" "kg kg-1" 179 state real qh2so4liq ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2SO4LIQ" "H2SO4 liq number mixing ratio" "kg kg-1" 180 146 181 147 182 … … 1516 1551 package newwatersca mars==12 - scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn 1517 1552 package radioac mars==20 - scalar:qtrac1 1553 package VenusPC mars==34 - scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq 1518 1554 ##### MARS OPTIONS 1519 1555 ##### MARS OPTIONS -
trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F
r2020 r2021 87 87 INTEGER, PARAMETER :: nl_max = 1000 88 88 REAL, DIMENSION(nl_max) :: zk, p_in, theta, rho, u, v, qv, pd_in 89 REAL*8, DIMENSION(nl_max) :: pd_in8 89 90 INTEGER :: nl_in 90 91 91 92 92 INTEGER :: icm,jcm, ii, im1, jj, jm1, loop, error, fid, nxc, nyc 93 93 REAL :: u_mean,v_mean, f0, p_surf, p_level, qvf, z_at_v, z_at_u 94 REAL*8 :: p_level8 94 95 REAL :: xrad, yrad, zrad, rad, delt, cof1, cof2 95 96 ! REAL, EXTERNAL :: interp_0 … … 121 122 REAL, DIMENSION(nl_max) :: venus_hrdyn 122 123 REAL, DIMENSION(nl_max) :: altitude 124 REAL*8, DIMENSION(nl_max) :: trac 123 125 !!MARS 124 126 … … 969 971 ENDIF 970 972 973 IF ((planet.eq."venus") .AND. ( config_flags%mars == 34 )) Then 974 pd_in8(:)=pd_in(:) 975 do i = 1,34 976 call read_tracer(trac,num_scalar,i,nl_in) 977 DO k=1,kte-1 978 p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top 979 p_level8=p_level 980 scalar(its:ite,k,jts:jte,i+1) = interp_0_log2( trac, pd_in8, p_level8, nl_in ) 981 ENDDO 982 ENDDO 983 !close(unit=22) 984 ENDIF 971 985 972 986 open(unit=21,file="altitude.txt",action="write") … … 1358 1372 end subroutine read_lsf 1359 1373 1374 subroutine read_tracer(trace,nq,qn,n) 1375 implicit none 1376 integer n,qn,nq ! qn : number of the tracer 1377 real*8 tra(nq-1,n+1) 1378 real*8 trace(n+1) !output 1379 logical end_of_file 1380 1381 integer k,j 1382 1383 ! first element is the surface 1384 open(unit=14,file='input_tracer',form='formatted',status='old') 1385 rewind(14) 1386 end_of_file = .false. 1387 DO k=1,n 1388 read(14,*) tra(:,k) 1389 write(*,*) k,tra(qn,k) 1390 ENDDO 1391 1392 close(14) 1393 trace(:)=tra(qn,:) 1394 end subroutine read_tracer 1395 1396 1360 1397 END MODULE module_initialize_ideal
Note: See TracChangeset
for help on using the changeset viewer.
