Changeset 3026 for trunk/LMDZ.MARS/libf/phymars/dyn1d/read_profile_mod.F90
- Timestamp:
- Aug 9, 2023, 11:56:36 AM (16 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/dyn1d/read_profile_mod.F90
r2991 r3026 61 61 !-------- 62 62 integer :: & 63 iq, & ! loop over nb_tracer64 ilayer, & ! loop over nb_layer65 ierr, & ! open file iostat66 indice_h2o_vap , & ! indice of h2o_vap tracer67 indice_h2o_ice , & ! indice of h2o_ice tracer68 indice_hdo_vap , & ! indice of hdo_vap tracer69 indice_hdo_ice ! indice of hdo_ice tracer63 iq, & ! loop over nb_tracer 64 ilayer, & ! loop over nb_layer 65 ierr, & ! open file iostat 66 indice_h2o_vap = 0, & ! indice of h2o_vap tracer 67 indice_h2o_ice = 0, & ! indice of h2o_ice tracer 68 indice_hdo_vap = 0, & ! indice of hdo_vap tracer 69 indice_hdo_ice = 0 ! indice of hdo_ice tracer 70 70 71 71 character(len=80), dimension(nb_tracer) :: & … … 148 148 ! the main isotopologue profile is already read before minors isotopologues 149 149 !---------------------------------------------------------------------------------------------------------------------! 150 if (hdo_vap .eqv..true. .and. indice_h2o_vap.ne.0) then150 if (hdo_vap .and. indice_h2o_vap /= 0) then 151 151 do ilayer = 1, nb_layer 152 152 q(ilayer, indice_hdo_vap) = q(ilayer, indice_h2o_vap)*2*155.76e-6*5 … … 154 154 end if 155 155 156 if (hdo_ice .eqv..true. .and. indice_h2o_ice.ne.0) then156 if (hdo_ice .and. indice_h2o_ice /= 0) then 157 157 qsurf(indice_hdo_ice) = qsurf(indice_h2o_ice) * 2*155.76e-6*5 158 158 do ilayer = 1, nb_layer
Note: See TracChangeset
for help on using the changeset viewer.