Changeset 5256 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
- Timestamp:
- Oct 22, 2024, 4:06:28 PM (4 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r5255 r5256 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_s_a_wk,d_dens_a_wk, d_dens_wk, & ! due to wakes331 d_s_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 ISO1417 #ifdef ISOVERIF1418 write(*,*) 'physiq 1421b: qx(1,1,:)=',qx(1,1,:)1419 #endif1420 #endif1421 1416 1422 1417 … … 4450 4445 dt_a, dq_a, cv_gen, & 4451 4446 sigd, cin, & 4452 wake_deltat, wake_deltaq, wake_s, awake_ s, wake_dens, awake_dens, &4447 wake_deltat, wake_deltaq, wake_s, awake_dens, wake_dens, & 4453 4448 wake_dth, wake_h, & 4454 4449 !! wake_pe, wake_fip, wake_gfl, & … … 4460 4455 wake_omg, wake_dp_deltomg, & 4461 4456 wake_spread, wake_Cstar, d_deltat_wk_gw, & 4462 d_deltat_wk, d_deltaq_wk, d_s_wk, d_ s_a_wk, d_dens_wk, d_dens_a_wk &4457 d_deltat_wk, d_deltaq_wk, d_s_wk, d_dens_a_wk, d_dens_wk & 4463 4458 #ifdef ISO 4464 4459 & ,xt_seri,dxt_dwn,dxt_a & … … 4500 4495 4501 4496 CALL add_wake_tend & 4502 (d_deltat_wk, d_deltaq_wk, d_s_wk, d_ s_a_wk, d_dens_wk, d_dens_a_wk, wake_k, &4497 (d_deltat_wk, d_deltaq_wk, d_s_wk, d_dens_a_wk, d_dens_wk, wake_k, & 4503 4498 'wake', abortphy & 4504 4499 #ifdef ISO … … 4775 4770 ,zqla,ztva & 4776 4771 #ifdef ISO 4777 & ,xt_ therm,d_xt_ajs &4772 & ,xt_env,d_xt_ajs & 4778 4773 #ifdef DIAGISO 4779 4774 & ,q_the,xt_the & … … 6562 6557 IF(flag_verbose_strataer) print *,'IN physiq_mod: min max d_q_emiss=',& 6563 6558 minval(d_q_emiss),maxval(d_q_emiss) 6564 #ifdef ISO6565 abort_message='isotopes pas encore dans emission volc H2O'6566 CALL abort_physic(modname,abort_message,1)6567 #endif6568 6559 6569 6560 CALL add_phys_tend(du0, dv0, dt0, d_q_emiss, dql0, dqi0, dqbs0, paprs, & … … 6957 6948 6958 6949 print*,'avt add phys rep',abortphy 6959 #ifdef ISO 6960 abort_message='isotopes pas encore dans rep' 6961 CALL abort_physic(modname,abort_message,1) 6962 #endif 6950 6963 6951 CALL add_phys_tend & 6964 6952 (du0,dv0,dt0,d_q_rep,d_ql_rep,d_qi_rep,dqbs0,paprs,&
Note: See TracChangeset
for help on using the changeset viewer.