Changeset 4295
- Timestamp:
- Oct 7, 2022, 4:06:19 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/wake.F90
r4294 r4295 291 291 REAL, DIMENSION(klon) :: awdens_in, wdens_in ! pour les prints 292 292 293 write(81) klon,klev,znatsurf,p,pi,ph,omgb,dtime,tenv0,qe0,dtdwn,dqdwn,amdwn,amup,dta,dqa,wgen,sigd_con,cin,deltatw,deltaqw,sigmaw,awdens,wdens294 293 ! ------------------------------------------------------------------------- 295 294 ! Initialisations … … 319 318 ! alpk = 0.5 320 319 ! alpk = 0.05 321 print *,'XXXX dtime input ', dtime320 !print *,'XXXX dtime input ', dtime 322 321 igout = klon/2+1/klon 323 322 … … 451 450 d_sig_bnd2(i) = sigmaw_targ - sigmaw(i) 452 451 d_sigmaw2(i) = sigmaw_targ - sigmaw(i) 453 print *,'XXXX1 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)452 ! print *,'XXXX1 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 454 453 sigmaw(i) = sigmaw_targ 455 454 !>jyg … … 836 835 d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i) 837 836 d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i) 838 print *,'XXXX2 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)837 ! print *,'XXXX2 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 839 838 sigmaw(i) = sigmaw_targ 840 839 !>jyg … … 950 949 d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i) 951 950 d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i) 952 print *,'XXXX3 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)951 ! print *,'XXXX3 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 953 952 sigmaw(i) = sigmaw_targ 954 953 !>jyg … … 1496 1495 sigmaw(i) = sigmaw(i) + d_sigmaw(i) 1497 1496 d_sigmaw2(i) = d_sigmaw2(i) + d_sigmaw(i) 1498 print *,'XXXX4 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)1497 ! print *,'XXXX4 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 1499 1498 END IF 1500 1499 END DO … … 1518 1517 d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i) 1519 1518 d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i) 1520 print *,'XXXX5 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)1519 ! print *,'XXXX5 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 1521 1520 sigmaw(i) = sigmaw_targ 1522 1521 END IF … … 1822 1821 d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i) 1823 1822 d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i) 1824 print *,'XXXX6 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)1823 ! print *,'XXXX6 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 1825 1824 sigmaw(i) = sigmaw_targ 1826 1825 !>jyg … … 2041 2040 d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i) 2042 2041 d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i) 2043 print *,'XXXX7 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)2042 ! print *,'XXXX7 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 2044 2043 sigmaw(i) = sigmaw_targ 2045 2044 !>jyg … … 2141 2140 !! d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i) 2142 2141 d_sigmaw2(i) = sigmaw_targ - sigmaw_in(i) ! _in = correction jyg 20220124 2143 print *,'XXXX8 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)2142 ! print *,'XXXX8 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i) 2144 2143 sigmaw(i) = sigmaw_targ 2145 2144 IF (iflag_wk_pop_dyn >= 1) THEN … … 2209 2208 d_sig_bnd2(i) = d_sig_bnd2(i)/dtime 2210 2209 d_sigmaw2(i) = d_sigmaw2(i)/dtime 2211 print *,'XXXX9 d_sigmaw2(i), sigmaw(i), dtime ', d_sigmaw2(i), sigmaw(i), dtime2210 ! print *,'XXXX9 d_sigmaw2(i), sigmaw(i), dtime ', d_sigmaw2(i), sigmaw(i), dtime 2212 2211 ! 2213 2212 d_dens_gen2(i) = d_dens_gen2(i)/dtime
Note: See TracChangeset
for help on using the changeset viewer.