/[lmdze]/trunk/phylmd/physiq.f
ViewVC logotype

Diff of /trunk/phylmd/physiq.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 190 by guez, Thu Apr 14 15:15:56 2016 UTC revision 202 by guez, Wed Jun 8 12:23:41 2016 UTC
# Line 16  contains Line 16  contains
16    
17      use aaam_bud_m, only: aaam_bud      use aaam_bud_m, only: aaam_bud
18      USE abort_gcm_m, ONLY: abort_gcm      USE abort_gcm_m, ONLY: abort_gcm
     use aeropt_m, only: aeropt  
19      use ajsec_m, only: ajsec      use ajsec_m, only: ajsec
20      use calltherm_m, only: calltherm      use calltherm_m, only: calltherm
21      USE clesphys, ONLY: cdhmax, cdmmax, ecrit_hf, ecrit_ins, ecrit_mth, &      USE clesphys, ONLY: cdhmax, cdmmax, ecrit_ins, ksta, ksta_ter, ok_kzmin, &
22           ecrit_reg, ecrit_tra, ksta, ksta_ter, ok_kzmin           ok_instan
23      USE clesphys2, ONLY: cycle_diurne, conv_emanuel, nbapp_rad, new_oliq, &      USE clesphys2, ONLY: cycle_diurne, conv_emanuel, nbapp_rad, new_oliq, &
24           ok_orodr, ok_orolf           ok_orodr, ok_orolf
25      USE clmain_m, ONLY: clmain      USE clmain_m, ONLY: clmain
# Line 28  contains Line 27  contains
27      use comconst, only: dtphys      use comconst, only: dtphys
28      USE comgeomphy, ONLY: airephy      USE comgeomphy, ONLY: airephy
29      USE concvl_m, ONLY: concvl      USE concvl_m, ONLY: concvl
30      USE conf_gcm_m, ONLY: offline, raz_date, day_step, iphysiq      USE conf_gcm_m, ONLY: offline, day_step, iphysiq, lmt_pas
31      USE conf_phys_m, ONLY: conf_phys      USE conf_phys_m, ONLY: conf_phys
32      use conflx_m, only: conflx      use conflx_m, only: conflx
33      USE ctherm, ONLY: iflag_thermals, nsplit_thermals      USE ctherm, ONLY: iflag_thermals, nsplit_thermals
# Line 43  contains Line 42  contains
42      USE fcttre, ONLY: foeew, qsatl, qsats, thermcep      USE fcttre, ONLY: foeew, qsatl, qsats, thermcep
43      use fisrtilp_m, only: fisrtilp      use fisrtilp_m, only: fisrtilp
44      USE hgardfou_m, ONLY: hgardfou      USE hgardfou_m, ONLY: hgardfou
45        USE histsync_m, ONLY: histsync
46        USE histwrite_phy_m, ONLY: histwrite_phy
47      USE indicesol, ONLY: clnsurf, epsfra, is_lic, is_oce, is_sic, is_ter, &      USE indicesol, ONLY: clnsurf, epsfra, is_lic, is_oce, is_sic, is_ter, &
48           nbsrf           nbsrf
49      USE ini_histins_m, ONLY: ini_histins      USE ini_histins_m, ONLY: ini_histins, nid_ins
50      use netcdf95, only: NF95_CLOSE      use netcdf95, only: NF95_CLOSE
51      use newmicro_m, only: newmicro      use newmicro_m, only: newmicro
52        use nr_util, only: assert
53      use nuage_m, only: nuage      use nuage_m, only: nuage
54      USE orbite_m, ONLY: orbite      USE orbite_m, ONLY: orbite
55      USE ozonecm_m, ONLY: ozonecm      USE ozonecm_m, ONLY: ozonecm
# Line 58  contains Line 60  contains
60      USE phytrac_m, ONLY: phytrac      USE phytrac_m, ONLY: phytrac
61      USE qcheck_m, ONLY: qcheck      USE qcheck_m, ONLY: qcheck
62      use radlwsw_m, only: radlwsw      use radlwsw_m, only: radlwsw
     use readsulfate_m, only: readsulfate  
     use readsulfate_preind_m, only: readsulfate_preind  
