- Timestamp:
- Jul 22, 2024, 9:29:09 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/add_phys_tend_mod.F90
r5087 r5099 1 ! 1 2 2 ! $Id$ 3 ! 4 ! 3 4 5 5 MODULE add_phys_tend_mod 6 6 … … 120 120 RETURN 121 121 END SUBROUTINE add_pbl_tend 122 ! 122 123 123 ! $Id$ 124 ! 124 125 125 SUBROUTINE add_phys_tend (zdu,zdv,zdt,zdq,zdql,zdqi,zdqbs,paprs,text, & 126 126 abortphy,flag_inhib_tend, itap, diag_mode & … … 175 175 INTEGER, INTENT(IN) :: diag_mode ! 0 -> normal effective mode 176 176 ! 1 -> only conservation stats are computed 177 ! 177 178 178 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: zdq 179 179 #ifdef ISO … … 204 204 REAL, DIMENSION(ntraciso,klon,klev) :: sav_xtl_seri, sav_xts_seri, sav_xt_seri 205 205 #endif 206 ! 206 207 207 INTEGER i, k,j, n 208 208 INTEGER jadrs(klon*klev), jbad … … 216 216 logical, save :: first=.true. 217 217 !$OMP THREADPRIVATE(first) 218 ! 218 219 219 !====================================================================== 220 220 ! Variables for energy conservation tests 221 221 !====================================================================== 222 !223 222 224 223 ! zh_col------- total enthalpy of vertical air column … … 234 233 ! zqbs_col------ total mass of blowing snow (kg/m2) 235 234 ! zek_col------ total kinetic energy (kg/m2) 236 ! 235 237 236 REAL zairm(klon, klev) ! layer air mass (kg/m2) 238 237 REAL zqw_col(klon,2) … … 273 272 first=.false. 274 273 endif 275 ! 274 276 275 ! print *,'add_phys_tend: paprs ',paprs 277 276 ! When in diagnostic mode, save initial values of out variables … … 422 421 ENDDO 423 422 ENDIF 424 ! 423 425 424 !===================================================================================== 426 425 ! Impression, warning et correction en cas de probleme moins important … … 531 530 #endif 532 531 #endif 533 !534 532 535 533 !IM ajout memes tests pour reverifier les jbad, jqbad beg … … 572 570 ENDDO 573 571 ENDIF 574 ! 572 575 573 IF (jqbad > 0) THEN 576 574 DO j = 1, jqbad … … 646 644 647 645 end if ! end if (fl_ebil .GT. 0) 648 ! 646 649 647 ! When in diagnostic mode, restore "out" variables to initial values. 650 648 IF (diag_mode == 1) THEN … … 706 704 REAL, DIMENSION(nlon,nlev) :: temp_n, qv_n, ql_n, qs_n, qbs_n 707 705 708 709 !710 706 INTEGER k, n 711 707 … … 713 709 logical, save :: first=.true. 714 710 !$OMP THREADPRIVATE(first) 715 ! 711 716 712 !====================================================================== 717 713 ! Variables for energy conservation tests 718 714 !====================================================================== 719 !720 715 721 716 ! zh_col------- total enthalpy of vertical air column … … 730 725 ! zqbs_col------ total mass of blowing snow (kg/m2) 731 726 ! zek_col------ total kinetic energy (kg/m2) 732 ! 727 733 728 REAL zairm(nlon, nlev) ! layer air mass (kg/m2) 734 729 REAL zqw_col(nlon,2) … … 754 749 first=.false. 755 750 endif 756 ! 751 757 752 ! print *,'add_phys_tend: paprs ',paprs 758 753 !====================================================================== … … 832 827 833 828 end if ! end if (fl_ebil .GT. 0) 834 !835 829 836 830 RETURN … … 978 972 if ( abs(bilq_error) > bilq_seuil) bilq_ok=1 979 973 if ( abs(bilh_error) > bilh_seuil) bilh_ok=1 980 ! 974 981 975 ! Print diagnostics 982 976 ! =================
Note: See TracChangeset
for help on using the changeset viewer.