Changeset 4143 for LMDZ6/trunk/libf/dyn3d/iniacademic.F90
- Timestamp:
- May 9, 2022, 12:35:40 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/iniacademic.F90
r4124 r4143 5 5 6 6 USE filtreg_mod, ONLY: inifilr 7 USE infotrac, ONLY: nqtot, niso, niso_possibles, ok_iso_verif, tnat, alpha_ideal, & 8 iqiso, tracers, iso_indnum, iso_num 7 USE infotrac, ONLY: nqtot, niso, tnat, alpha_ideal, iqIsoPha, tracers 9 8 USE control_mod, ONLY: day_step,planet_type 10 9 use exner_hyb_m, only: exner_hyb … … 282 281 ! CRisi: init des isotopes 283 282 ! distill de Rayleigh très simplifiée 284 ! iName = tracers(iq)%iso_iName ! (next commit) 285 iName = iso_num(iq) 283 iName = tracers(iq)%iso_iName 286 284 if (niso <= 0 .OR. iName <= 0) CYCLE 287 285 iPhase = tracers(iq)%iso_iPhase … … 290 288 q(:,:,iq) = q(:,:,iqParent)*tnat(iName)*(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.) 291 289 ELSE 292 q(:,:,iq) = q(:,:,iq iso(iso_indnum(iq),iPhase))290 q(:,:,iq) = q(:,:,iqIsoPha(iName,iPhase)) 293 291 END IF 294 292 enddo … … 297 295 endif ! of if (planet_type=="earth") 298 296 299 if (ok_iso_verif)call check_isotopes_seq(q,1,ip1jmp1,'iniacademic_loc')297 call check_isotopes_seq(q,1,ip1jmp1,'iniacademic_loc') 300 298 301 299 ! add random perturbation to temperature
Note: See TracChangeset
for help on using the changeset viewer.