- Timestamp:
- Jul 20, 2024, 12:00:23 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/add_phys_tend_mod.F90
r5081 r5087 151 151 USE cmp_seri_mod 152 152 USE phys_output_var_mod, ONLY : d_qw_col, d_ql_col, d_qs_col, d_qbs_col, d_qt_col, d_ek_col, d_h_dair_col & 153 &, d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col153 , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col 154 154 155 155 #ifdef ISO … … 330 330 #ifdef ISOTRAC 331 331 call iso_verif_traceur_vect(xt_seri,klon,klev, & 332 &'add_phys_tend 328a: '//text)332 'add_phys_tend 328a: '//text) 333 333 call iso_verif_traceur_vect(xtl_seri,klon,klev, & 334 &'add_phys_tend 328b: '//text)334 'add_phys_tend 328b: '//text) 335 335 #endif 336 336 #endif … … 396 396 #ifdef ISOTRAC 397 397 call iso_verif_traceur_vect(xt_seri,klon,klev, & 398 &'add_phys_tend 374a: '//text)398 'add_phys_tend 374a: '//text) 399 399 call iso_verif_traceur_vect(xtl_seri,klon,klev, & 400 &'add_phys_tend 374b: '//text)400 'add_phys_tend 374b: '//text) 401 401 #endif 402 402 #endif … … 527 527 #ifdef ISOTRAC 528 528 call iso_verif_traceur_vect(xt_seri,klon,klev, & 529 &'add_phys_tend 499'//text)529 'add_phys_tend 499'//text) 530 530 #endif 531 531 #endif … … 562 562 print*,'PLANTAGE2 POUR LE POINT i itap lon lat txt jbad zdt t',& 563 563 i,itap,longitude_deg(i),latitude_deg(i),text,jbad, & 564 &zdt(i,k),t_seri(i,k)-zdt(i,k)564 zdt(i,k),t_seri(i,k)-zdt(i,k) 565 565 !!! if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN 566 566 print*,'l T dT Q dQ ' … … 580 580 print*,'WARNING : EAU2 POUR LE POINT i itap lon lat txt jqbad zdq q zdql ql',& 581 581 i,itap,longitude_deg(i),latitude_deg(i),text,jqbad,& 582 &zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)582 zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k) 583 583 !!! if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN 584 584 print*,'l T dT Q dQ ' … … 686 686 USE cmp_seri_mod 687 687 USE phys_output_var_mod, ONLY : d_qw_col, d_ql_col, d_qs_col, d_qbs_col, d_qt_col, d_ek_col, d_h_dair_col & 688 &, d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col688 , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col 689 689 IMPLICIT none 690 690 include "YOMCST.h" … … 918 918 USE cmp_seri_mod 919 919 USE phys_output_var_mod, ONLY : d_qw_col, d_ql_col, d_qs_col, d_qbs_col, d_qt_col, d_ek_col, d_h_dair_col & 920 &, d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col920 , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col 921 921 USE phys_local_var_mod, ONLY: evap, sens 922 922 USE phys_local_var_mod, ONLY: u_seri, v_seri, ql_seri, qs_seri, qbs_seri, q_seri, t_seri & 923 &, rain_lsc, snow_lsc923 , rain_lsc, snow_lsc 924 924 USE climb_hq_mod, ONLY : d_h_col_vdf, f_h_bnd 925 925 IMPLICIT none … … 956 956 bilq_bnd = - rain_lsc(1) - snow_lsc(1) 957 957 bilh_bnd = (-(rcw-rcpd)*t_seri(1,1) + rlvtt) * rain_lsc(1) & 958 &+ (-(rcs-rcpd)*t_seri(1,1) + rlstt) * snow_lsc(1)958 + (-(rcs-rcpd)*t_seri(1,1) + rlstt) * snow_lsc(1) 959 959 CASE("bsss") param 960 960 bilq_bnd = - bs_fall(1) … … 963 963 bilq_bnd = - rain_con(1) - snow_con(1) 964 964 bilh_bnd = (-(rcw-rcpd)*t_seri(1,1) + rlvtt) * rain_con(1) & 965 &+ (-(rcs-rcpd)*t_seri(1,1) + rlstt) * snow_con(1)965 + (-(rcs-rcpd)*t_seri(1,1) + rlstt) * snow_con(1) 966 966 CASE("SW") param 967 967 bilh_bnd = topsw(1) - solsw(1)
Note: See TracChangeset
for help on using the changeset viewer.