63      use yoegwd, only: sugwd      use yoegwd, only: sugwd
64      USE suphec_m, ONLY: rcpd, retv, rg, rlvtt, romega, rsigma, rtt      USE suphec_m, ONLY: rcpd, retv, rg, rlvtt, romega, rsigma, rtt
65        use time_phylmdz, only: itap, increment_itap
66      use transp_m, only: transp      use transp_m, only: transp
67      use transp_lay_m, only: transp_lay      use transp_lay_m, only: transp_lay
68      use unit_nml_m, only: unit_nml      use unit_nml_m, only: unit_nml
# Line 88  contains Line 89  contains
89      REAL, intent(in):: pphis(:) ! (klon) géopotentiel du sol      REAL, intent(in):: pphis(:) ! (klon) géopotentiel du sol
90    
91      REAL, intent(in):: u(:, :) ! (klon, llm)      REAL, intent(in):: u(:, :) ! (klon, llm)
92      ! vitesse dans la direction X (de O a E) en m/s      ! vitesse dans la direction X (de O a E) en m / s
93    
94      REAL, intent(in):: v(:, :) ! (klon, llm) vitesse Y (de S a N) en m/s      REAL, intent(in):: v(:, :) ! (klon, llm) vitesse Y (de S a N) en m / s
95      REAL, intent(in):: t(:, :) ! (klon, llm) temperature (K)      REAL, intent(in):: t(:, :) ! (klon, llm) temperature (K)
96    
97      REAL, intent(in):: qx(:, :, :) ! (klon, llm, nqmx)      REAL, intent(in):: qx(:, :, :) ! (klon, llm, nqmx)
98      ! (humidit\'e sp\'ecifique et fractions massiques des autres traceurs)      ! (humidit\'e sp\'ecifique et fractions massiques des autres traceurs)
99    
100      REAL, intent(in):: omega(:, :) ! (klon, llm) vitesse verticale en Pa/s      REAL, intent(in):: omega(:, :) ! (klon, llm) vitesse verticale en Pa / s
101      REAL, intent(out):: d_u(:, :) ! (klon, llm) tendance physique de "u" (m s-2)      REAL, intent(out):: d_u(:, :) ! (klon, llm) tendance physique de "u" (m s-2)
102      REAL, intent(out):: d_v(:, :) ! (klon, llm) tendance physique de "v" (m s-2)      REAL, intent(out):: d_v(:, :) ! (klon, llm) tendance physique de "v" (m s-2)
103      REAL, intent(out):: d_t(:, :) ! (klon, llm) tendance physique de "t" (K/s)      REAL, intent(out):: d_t(:, :) ! (klon, llm) tendance physique de "t" (K / s)
104    
105      REAL, intent(out):: d_qx(:, :, :) ! (klon, llm, nqmx)      REAL, intent(out):: d_qx(:, :, :) ! (klon, llm, nqmx)
106      ! tendance physique de "qx" (s-1)      ! tendance physique de "qx" (s-1)
# Line 114  contains Line 115  contains
115      LOGICAL, PARAMETER:: ok_stratus = .FALSE.      LOGICAL, PARAMETER:: ok_stratus = .FALSE.
116      ! Ajouter artificiellement les stratus      ! Ajouter artificiellement les stratus
117    
118      logical:: ok_journe = .false., ok_mensuel = .true., ok_instan = .false.      ! pour phystoke avec thermiques
     ! sorties journalieres, mensuelles et instantanees dans les  
     ! fichiers histday, histmth et histins  
   
     LOGICAL ok_region ! sortir le fichier regional  
     PARAMETER (ok_region = .FALSE.)  
   
     ! pour phsystoke avec thermiques  
119      REAL fm_therm(klon, llm + 1)      REAL fm_therm(klon, llm + 1)
120      REAL entr_therm(klon, llm)      REAL entr_therm(klon, llm)
121      real, save:: q2(klon, llm + 1, nbsrf)      real, save:: q2(klon, llm + 1, nbsrf)
# Line 132  contains Line 126  contains
126      REAL, save:: t_ancien(klon, llm), q_ancien(klon, llm)      REAL, save:: t_ancien(klon, llm), q_ancien(klon, llm)
127      LOGICAL, save:: ancien_ok      LOGICAL, save:: ancien_ok
128    
129      REAL d_t_dyn(klon, llm) ! tendance dynamique pour "t" (K/s)      REAL d_t_dyn(klon, llm) ! tendance dynamique pour "t" (K / s)
130      REAL d_q_dyn(klon, llm) ! tendance dynamique pour "q" (kg/kg/s)      REAL d_q_dyn(klon, llm) ! tendance dynamique pour "q" (kg / kg / s)
131    
132      real da(klon, llm), phi(klon, llm, llm), mp(klon, llm)      real da(klon, llm), phi(klon, llm, llm), mp(klon, llm)
133    
# Line 148  contains Line 142  contains
142      ! prw: precipitable water      ! prw: precipitable water
143      real prw(klon)      real prw(klon)
144    
145      ! flwp, fiwp = Liquid Water Path & Ice Water Path (kg/m2)      ! flwp, fiwp = Liquid Water Path & Ice Water Path (kg / m2)
146      ! flwc, fiwc = Liquid Water Content & Ice Water Content (kg/kg)      ! flwc, fiwc = Liquid Water Content & Ice Water Content (kg / kg)
147      REAL flwp(klon), fiwp(klon)      REAL flwp(klon), fiwp(klon)
148      REAL flwc(klon, llm), fiwc(klon, llm)      REAL flwc(klon, llm), fiwc(klon, llm)
149    
# Line 162  contains Line 156  contains
156      REAL radsol(klon)      REAL radsol(klon)
157      SAVE radsol ! bilan radiatif au sol calcule par code radiatif      SAVE radsol ! bilan radiatif au sol calcule par code radiatif
158    
     INTEGER:: itap = 0 ! number of calls to "physiq"  
   
159      REAL, save:: ftsol(klon, nbsrf) ! skin temperature of surface fraction      REAL, save:: ftsol(klon, nbsrf) ! skin temperature of surface fraction
160    
161      REAL, save:: ftsoil(klon, nsoilmx, nbsrf)      REAL, save:: ftsoil(klon, nsoilmx, nbsrf)
# Line 210  contains Line 202  contains
202      REAL ycoefh(klon, llm) ! coef d'echange pour phytrac      REAL ycoefh(klon, llm) ! coef d'echange pour phytrac
203      REAL yu1(klon) ! vents dans la premiere couche U      REAL yu1(klon) ! vents dans la premiere couche U
204      REAL yv1(klon) ! vents dans la premiere couche V      REAL yv1(klon) ! vents dans la premiere couche V
205      REAL ffonte(klon, nbsrf) !Flux thermique utilise pour fondre la neige      REAL ffonte(klon, nbsrf) ! flux thermique utilise pour fondre la neige
206      REAL fqcalving(klon, nbsrf) !Flux d'eau "perdue" par la surface  
207      ! !et necessaire pour limiter la      REAL fqcalving(klon, nbsrf)
208      ! !hauteur de neige, en kg/m2/s      ! flux d'eau "perdue" par la surface et necessaire pour limiter la
209        ! hauteur de neige, en kg / m2 / s
210    
211      REAL zxffonte(klon), zxfqcalving(klon)      REAL zxffonte(klon), zxfqcalving(klon)
212    
213      REAL pfrac_impa(klon, llm)! Produits des coefs lessivage impaction      REAL pfrac_impa(klon, llm)! Produits des coefs lessivage impaction
# Line 226  contains Line 220  contains
220      REAL frac_nucl(klon, llm) ! idem (nucleation)      REAL frac_nucl(klon, llm) ! idem (nucleation)
221    
222      REAL, save:: rain_fall(klon)      REAL, save:: rain_fall(klon)
223      ! liquid water mass flux (kg/m2/s), positive down      ! liquid water mass flux (kg / m2 / s), positive down
224    
225      REAL, save:: snow_fall(klon)      REAL, save:: snow_fall(klon)
226      ! solid water mass flux (kg/m2/s), positive down      ! solid water mass flux (kg / m2 / s), positive down
227    
228      REAL rain_tiedtke(klon), snow_tiedtke(klon)      REAL rain_tiedtke(klon), snow_tiedtke(klon)
229    
# Line 250  contains Line 244  contains
244      ! Conditions aux limites      ! Conditions aux limites
245    
246      INTEGER julien      INTEGER julien
     INTEGER, SAVE:: lmt_pas ! number of time steps of "physics" per day  
247      REAL, save:: pctsrf(klon, nbsrf) ! percentage of surface      REAL, save:: pctsrf(klon, nbsrf) ! percentage of surface
     REAL pctsrf_new(klon, nbsrf) ! pourcentage surfaces issus d'ORCHIDEE  
248      REAL, save:: albsol(klon) ! albedo du sol total visible      REAL, save:: albsol(klon) ! albedo du sol total visible
249      REAL, SAVE:: wo(klon, llm) ! column density of ozone in a cell, in kDU      REAL, SAVE:: wo(klon, llm) ! column density of ozone in a cell, in kDU
250    
# Line 288  contains Line 280  contains
280      real, save:: sollwdown(klon) ! downward LW flux at surface      real, save:: sollwdown(klon) ! downward LW flux at surface
281      REAL, save:: topsw0(klon), toplw0(klon), solsw0(klon), sollw0(klon)      REAL, save:: topsw0(klon), toplw0(klon), solsw0(klon), sollw0(klon)
282      REAL, save:: albpla(klon)      REAL, save:: albpla(klon)
283      REAL fsollw(klon, nbsrf) ! bilan flux IR pour chaque sous surface      REAL fsollw(klon, nbsrf) ! bilan flux IR pour chaque sous-surface
284      REAL fsolsw(klon, nbsrf) ! flux solaire absorb. pour chaque sous surface      REAL fsolsw(klon, nbsrf) ! flux solaire absorb\'e pour chaque sous-surface
285    
286      REAL conv_q(klon, llm) ! convergence de l'humidite (kg/kg/s)      REAL conv_q(klon, llm) ! convergence de l'humidite (kg / kg / s)
287      REAL conv_t(klon, llm) ! convergence of temperature (K/s)      REAL conv_t(klon, llm) ! convergence of temperature (K / s)
288    
289      REAL cldl(klon), cldm(klon), cldh(klon) !nuages bas, moyen et haut      REAL cldl(klon), cldm(klon), cldh(klon) ! nuages bas, moyen et haut
290      REAL cldt(klon), cldq(klon) !nuage total, eau liquide integree      REAL cldt(klon), cldq(klon) ! nuage total, eau liquide integree
291    
292      REAL zxtsol(klon), zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)      REAL zxtsol(klon), zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)
293    
# Line 309  contains Line 301  contains
301      REAL, PARAMETER:: t_coup = 234.      REAL, PARAMETER:: t_coup = 234.
302      REAL zphi(klon, llm)      REAL zphi(klon, llm)
303    
304      ! cf. Anne Mathieu variables pour la couche limite atmosphérique (hbtm)      ! cf. Anne Mathieu, variables pour la couche limite atmosphérique (hbtm)
305    
306      REAL, SAVE:: pblh(klon, nbsrf) ! Hauteur de couche limite      REAL, SAVE:: pblh(klon, nbsrf) ! Hauteur de couche limite
307      REAL, SAVE:: plcl(klon, nbsrf) ! Niveau de condensation de la CLA      REAL, SAVE:: plcl(klon, nbsrf) ! Niveau de condensation de la CLA
# Line 401  contains Line 393  contains
393      REAL zustrph(klon), zvstrph(klon)      REAL zustrph(klon), zvstrph(klon)
394      REAL aam, torsfc      REAL aam, torsfc
395    
     INTEGER, SAVE:: nid_ins  
   
