47 |
use netcdf95, only: NF95_CLOSE |
use netcdf95, only: NF95_CLOSE |
48 |
use newmicro_m, only: newmicro |
use newmicro_m, only: newmicro |
49 |
use nr_util, only: assert |
use nr_util, only: assert |
|
use nuage_m, only: nuage |
|
50 |
USE orbite_m, ONLY: orbite |
USE orbite_m, ONLY: orbite |
51 |
USE ozonecm_m, ONLY: ozonecm |
USE ozonecm_m, ONLY: ozonecm |
52 |
USE phyetat0_m, ONLY: phyetat0 |
USE phyetat0_m, ONLY: phyetat0 |
351 |
! Param\`etres li\'es au nouveau sch\'ema de nuages : |
! Param\`etres li\'es au nouveau sch\'ema de nuages : |
352 |
real:: fact_cldcon = 0.375 |
real:: fact_cldcon = 0.375 |
353 |
real:: facttemps = 1.e-4 |
real:: facttemps = 1.e-4 |
|
logical:: ok_newmicro = .true. |
|
354 |
real facteur |
real facteur |
355 |
|
|
356 |
integer:: iflag_cldcon = 1 |
integer:: iflag_cldcon = 1 |
404 |
|
|
405 |
integer, save:: ncid_startphy |
integer, save:: ncid_startphy |
406 |
|
|
407 |
namelist /physiq_nml/ fact_cldcon, facttemps, ok_newmicro, iflag_cldcon, & |
namelist /physiq_nml/ fact_cldcon, facttemps, iflag_cldcon, ratqsbas, & |
408 |
ratqsbas, ratqshaut, ok_ade, bl95_b0, bl95_b1 |
ratqshaut, ok_ade, bl95_b0, bl95_b1 |
409 |
|
|
410 |
!---------------------------------------------------------------- |
!---------------------------------------------------------------- |
411 |
|
|
469 |
ENDIF |
ENDIF |
470 |
|
|
471 |
! Initialisation des sorties |
! Initialisation des sorties |
472 |
call ini_histins(ok_newmicro) |
call ini_histins |
473 |
CALL phyredem0 |
CALL phyredem0 |
474 |
call conf_interface |
call conf_interface |
475 |
ENDIF test_firstcal |
ENDIF test_firstcal |
818 |
|
|
819 |
! Param\`etres optiques des nuages et quelques param\`etres pour |
! Param\`etres optiques des nuages et quelques param\`etres pour |
820 |
! diagnostics : |
! diagnostics : |
821 |
if (ok_newmicro) then |
CALL newmicro(paprs, play, t_seri, cldliq, cldfra, cldtau, cldemi, cldh, & |
822 |
CALL newmicro(paprs, play, t_seri, cldliq, cldfra, cldtau, cldemi, & |
cldl, cldm, cldt, cldq, flwp, fiwp, flwc, fiwc) |
|
cldh, cldl, cldm, cldt, cldq, flwp, fiwp, flwc, fiwc) |
|
|
else |
|
|
CALL nuage(paprs, play, t_seri, cldliq, cldfra, cldtau, cldemi, cldh, & |
|
|
cldl, cldm, cldt, cldq) |
|
|
endif |
|
823 |
|
|
824 |
IF (MOD(itap - 1, radpas) == 0) THEN |
IF (MOD(itap - 1, radpas) == 0) THEN |
825 |
wo = ozonecm(REAL(julien), paprs) |
wo = ozonecm(REAL(julien), paprs) |