Changeset 5255 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
- Timestamp:
- Oct 22, 2024, 3:35:01 PM (4 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r5253 r5255 329 329 d_deltat_ajs_cv, d_deltaq_ajs_cv, & ! due to dry adjustment of (w) before convection 330 330 ! tendencies of wake fractional area and wake number per unit area: 331 d_s_wk, d_dens_a_wk, d_dens_wk, & ! due to wakes331 d_s_wk, d_s_a_wk, d_dens_a_wk, d_dens_wk, & ! due to wakes 332 332 !!! d_s_vdf, d_dens_a_vdf, d_dens_vdf, & ! due to vertical diffusion 333 333 !!! d_s_the, d_dens_a_the, d_dens_the, & ! due to thermals … … 1414 1414 endif ! 1415 1415 !======================================================================! 1416 #ifdef ISO 1417 #ifdef ISOVERIF 1418 write(*,*) 'physiq 1421b: qx(1,1,:)=',qx(1,1,:) 1419 #endif 1420 #endif 1416 1421 1417 1422 … … 4445 4450 dt_a, dq_a, cv_gen, & 4446 4451 sigd, cin, & 4447 wake_deltat, wake_deltaq, wake_s, awake_ dens,wake_dens, &4452 wake_deltat, wake_deltaq, wake_s, awake_s, wake_dens, awake_dens, & 4448 4453 wake_dth, wake_h, & 4449 4454 !! wake_pe, wake_fip, wake_gfl, & … … 4455 4460 wake_omg, wake_dp_deltomg, & 4456 4461 wake_spread, wake_Cstar, d_deltat_wk_gw, & 4457 d_deltat_wk, d_deltaq_wk, d_s_wk, d_ dens_a_wk, d_dens_wk &4462 d_deltat_wk, d_deltaq_wk, d_s_wk, d_s_a_wk, d_dens_wk, d_dens_a_wk & 4458 4463 #ifdef ISO 4459 4464 & ,xt_seri,dxt_dwn,dxt_a & … … 4495 4500 4496 4501 CALL add_wake_tend & 4497 (d_deltat_wk, d_deltaq_wk, d_s_wk, d_ dens_a_wk, d_dens_wk, wake_k, &4502 (d_deltat_wk, d_deltaq_wk, d_s_wk, d_s_a_wk, d_dens_wk, d_dens_a_wk, wake_k, & 4498 4503 'wake', abortphy & 4499 4504 #ifdef ISO … … 4770 4775 ,zqla,ztva & 4771 4776 #ifdef ISO 4772 & ,xt_ env,d_xt_ajs &4777 & ,xt_therm,d_xt_ajs & 4773 4778 #ifdef DIAGISO 4774 4779 & ,q_the,xt_the & … … 6557 6562 IF(flag_verbose_strataer) print *,'IN physiq_mod: min max d_q_emiss=',& 6558 6563 minval(d_q_emiss),maxval(d_q_emiss) 6564 #ifdef ISO 6565 abort_message='isotopes pas encore dans emission volc H2O' 6566 CALL abort_physic(modname,abort_message,1) 6567 #endif 6559 6568 6560 6569 CALL add_phys_tend(du0, dv0, dt0, d_q_emiss, dql0, dqi0, dqbs0, paprs, & … … 6948 6957 6949 6958 print*,'avt add phys rep',abortphy 6950 6959 #ifdef ISO 6960 abort_message='isotopes pas encore dans rep' 6961 CALL abort_physic(modname,abort_message,1) 6962 #endif 6951 6963 CALL add_phys_tend & 6952 6964 (du0,dv0,dt0,d_q_rep,d_ql_rep,d_qi_rep,dqbs0,paprs,&
Note: See TracChangeset
for help on using the changeset viewer.