396      REAL ve_lay(klon, llm) ! transport meri. de l'energie a chaque niveau vert.      REAL ve_lay(klon, llm) ! transport meri. de l'energie a chaque niveau vert.
397      REAL vq_lay(klon, llm) ! transport meri. de l'eau a chaque niveau vert.      REAL vq_lay(klon, llm) ! transport meri. de l'eau a chaque niveau vert.
398      REAL ue_lay(klon, llm) ! transport zonal de l'energie a chaque niveau vert.      REAL ue_lay(klon, llm) ! transport zonal de l'energie a chaque niveau vert.
# Line 419  contains Line 409  contains
409      INTEGER:: ip_ebil = 0 ! print level for energy conservation diagnostics      INTEGER:: ip_ebil = 0 ! print level for energy conservation diagnostics
410      INTEGER:: if_ebil = 0 ! verbosity for diagnostics of energy conservation      INTEGER:: if_ebil = 0 ! verbosity for diagnostics of energy conservation
411    
412      REAL d_t_ec(klon, llm) ! tendance due \`a la conversion Ec -> E thermique      REAL d_t_ec(klon, llm)
413        ! tendance due \`a la conversion Ec en énergie thermique
414    
415      REAL ZRCPD      REAL ZRCPD
416    
417      REAL t2m(klon, nbsrf), q2m(klon, nbsrf) ! temperature and humidity at 2 m      REAL t2m(klon, nbsrf), q2m(klon, nbsrf) ! temperature and humidity at 2 m
# Line 429  contains Line 421  contains
421    
422      ! Aerosol effects:      ! Aerosol effects:
423    
424      REAL sulfate(klon, llm) ! SO4 aerosol concentration (micro g/m3)      REAL sulfate(klon, llm) ! SO4 aerosol concentration (micro g / m3)
425    
426      REAL, save:: sulfate_pi(klon, llm)      REAL, save:: sulfate_pi(klon, llm)
427      ! SO4 aerosol concentration, in \mu g/m3, pre-industrial value      ! SO4 aerosol concentration, in \mu g / m3, pre-industrial value
428    
429      REAL cldtaupi(klon, llm)      REAL cldtaupi(klon, llm)
430      ! cloud optical thickness for pre-industrial (pi) aerosols      ! cloud optical thickness for pre-industrial aerosols
431    
432      REAL re(klon, llm) ! Cloud droplet effective radius      REAL re(klon, llm) ! Cloud droplet effective radius
433      REAL fl(klon, llm) ! denominator of re      REAL fl(klon, llm) ! denominator of re
# Line 447  contains Line 439  contains
439      REAL topswad(klon), solswad(klon) ! aerosol direct effect      REAL topswad(klon), solswad(klon) ! aerosol direct effect
440      REAL topswai(klon), solswai(klon) ! aerosol indirect effect      REAL topswai(klon), solswai(klon) ! aerosol indirect effect
441    
     REAL aerindex(klon) ! POLDER aerosol index  
   
442      LOGICAL:: ok_ade = .false. ! apply aerosol direct effect      LOGICAL:: ok_ade = .false. ! apply aerosol direct effect
443      LOGICAL:: ok_aie = .false. ! apply aerosol indirect effect      LOGICAL:: ok_aie = .false. ! apply aerosol indirect effect
444    
# Line 474  contains Line 464  contains
464      real zmasse(klon, llm)      real zmasse(klon, llm)
465      ! (column-density of mass of air in a cell, in kg m-2)      ! (column-density of mass of air in a cell, in kg m-2)
466    
467      integer, save:: ncid_startphy, itau_phy      integer, save:: ncid_startphy
468    
469      namelist /physiq_nml/ ok_journe, ok_mensuel, ok_instan, fact_cldcon, &      namelist /physiq_nml/ fact_cldcon, facttemps, ok_newmicro, &
470           facttemps, ok_newmicro, iflag_cldcon, ratqsbas, ratqshaut, if_ebil, &           iflag_cldcon, ratqsbas, ratqshaut, if_ebil, ok_ade, ok_aie, bl95_b0, &
471           ok_ade, ok_aie, bl95_b0, bl95_b1, iflag_thermals, nsplit_thermals           bl95_b1, iflag_thermals, nsplit_thermals
472    
473      !----------------------------------------------------------------      !----------------------------------------------------------------
474    
# Line 535  contains Line 525  contains
525         ! Initialiser les compteurs:         ! Initialiser les compteurs:
526    
527         frugs = 0.         frugs = 0.
528         CALL phyetat0(pctsrf, ftsol, ftsoil, fqsurf, qsol, &         CALL phyetat0(pctsrf, ftsol, ftsoil, fqsurf, qsol, fsnow, falbe, &
529              fsnow, falbe, fevap, rain_fall, snow_fall, solsw, sollw, dlw, &              fevap, rain_fall, snow_fall, solsw, sollw, dlw, radsol, frugs, &
530              radsol, frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, &              agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &
531              t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, clwcon, &              q_ancien, ancien_ok, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, &
532              run_off_lic_0, sig1, w01, ncid_startphy, itau_phy)              w01, ncid_startphy)
533    
534         ! ATTENTION : il faudra a terme relire q2 dans l'etat initial         ! ATTENTION : il faudra a terme relire q2 dans l'etat initial
535         q2 = 1e-8         q2 = 1e-8
536    
        lmt_pas = day_step / iphysiq  
        print *, 'Number of time steps of "physics" per day: ', lmt_pas  
   
537         radpas = lmt_pas / nbapp_rad         radpas = lmt_pas / nbapp_rad
538           print *, "radpas = ", radpas
        ! On remet le calendrier a zero  
        IF (raz_date) itau_phy = 0  
   
        CALL printflag(radpas, ok_journe, ok_instan, ok_region)  
539    
540         ! Initialisation pour le sch\'ema de convection d'Emanuel :         ! Initialisation pour le sch\'ema de convection d'Emanuel :
541         IF (conv_emanuel) THEN         IF (conv_emanuel) THEN
# Line 567  contains Line 550  contains
550            rugoro = 0.            rugoro = 0.
551         ENDIF         ENDIF
552    
553         ecrit_ins = NINT(ecrit_ins/dtphys)         ecrit_ins = NINT(ecrit_ins / dtphys)
        ecrit_hf = NINT(ecrit_hf/dtphys)  
        ecrit_mth = NINT(ecrit_mth/dtphys)  
        ecrit_tra = NINT(86400.*ecrit_tra/dtphys)  
        ecrit_reg = NINT(ecrit_reg/dtphys)  
554    
555         ! Initialisation des sorties         ! Initialisation des sorties
556    
557         call ini_histins(dtphys, ok_instan, nid_ins, itau_phy)         call ini_histins(dtphys)
558         CALL ymds2ju(annee_ref, 1, day_ref, 0., date0)         CALL ymds2ju(annee_ref, 1, day_ref, 0., date0)
559         ! Positionner date0 pour initialisation de ORCHIDEE         ! Positionner date0 pour initialisation de ORCHIDEE
560         print *, 'physiq date0: ', date0         print *, 'physiq date0: ', date0
561         CALL phyredem0(lmt_pas, itau_phy)         CALL phyredem0
562      ENDIF test_firstcal      ENDIF test_firstcal
563    
564      ! We will modify variables *_seri and we will not touch variables      ! We will modify variables *_seri and we will not touch variables
# Line 635  contains Line 614  contains
614      ! Check temperatures:      ! Check temperatures:
615      CALL hgardfou(t_seri, ftsol)      CALL hgardfou(t_seri, ftsol)
616    
617      ! Incrémenter le compteur de la physique      call increment_itap
     itap = itap + 1  
618      julien = MOD(dayvrai, 360)      julien = MOD(dayvrai, 360)
619      if (julien == 0) julien = 360      if (julien == 0) julien = 360
620    
# Line 667  contains Line 645  contains
645      frugs = MAX(frugs, 0.000015)      frugs = MAX(frugs, 0.000015)
646      zxrugs = sum(frugs * pctsrf, dim = 2)      zxrugs = sum(frugs * pctsrf, dim = 2)
647    
648      ! Calculs nécessaires au calcul de l'albedo dans l'interface avec      ! Calculs n\'ecessaires au calcul de l'albedo dans l'interface avec
649      ! la surface.      ! la surface.
650    
651      CALL orbite(REAL(julien), longi, dist)      CALL orbite(REAL(julien), longi, dist)
# Line 693  contains Line 671  contains
671    
672      ! Couche limite:      ! Couche limite:
673    
674      CALL clmain(dtphys, itap, pctsrf, pctsrf_new, t_seri, q_seri, u_seri, &      CALL clmain(dtphys, pctsrf, t_seri, q_seri, u_seri, v_seri, julien, mu0, &
675           v_seri, julien, mu0, ftsol, cdmmax, cdhmax, ksta, ksta_ter, &           ftsol, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ftsoil, qsol, &
676           ok_kzmin, ftsoil, qsol, paprs, play, fsnow, fqsurf, fevap, falbe, &           paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, &
677           fluxlat, rain_fall, snow_fall, fsolsw, fsollw, fder, rlat, frugs, &           snow_fall, fsolsw, fsollw, fder, rlat, frugs, agesno, rugoro, &
678           firstcal, agesno, rugoro, d_t_vdf, d_q_vdf, d_u_vdf, d_v_vdf, d_ts, &           d_t_vdf, d_q_vdf, d_u_vdf, d_v_vdf, d_ts, fluxt, fluxq, fluxu, &
679           fluxt, fluxq, fluxu, fluxv, cdragh, cdragm, q2, dsens, devap, &           fluxv, cdragh, cdragm, q2, dsens, devap, ycoefh, yu1, yv1, t2m, q2m, &
680           ycoefh, yu1, yv1, t2m, q2m, u10m, v10m, pblh, capCL, oliqCL, cteiCL, &           u10m, v10m, pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, &
681           pblT, therm, trmb1, trmb2, trmb3, plcl, fqcalving, ffonte, &           trmb3, plcl, fqcalving, ffonte, run_off_lic_0)
          run_off_lic_0)  
682    
683      ! Incr\'ementation des flux      ! Incr\'ementation des flux
684    
# Line 745  contains Line 722  contains
722      ! Update surface temperature:      ! Update surface temperature:
723    
724      DO i = 1, klon      DO i = 1, klon
        zxtsol(i) = 0.  
725         zxfluxlat(i) = 0.         zxfluxlat(i) = 0.
726    
727         zt2m(i) = 0.         zt2m(i) = 0.
# Line 765  contains Line 741  contains
741         s_trmb1(i) = 0.         s_trmb1(i) = 0.
742         s_trmb2(i) = 0.         s_trmb2(i) = 0.
743         s_trmb3(i) = 0.         s_trmb3(i) = 0.
   
        IF (abs(pctsrf(i, is_ter) + pctsrf(i, is_lic) + pctsrf(i, is_oce) &  
             + pctsrf(i, is_sic) - 1.) > EPSFRA) print *, &  
             'physiq : probl\`eme sous surface au point ', i, &  
             pctsrf(i, 1 : nbsrf)  
744      ENDDO      ENDDO
745    
746        call assert(abs(sum(pctsrf, dim = 2) - 1.) <= EPSFRA, 'physiq: pctsrf')
747    
748        ftsol = ftsol + d_ts
749        zxtsol = sum(ftsol * pctsrf, dim = 2)
750      DO nsrf = 1, nbsrf      DO nsrf = 1, nbsrf
751         DO i = 1, klon         DO i = 1, klon
752            ftsol(i, nsrf) = ftsol(i, nsrf) + d_ts(i, nsrf)            zxfluxlat(i) = zxfluxlat(i) + fluxlat(i, nsrf) * pctsrf(i, nsrf)
753            zxtsol(i) = zxtsol(i) + ftsol(i, nsrf)*pctsrf(i, nsrf)  
754            zxfluxlat(i) = zxfluxlat(i) + fluxlat(i, nsrf)*pctsrf(i, nsrf)            zt2m(i) = zt2m(i) + t2m(i, nsrf) * pctsrf(i, nsrf)
755              zq2m(i) = zq2m(i) + q2m(i, nsrf) * pctsrf(i, nsrf)
756            zt2m(i) = zt2m(i) + t2m(i, nsrf)*pctsrf(i, nsrf)            zu10m(i) = zu10m(i) + u10m(i, nsrf) * pctsrf(i, nsrf)
757            zq2m(i) = zq2m(i) + q2m(i, nsrf)*pctsrf(i, nsrf)            zv10m(i) = zv10m(i) + v10m(i, nsrf) * pctsrf(i, nsrf)
758            zu10m(i) = zu10m(i) + u10m(i, nsrf)*pctsrf(i, nsrf)            zxffonte(i) = zxffonte(i) + ffonte(i, nsrf) * pctsrf(i, nsrf)
           zv10m(i) = zv10m(i) + v10m(i, nsrf)*pctsrf(i, nsrf)  
           zxffonte(i) = zxffonte(i) + ffonte(i, nsrf)*pctsrf(i, nsrf)  
759            zxfqcalving(i) = zxfqcalving(i) + &            zxfqcalving(i) = zxfqcalving(i) + &
760                 fqcalving(i, nsrf)*pctsrf(i, nsrf)                 fqcalving(i, nsrf) * pctsrf(i, nsrf)
761            s_pblh(i) = s_pblh(i) + pblh(i, nsrf)*pctsrf(i, nsrf)            s_pblh(i) = s_pblh(i) + pblh(i, nsrf) * pctsrf(i, nsrf)
762            s_lcl(i) = s_lcl(i) + plcl(i, nsrf)*pctsrf(i, nsrf)            s_lcl(i) = s_lcl(i) + plcl(i, nsrf) * pctsrf(i, nsrf)
763            s_capCL(i) = s_capCL(i) + capCL(i, nsrf) *pctsrf(i, nsrf)            s_capCL(i) = s_capCL(i) + capCL(i, nsrf) * pctsrf(i, nsrf)
764            s_oliqCL(i) = s_oliqCL(i) + oliqCL(i, nsrf) *pctsrf(i, nsrf)            s_oliqCL(i) = s_oliqCL(i) + oliqCL(i, nsrf) * pctsrf(i, nsrf)
765            s_cteiCL(i) = s_cteiCL(i) + cteiCL(i, nsrf) *pctsrf(i, nsrf)            s_cteiCL(i) = s_cteiCL(i) + cteiCL(i, nsrf) * pctsrf(i, nsrf)
766            s_pblT(i) = s_pblT(i) + pblT(i, nsrf) *pctsrf(i, nsrf)            s_pblT(i) = s_pblT(i) + pblT(i, nsrf) * pctsrf(i, nsrf)
767            s_therm(i) = s_therm(i) + therm(i, nsrf) *pctsrf(i, nsrf)            s_therm(i) = s_therm(i) + therm(i, nsrf) * pctsrf(i, nsrf)
768            s_trmb1(i) = s_trmb1(i) + trmb1(i, nsrf) *pctsrf(i, nsrf)            s_trmb1(i) = s_trmb1(i) + trmb1(i, nsrf) * pctsrf(i, nsrf)
769            s_trmb2(i) = s_trmb2(i) + trmb2(i, nsrf) *pctsrf(i, nsrf)            s_trmb2(i) = s_trmb2(i) + trmb2(i, nsrf) * pctsrf(i, nsrf)
770            s_trmb3(i) = s_trmb3(i) + trmb3(i, nsrf) *pctsrf(i, nsrf)            s_trmb3(i) = s_trmb3(i) + trmb3(i, nsrf) * pctsrf(i, nsrf)
771         ENDDO         ENDDO
772      ENDDO      ENDDO
773    
# Line 836  contains Line 810  contains
810         da = 0.         da = 0.
811         mp = 0.         mp = 0.
812         phi = 0.         phi = 0.
813         CALL concvl(dtphys, paprs, play, t_seri, q_seri, u_seri, v_seri, sig1, &         CALL concvl(paprs, play, t_seri, q_seri, u_seri, v_seri, sig1, w01, &
814              w01, d_t_con, d_q_con, d_u_con, d_v_con, rain_con, ibas_con, &              d_t_con, d_q_con, d_u_con, d_v_con, rain_con, ibas_con, itop_con, &
815              itop_con, upwd, dnwd, dnwd0, Ma, cape, iflagctrl, qcondc, pmflxr, &              upwd, dnwd, dnwd0, Ma, cape, iflagctrl, qcondc, pmflxr, da, phi, mp)
             da, phi, mp)  
816         snow_con = 0.         snow_con = 0.
817         clwcon0 = qcondc         clwcon0 = qcondc
818         mfu = upwd + dnwd         mfu = upwd + dnwd
# Line 900  contains Line 873  contains
873         zx_t = 0.         zx_t = 0.
874         za = 0.         za = 0.
875         DO i = 1, klon         DO i = 1, klon
876            za = za + airephy(i)/REAL(klon)            za = za + airephy(i) / REAL(klon)
877            zx_t = zx_t + (rain_con(i)+ &            zx_t = zx_t + (rain_con(i)+ &
878                 snow_con(i))*airephy(i)/REAL(klon)                 snow_con(i)) * airephy(i) / REAL(klon)
879         ENDDO         ENDDO
880         zx_t = zx_t/za*dtphys         zx_t = zx_t / za * dtphys
881         print *, "Precip = ", zx_t         print *, "Precip = ", zx_t
882      ENDIF      ENDIF
883    
# Line 935  contains Line 908  contains
908         t_seri = t_seri + d_t_ajs         t_seri = t_seri + d_t_ajs
909         q_seri = q_seri + d_q_ajs         q_seri = q_seri + d_q_ajs
910      else      else
        ! Thermiques  
911         call calltherm(dtphys, play, paprs, pphi, u_seri, v_seri, t_seri, &         call calltherm(dtphys, play, paprs, pphi, u_seri, v_seri, t_seri, &
912              q_seri, d_u_ajs, d_v_ajs, d_t_ajs, d_q_ajs, fm_therm, entr_therm)              q_seri, d_u_ajs, d_v_ajs, d_t_ajs, d_q_ajs, fm_therm, entr_therm)
913      endif      endif
# Line 1006  contains Line 978  contains
978         zx_t = 0.         zx_t = 0.
979         za = 0.         za = 0.
980         DO i = 1, klon         DO i = 1, klon
981            za = za + airephy(i)/REAL(klon)            za = za + airephy(i) / REAL(klon)
982            zx_t = zx_t + (rain_lsc(i) &            zx_t = zx_t + (rain_lsc(i) &
983                 + snow_lsc(i))*airephy(i)/REAL(klon)                 + snow_lsc(i)) * airephy(i) / REAL(klon)
984         ENDDO         ENDDO
985         zx_t = zx_t/za*dtphys         zx_t = zx_t / za * dtphys
986         print *, "Precip = ", zx_t         print *, "Precip = ", zx_t
987      ENDIF      ENDIF
988    
# Line 1036  contains Line 1008  contains
1008            do k = 1, llm            do k = 1, llm
1009               do i = 1, klon               do i = 1, klon
1010                  if (d_q_con(i, k) < 0.) then                  if (d_q_con(i, k) < 0.) then
1011                     rain_tiedtke(i) = rain_tiedtke(i) - d_q_con(i, k)/dtphys &                     rain_tiedtke(i) = rain_tiedtke(i) - d_q_con(i, k) / dtphys &
1012                          *zmasse(i, k)                          * zmasse(i, k)
1013                  endif                  endif
1014               enddo               enddo
1015            enddo            enddo
# Line 1072  contains Line 1044  contains
1044    
1045         ! On prend la somme des fractions nuageuses et des contenus en eau         ! On prend la somme des fractions nuageuses et des contenus en eau
1046         cldfra = min(max(cldfra, rnebcon), 1.)         cldfra = min(max(cldfra, rnebcon), 1.)
1047         cldliq = cldliq + rnebcon*clwcon         cldliq = cldliq + rnebcon * clwcon
1048      ENDIF      ENDIF
1049    
1050      ! 2. Nuages stratiformes      ! 2. Nuages stratiformes
# Line 1104  contains Line 1076  contains
1076         DO i = 1, klon         DO i = 1, klon
1077            zx_t = t_seri(i, k)            zx_t = t_seri(i, k)
1078            IF (thermcep) THEN            IF (thermcep) THEN
1079               zx_qs = r2es * FOEEW(zx_t, rtt >= zx_t)/play(i, k)               zx_qs = r2es * FOEEW(zx_t, rtt >= zx_t) / play(i, k)
1080               zx_qs = MIN(0.5, zx_qs)               zx_qs = MIN(0.5, zx_qs)
1081               zcor = 1./(1. - retv*zx_qs)               zcor = 1. / (1. - retv * zx_qs)
1082               zx_qs = zx_qs*zcor               zx_qs = zx_qs * zcor
1083            ELSE            ELSE
1084               IF (zx_t < t_coup) THEN               IF (zx_t < t_coup) THEN
1085                  zx_qs = qsats(zx_t)/play(i, k)                  zx_qs = qsats(zx_t) / play(i, k)
1086               ELSE               ELSE
1087                  zx_qs = qsatl(zx_t)/play(i, k)                  zx_qs = qsatl(zx_t) / play(i, k)
1088               ENDIF               ENDIF
1089            ENDIF            ENDIF
1090            zx_rh(i, k) = q_seri(i, k)/zx_qs            zx_rh(i, k) = q_seri(i, k) / zx_qs
1091            zqsat(i, k) = zx_qs            zqsat(i, k) = zx_qs
1092         ENDDO         ENDDO
1093      ENDDO      ENDDO
1094    
1095      ! Introduce the aerosol direct and first indirect radiative forcings:      ! Introduce the aerosol direct and first indirect radiative forcings:
1096      IF (ok_ade .OR. ok_aie) THEN      tau_ae = 0.
1097         ! Get sulfate aerosol distribution :      piz_ae = 0.
1098         CALL readsulfate(dayvrai, time, firstcal, sulfate)      cg_ae = 0.
        CALL readsulfate_preind(dayvrai, time, firstcal, sulfate_pi)  
   
        CALL aeropt(play, paprs, t_seri, sulfate, rhcl, tau_ae, piz_ae, cg_ae, &  
             aerindex)  
     ELSE  
        tau_ae = 0.  
        piz_ae = 0.  
        cg_ae = 0.  
     ENDIF  
1099    
1100      ! Param\`etres optiques des nuages et quelques param\`etres pour      ! Param\`etres optiques des nuages et quelques param\`etres pour
1101      ! diagnostics :      ! diagnostics :
# Line 1164  contains Line 1127  contains
1127    
1128      DO k = 1, llm      DO k = 1, llm
1129         DO i = 1, klon         DO i = 1, klon
1130            t_seri(i, k) = t_seri(i, k) + (heat(i, k) - cool(i, k)) * dtphys/86400.            t_seri(i, k) = t_seri(i, k) + (heat(i, k) - cool(i, k)) * dtphys &
1131                   / 86400.
1132         ENDDO         ENDDO
1133      ENDDO      ENDDO
1134    
# Line 1183  contains Line 1147  contains
1147      ENDDO      ENDDO
1148      DO nsrf = 1, nbsrf      DO nsrf = 1, nbsrf
1149         DO i = 1, klon         DO i = 1, klon
1150            zxqsurf(i) = zxqsurf(i) + fqsurf(i, nsrf)*pctsrf(i, nsrf)            zxqsurf(i) = zxqsurf(i) + fqsurf(i, nsrf) * pctsrf(i, nsrf)
1151            zxsnow(i) = zxsnow(i) + fsnow(i, nsrf)*pctsrf(i, nsrf)            zxsnow(i) = zxsnow(i) + fsnow(i, nsrf) * pctsrf(i, nsrf)
1152         ENDDO         ENDDO
1153      ENDDO      ENDDO
1154    
# Line 1269  contains Line 1233  contains
1233           d_qt, d_ec)           d_qt, d_ec)
1234    
1235      ! Calcul des tendances traceurs      ! Calcul des tendances traceurs
1236      call phytrac(itap, lmt_pas, julien, time, firstcal, lafin, dtphys, t, &      call phytrac(julien, time, firstcal, lafin, dtphys, t, paprs, play, mfu, &
1237           paprs, play, mfu, mfd, pde_u, pen_d, ycoefh, fm_therm, entr_therm, &           mfd, pde_u, pen_d, ycoefh, fm_therm, entr_therm, yu1, yv1, ftsol, &
1238           yu1, yv1, ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, &           pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, tr_seri, &
1239           dnwd, tr_seri, zmasse, ncid_startphy, nid_ins, itau_phy)           zmasse, ncid_startphy)
1240    
1241      IF (offline) call phystokenc(dtphys, t, mfu, mfd, pen_u, pde_u, pen_d, &      IF (offline) call phystokenc(dtphys, t, mfu, mfd, pen_u, pde_u, pen_d, &
1242           pde_d, fm_therm, entr_therm, ycoefh, yu1, yv1, ftsol, pctsrf, &           pde_d, fm_therm, entr_therm, ycoefh, yu1, yv1, ftsol, pctsrf, &
1243           frac_impa, frac_nucl, pphis, airephy, dtphys, itap)           frac_impa, frac_nucl, pphis, airephy, dtphys)
1244    
1245      ! Calculer le transport de l'eau et de l'energie (diagnostique)      ! Calculer le transport de l'eau et de l'energie (diagnostique)
1246      CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq)      CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq)
# Line 1288  contains Line 1252  contains
1252    
1253      ! Accumuler les variables a stocker dans les fichiers histoire:      ! Accumuler les variables a stocker dans les fichiers histoire:
1254    
1255      ! conversion Ec -> E thermique      ! conversion Ec en énergie thermique
1256      DO k = 1, llm      DO k = 1, llm
1257         DO i = 1, klon         DO i = 1, klon
1258            ZRCPD = RCPD * (1. + RVTMP2 * q_seri(i, k))            ZRCPD = RCPD * (1. + RVTMP2 * q_seri(i, k))
# Line 1318  contains Line 1282  contains
1282      DO i = 1, klon      DO i = 1, klon
1283         prw(i) = 0.         prw(i) = 0.
1284         DO k = 1, llm         DO k = 1, llm
1285            prw(i) = prw(i) + q_seri(i, k)*zmasse(i, k)            prw(i) = prw(i) + q_seri(i, k) * zmasse(i, k)
1286         ENDDO         ENDDO
1287      ENDDO      ENDDO
1288    
# Line 1350  contains Line 1314  contains
1314         ENDDO         ENDDO
1315      ENDDO      ENDDO
1316    
1317      call write_histins      CALL histwrite_phy("phis", pphis)
1318        CALL histwrite_phy("aire", airephy)
1319        CALL histwrite_phy("psol", paprs(:, 1))
1320        CALL histwrite_phy("precip", rain_fall + snow_fall)
1321        CALL histwrite_phy("plul", rain_lsc + snow_lsc)
1322        CALL histwrite_phy("pluc", rain_con + snow_con)
1323        CALL histwrite_phy("tsol", zxtsol)
1324        CALL histwrite_phy("t2m", zt2m)
1325        CALL histwrite_phy("q2m", zq2m)
1326        CALL histwrite_phy("u10m", zu10m)
1327        CALL histwrite_phy("v10m", zv10m)
1328        CALL histwrite_phy("snow", snow_fall)
1329        CALL histwrite_phy("cdrm", cdragm)
1330        CALL histwrite_phy("cdrh", cdragh)
1331        CALL histwrite_phy("topl", toplw)
1332        CALL histwrite_phy("evap", evap)
1333        CALL histwrite_phy("sols", solsw)
1334        CALL histwrite_phy("soll", sollw)
1335        CALL histwrite_phy("solldown", sollwdown)
1336        CALL histwrite_phy("bils", bils)
1337        CALL histwrite_phy("sens", - sens)
1338        CALL histwrite_phy("fder", fder)
1339        CALL histwrite_phy("dtsvdfo", d_ts(:, is_oce))
1340        CALL histwrite_phy("dtsvdft", d_ts(:, is_ter))
1341        CALL histwrite_phy("dtsvdfg", d_ts(:, is_lic))
1342        CALL histwrite_phy("dtsvdfi", d_ts(:, is_sic))
1343    
1344        DO nsrf = 1, nbsrf
1345           CALL histwrite_phy("pourc_"//clnsurf(nsrf), pctsrf(:, nsrf) * 100.)
1346           CALL histwrite_phy("fract_"//clnsurf(nsrf), pctsrf(:, nsrf))
1347           CALL histwrite_phy("sens_"//clnsurf(nsrf), fluxt(:, 1, nsrf))
1348           CALL histwrite_phy("lat_"//clnsurf(nsrf), fluxlat(:, nsrf))
1349           CALL histwrite_phy("tsol_"//clnsurf(nsrf), ftsol(:, nsrf))
1350           CALL histwrite_phy("taux_"//clnsurf(nsrf), fluxu(:, 1, nsrf))
1351           CALL histwrite_phy("tauy_"//clnsurf(nsrf), fluxv(:, 1, nsrf))
1352           CALL histwrite_phy("rugs_"//clnsurf(nsrf), frugs(:, nsrf))
1353           CALL histwrite_phy("albe_"//clnsurf(nsrf), falbe(:, nsrf))
1354        END DO
1355    
1356        CALL histwrite_phy("albs", albsol)
1357        CALL histwrite_phy("rugs", zxrugs)
1358        CALL histwrite_phy("s_pblh", s_pblh)
1359        CALL histwrite_phy("s_pblt", s_pblt)
1360        CALL histwrite_phy("s_lcl", s_lcl)
1361        CALL histwrite_phy("s_capCL", s_capCL)
1362        CALL histwrite_phy("s_oliqCL", s_oliqCL)
1363        CALL histwrite_phy("s_cteiCL", s_cteiCL)
1364        CALL histwrite_phy("s_therm", s_therm)
1365        CALL histwrite_phy("s_trmb1", s_trmb1)
1366        CALL histwrite_phy("s_trmb2", s_trmb2)
1367        CALL histwrite_phy("s_trmb3", s_trmb3)
1368        if (conv_emanuel) CALL histwrite_phy("ptop", ema_pct)
1369        CALL histwrite_phy("temp", t_seri)
1370        CALL histwrite_phy("vitu", u_seri)
1371        CALL histwrite_phy("vitv", v_seri)
1372        CALL histwrite_phy("geop", zphi)
1373        CALL histwrite_phy("pres", play)
1374        CALL histwrite_phy("dtvdf", d_t_vdf)
1375        CALL histwrite_phy("dqvdf", d_q_vdf)
1376        CALL histwrite_phy("rhum", zx_rh)
1377    
1378        if (ok_instan) call histsync(nid_ins)
1379    
1380      IF (lafin) then      IF (lafin) then
1381         call NF95_CLOSE(ncid_startphy)         call NF95_CLOSE(ncid_startphy)
# Line 1363  contains Line 1388  contains
1388    
1389      firstcal = .FALSE.      firstcal = .FALSE.
1390    
   contains  
   
     subroutine write_histins  
   
       ! From phylmd/write_histins.h, version 1.2 2005/05/25 13:10:09  
   
       ! Ecriture des sorties  
   
       use gr_phy_write_m, only: gr_phy_write  
       USE histsync_m, ONLY: histsync  
       USE histwrite_m, ONLY: histwrite  
   
       integer itau_w ! pas de temps d'\'ecriture  
   
       !--------------------------------------------------  
   
       IF (ok_instan) THEN  
          itau_w = itau_phy + itap  
          CALL histwrite(nid_ins, "phis", itau_w, gr_phy_write(pphis))  
          CALL histwrite(nid_ins, "aire", itau_w, gr_phy_write(airephy))  
          CALL histwrite(nid_ins, "psol", itau_w, gr_phy_write(paprs(:, 1)))  
          CALL histwrite(nid_ins, "precip", itau_w, &  
               gr_phy_write(rain_fall + snow_fall))  
          CALL histwrite(nid_ins, "plul", itau_w, &  
               gr_phy_write(rain_lsc + snow_lsc))  
          CALL histwrite(nid_ins, "pluc", itau_w, &  
               gr_phy_write(rain_con + snow_con))  
          CALL histwrite(nid_ins, "tsol", itau_w, gr_phy_write(zxtsol))  
          CALL histwrite(nid_ins, "t2m", itau_w, gr_phy_write(zt2m))  
          CALL histwrite(nid_ins, "q2m", itau_w, gr_phy_write(zq2m))  
          CALL histwrite(nid_ins, "u10m", itau_w, gr_phy_write(zu10m))  
          CALL histwrite(nid_ins, "v10m", itau_w, gr_phy_write(zv10m))  
          CALL histwrite(nid_ins, "snow", itau_w, gr_phy_write(snow_fall))  
          CALL histwrite(nid_ins, "cdrm", itau_w, gr_phy_write(cdragm))  
          CALL histwrite(nid_ins, "cdrh", itau_w, gr_phy_write(cdragh))  
          CALL histwrite(nid_ins, "topl", itau_w, gr_phy_write(toplw))  
          CALL histwrite(nid_ins, "evap", itau_w, gr_phy_write(evap))  
          CALL histwrite(nid_ins, "sols", itau_w, gr_phy_write(solsw))  
          CALL histwrite(nid_ins, "soll", itau_w, gr_phy_write(sollw))  
          CALL histwrite(nid_ins, "solldown", itau_w, gr_phy_write(sollwdown))  
          CALL histwrite(nid_ins, "bils", itau_w, gr_phy_write(bils))  
          CALL histwrite(nid_ins, "sens", itau_w, gr_phy_write(- sens))  
          CALL histwrite(nid_ins, "fder", itau_w, gr_phy_write(fder))  
          CALL histwrite(nid_ins, "dtsvdfo", itau_w, &  
               gr_phy_write(d_ts(:, is_oce)))  
          CALL histwrite(nid_ins, "dtsvdft", itau_w, &  
               gr_phy_write(d_ts(:, is_ter)))  
          CALL histwrite(nid_ins, "dtsvdfg", itau_w, &  
               gr_phy_write(d_ts(:, is_lic)))  
          CALL histwrite(nid_ins, "dtsvdfi", itau_w, &  
               gr_phy_write(d_ts(:, is_sic)))  
   
          DO nsrf = 1, nbsrf  
             CALL histwrite(nid_ins, "pourc_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(pctsrf(:, nsrf)*100.))  
             CALL histwrite(nid_ins, "fract_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(pctsrf(:, nsrf)))  
             CALL histwrite(nid_ins, "sens_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(fluxt(:, 1, nsrf)))  
             CALL histwrite(nid_ins, "lat_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(fluxlat(:, nsrf)))  
             CALL histwrite(nid_ins, "tsol_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(ftsol(:, nsrf)))  
             CALL histwrite(nid_ins, "taux_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(fluxu(:, 1, nsrf)))  
             CALL histwrite(nid_ins, "tauy_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(fluxv(:, 1, nsrf)))  
             CALL histwrite(nid_ins, "rugs_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(frugs(:, nsrf)))  
             CALL histwrite(nid_ins, "albe_"//clnsurf(nsrf), itau_w, &  
                  gr_phy_write(falbe(:, nsrf)))  
          END DO  
   
          CALL histwrite(nid_ins, "albs", itau_w, gr_phy_write(albsol))  
          CALL histwrite(nid_ins, "rugs", itau_w, gr_phy_write(zxrugs))  
          CALL histwrite(nid_ins, "s_pblh", itau_w, gr_phy_write(s_pblh))  
          CALL histwrite(nid_ins, "s_pblt", itau_w, gr_phy_write(s_pblt))  
          CALL histwrite(nid_ins, "s_lcl", itau_w, gr_phy_write(s_lcl))  
          CALL histwrite(nid_ins, "s_capCL", itau_w, gr_phy_write(s_capCL))  
          CALL histwrite(nid_ins, "s_oliqCL", itau_w, gr_phy_write(s_oliqCL))  
          CALL histwrite(nid_ins, "s_cteiCL", itau_w, gr_phy_write(s_cteiCL))  
          CALL histwrite(nid_ins, "s_therm", itau_w, gr_phy_write(s_therm))  
          CALL histwrite(nid_ins, "s_trmb1", itau_w, gr_phy_write(s_trmb1))  
          CALL histwrite(nid_ins, "s_trmb2", itau_w, gr_phy_write(s_trmb2))  
          CALL histwrite(nid_ins, "s_trmb3", itau_w, gr_phy_write(s_trmb3))  
          if (conv_emanuel) CALL histwrite(nid_ins, "ptop", itau_w, &  
               gr_phy_write(ema_pct))  
          CALL histwrite(nid_ins, "temp", itau_w, gr_phy_write(t_seri))  
          CALL histwrite(nid_ins, "vitu", itau_w, gr_phy_write(u_seri))  
          CALL histwrite(nid_ins, "vitv", itau_w, gr_phy_write(v_seri))  
          CALL histwrite(nid_ins, "geop", itau_w, gr_phy_write(zphi))  
          CALL histwrite(nid_ins, "pres", itau_w, gr_phy_write(play))  
          CALL histwrite(nid_ins, "dtvdf", itau_w, gr_phy_write(d_t_vdf))  
          CALL histwrite(nid_ins, "dqvdf", itau_w, gr_phy_write(d_q_vdf))  
          CALL histwrite(nid_ins, "rhum", itau_w, gr_phy_write(zx_rh))  
          call histsync(nid_ins)  
       ENDIF  
   
     end subroutine write_histins  
   
1391    END SUBROUTINE physiq    END SUBROUTINE physiq
1392    
1393  end module physiq_m  end module physiq_m
Legend:
Removed from v.190  
changed lines
  Added in v.202
  ViewVC Help
Powered by ViewVC 1.1.21