Changeset 4088
- Timestamp:
- Feb 27, 2026, 8:46:42 AM (10 days ago)
- Location:
- trunk/LMDZ.VENUS/libf/phyvenus
- Files:
-
- 1 deleted
- 2 edited
-
conf_phys.F90 (modified) (4 diffs)
-
diagphy.F (deleted)
-
physiq_mod.F (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90
r4080 r4088 5 5 contains 6 6 7 subroutine conf_phys(ok_journe, ok_mensuel, ok_instan, & 8 & if_ebil) 7 subroutine conf_phys(ok_journe, ok_mensuel, ok_instan) 9 8 use init_print_control_mod, only: init_print_control 10 9 use print_control_mod, only: lunout … … 25 24 ! Sortie: 26 25 logical,intent(out) :: ok_journe, ok_mensuel, ok_instan 27 integer,intent(out) :: if_ebil28 26 29 27 … … 151 149 ! 152 150 ! 153 !Config Key = if_ebil 154 !Config Desc = Niveau de sortie pour les diags bilan d'energie 155 !Config Def = 0 156 !Config Help = 157 ! 158 ! 159 if_ebil = 0 160 call getin_p('if_ebil', if_ebil) 151 161 152 !! 162 153 !! Parametres orbitaux … … 746 737 write(lunout,*)' Sortie instantanee = ', ok_instan 747 738 write(lunout,*)' frequence sorties = ', ecriphy 748 write(lunout,*)' Sortie bilan d''energie, if_ebil =', if_ebil749 739 write(lunout,*)' Excentricite = ',R_ecc 750 740 write(lunout,*)' Equinoxe = ',R_peri -
trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F
r4087 r4088 372 372 REAL, SAVE :: h_qs_tot, qw_tot, ql_tot, qs_tot , ec_tot 373 373 !$OMP THREADPRIVATE(h_qs_tot,qw_tot,ql_tot,qs_tot,ec_tot) 374 REAL d_h_vcol, d_h_dair, d_qt, d_qw, d_ql, d_qs, d_ec375 REAL, SAVE :: d_h_vcol_phy376 !$OMP THREADPRIVATE(d_h_vcol_phy)377 REAL fs_bound, fq_bound378 REAL zero_v(klon),zero_v2(klon,klev)379 CHARACTER*15 ztit380 INTEGER, SAVE :: ip_ebil ! PRINT level for energy conserv. diag.381 !$OMP THREADPRIVATE(ip_ebil)382 DATA ip_ebil/2/383 INTEGER, SAVE :: if_ebil ! level for energy conserv. dignostics384 !$OMP THREADPRIVATE(if_ebil)385 374 c+jld ec_conser 386 375 REAL d_t_ec(klon,klev) ! tendance du a la conversion Ec -> E thermique … … 416 405 endif 417 406 #endif 418 IF (if_ebil.ge.1) THEN419 DO i=1,klon420 zero_v(i)=0.421 END DO422 DO i=1,klon423 DO j=1,klev424 zero_v2(i,j)=0.425 END DO426 END DO427 END IF428 407 429 408 c====================================================================== … … 446 425 ! CALL suphec ! initialiser constantes et parametres phys. 447 426 448 IF (if_ebil.ge.1) d_h_vcol_phy=0.449 427 ! 450 428 ! load flag and parameter values from physiq.def 451 429 ! 452 430 call conf_phys(ok_journe, ok_mensuel, 453 . ok_instan, 454 . if_ebil) 431 & ok_instan) 455 432 456 433 #ifndef MESOSCALE … … 811 788 ENDDO 812 789 C 813 IF (if_ebil.ge.1) THEN814 ztit='after dynamic'815 CALL diagetpq(cell_area,ztit,ip_ebil,1,1,dtime816 e , t_seri,zero_v2,zero_v2,zero_v2,u_seri,v_seri,paprs,pplay817 s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)818 C Comme les tendances de la physique sont ajoute dans la dynamique,819 C on devrait avoir que la variation d'entalpie par la dynamique820 C est egale a la variation de la physique au pas de temps precedent.821 C Donc la somme de ces 2 variations devrait etre nulle.822 call diagphy(cell_area,ztit,ip_ebil823 e , zero_v, zero_v, zero_v, zero_v, zero_v824 e , zero_v, zero_v, zero_v, ztsol825 e , d_h_vcol+d_h_vcol_phy, d_qt, 0.826 s , fs_bound, fq_bound )827 END IF828 790 829 791 c==================================================================== … … 1441 1403 endif 1442 1404 1443 IF (if_ebil.ge.2) THEN1444 ztit='after clmain'1445 CALL diagetpq(cell_area,ztit,ip_ebil,2,1,dtime1446 e , t_seri,zero_v2,zero_v2,zero_v2,u_seri,v_seri,paprs,pplay1447 s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)1448 call diagphy(cell_area,ztit,ip_ebil1449 e , zero_v, zero_v, zero_v, zero_v, sens1450 e , zero_v, zero_v, zero_v, ztsol1451 e , d_h_vcol, d_qt, d_ec1452 s , fs_bound, fq_bound )1453 END IF1454 1405 C 1455 1406 c … … 1525 1476 endif 1526 1477 endif 1527 1528 IF (if_ebil.ge.2) THEN1529 ztit='after dry_adjust'1530 CALL diagetpq(cell_area,ztit,ip_ebil,2,2,dtime1531 e , t_seri,zero_v2,zero_v2,zero_v2,u_seri,v_seri,paprs,pplay1532 s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)1533 call diagphy(cell_area,ztit,ip_ebil1534 e , zero_v, zero_v, zero_v, zero_v, sens1535 e , zero_v, zero_v, zero_v, ztsol1536 e , d_h_vcol, d_qt, d_ec1537 s , fs_bound, fq_bound )1538 END IF1539 1478 1540 1479 c==================================================================== … … 1804 1743 ! call writefield_phy('physiq_dtrad',dtrad,klev) 1805 1744 1806 IF (if_ebil.ge.2) THEN1807 ztit='after rad'1808 CALL diagetpq(cell_area,ztit,ip_ebil,2,2,dtime1809 e , t_seri,zero_v2,zero_v2,zero_v2,u_seri,v_seri,paprs,pplay1810 s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)1811 call diagphy(cell_area,ztit,ip_ebil1812 e , topsw, toplw, solsw, sollw, zero_v1813 e , zero_v, zero_v, zero_v, ztsol1814 e , d_h_vcol, d_qt, d_ec1815 s , fs_bound, fq_bound )1816 END IF1817 c1818 1819 1745 c==================================================================== 1820 1746 c Calcul des gravity waves FLOTT … … 2063 1989 enddo 2064 1990 2065 c-jld ec_conser2066 c====================================================================2067 IF (if_ebil.ge.1) THEN2068 ztit='after physic'2069 CALL diagetpq(cell_area,ztit,ip_ebil,1,1,dtime2070 e , t_seri,zero_v2,zero_v2,zero_v2,u_seri,v_seri,paprs,pplay2071 s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)2072 C Comme les tendances de la physique sont ajoute dans la dynamique,2073 C on devrait avoir que la variation d'entalpie par la dynamique2074 C est egale a la variation de la physique au pas de temps precedent.2075 C Donc la somme de ces 2 variations devrait etre nulle.2076 call diagphy(cell_area,ztit,ip_ebil2077 e , topsw, toplw, solsw, sollw, sens2078 e , zero_v, zero_v, zero_v, ztsol2079 e , d_h_vcol, d_qt, d_ec2080 s , fs_bound, fq_bound )2081 C2082 d_h_vcol_phy=d_h_vcol2083 C2084 END IF2085 1991 C 2086 1992 c=======================================================================
Note: See TracChangeset
for help on using the changeset viewer.
