Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/lmdz_scm.F90
- Timestamp:
- Aug 3, 2024, 2:56:58 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/lmdz_scm.F90
r5159 r5160 347 347 fnday = -nday / float(day_step) 348 348 endif 349 print*, 'fnday=', fnday349 PRINT *, 'fnday=', fnday 350 350 ! start_time doit etre en FRACTION DE JOUR 351 351 start_time = time_ini / 24. … … 372 372 ! Print out the actual date of the beginning of the simulation : 373 373 CALL ju2ymds(daytime, year_print, month_print, day_print, sec_print) 374 print*, ' Time of beginning : ', &374 PRINT *, ' Time of beginning : ', & 375 375 year_print, month_print, day_print, sec_print 376 376 … … 429 429 fsens = -wtsurf * rcpd * rho(1) 430 430 flat = -wqsurf * rlvtt * rho(1) 431 print*, 'Flux: ok_flux wtsurf wqsurf', ok_flux_surf, wtsurf, wqsurf431 PRINT *, 'Flux: ok_flux wtsurf wqsurf', ok_flux_surf, wtsurf, wqsurf 432 432 ENDIF 433 433 PRINT*, 'Flux sol ', fsens, flat … … 440 440 IF (ok_old_disvert) THEN 441 441 CALL disvert0(pa, preff, ap, bp, dpres, presnivs, nivsigs, nivsig) 442 print*, 'On utilise disvert0'442 PRINT *, 'On utilise disvert0' 443 443 aps(1:llm) = 0.5 * (ap(1:llm) + ap(2:llm + 1)) 444 444 bps(1:llm) = 0.5 * (bp(1:llm) + bp(2:llm + 1)) … … 447 447 ELSE 448 448 CALL disvert() 449 print*, 'On utilise disvert'449 PRINT *, 'On utilise disvert' 450 450 ! Nouvelle version disvert permettant d imposer ap,bp (modif L.Guez) MPL 18092012 451 451 ! Dans ce cas, on lit ap,bp dans le fichier hybrid.txt … … 464 464 IF (trouve_700 .AND. play(l)<=70000) THEN 465 465 llm700 = l 466 print*, 'llm700,play=', llm700, play(l) / 100.466 PRINT *, 'llm700,play=', llm700, play(l) / 100. 467 467 trouve_700 = .FALSE. 468 468 endif … … 750 750 ! Print out the actual date of the beginning of the simulation : 751 751 CALL ju2ymds(daytime, an, mois, jour, heure) 752 print*, ' Time of beginning : y m d h', an, mois,jour, heure/3600.752 PRINT *, ' Time of beginning : y m d h', an, mois,jour, heure/3600. 753 753 754 754 day = int(daytime) … … 844 844 !--------------------------------------------------------------------- 845 845 IF (prt_level>=1) THEN 846 print*, ' avant physiq : -------- day time ', day, time846 PRINT *, ' avant physiq : -------- day time ', day, time 847 847 WRITE(*, *) 'firstcall,lastcall,phis', & 848 848 firstcall, lastcall, phis … … 1019 1019 1020 1020 IF (prt_level>=3) THEN 1021 print*, &1021 PRINT *, & 1022 1022 'physiq-> temp(1),dt_phys(1),d_t_adv(1),dt_cooling(1) ', & 1023 1023 temp(1), dt_phys(1), d_t_adv(1), dt_cooling(1)
Note: See TracChangeset
for help on using the changeset viewer.