source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_paramLMDZ_phy.h @ 1533

Last change on this file since 1533 was 1533, checked in by musat, 13 years ago

Ajout fichier 1D mensuel paramLMDZ_phy.nc contenant

  • les parametres orbitaux
  • les taux des GES
  • les moyennes globales ((bils, evap, evap_land, flat, nettop0, nettop, precip, tsol, t2m, prw)

IM

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 3.3 KB
RevLine 
[683]1c
[1533]2c calcul moyennes globales
[683]3c
[1533]4       zx_tmp_fi2d=bils*airephy
5       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gbils)
6       zx_tmp_fi2d=evap*airephy
7       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gevap)
8       zx_tmp_fi2d(:)=fevap(:, is_ter)*airephy(:)
9       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gevapt)
10       zx_tmp_fi2d=zxfluxlat*airephy
11       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,glat)
12       zx_tmp_fi2d=(topsw0-toplw0)*airephy
13       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gnet0)
14       zx_tmp_fi2d=(topsw-toplw)*airephy
15       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gnet)
16       zx_tmp_fi2d=(rain_fall+snow_fall)*airephy
17       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,grain)
18       zx_tmp_fi2d=zxtsol*airephy
19       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gtsol)
20       zx_tmp_fi2d=zt2m*airephy
21       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gt2m)
22       zx_tmp_fi2d=prw*airephy
23       CALL global_mean(zx_tmp_fi2d,airephy,.TRUE.,gprw)
[683]24c
[1533]25c$OMP MASTER
26      if (is_mpi_root) then     
[683]27c
[1533]28      ndex2d = 0
29      itau_w = itau_phy + itap
[683]30c
[1533]31c Variables globales
[683]32c
[1533]33      zx_tmp_0d=R_ecc
[683]34      CALL histwrite(nid_ctesGCM,"R_ecc",itau_w,
[1533]35     .               zx_tmp_0d,np,ndex2d)
[683]36c
[1533]37      zx_tmp_0d=R_peri
[683]38      CALL histwrite(nid_ctesGCM,"R_peri",itau_w,
[1533]39     .               zx_tmp_0d,np,ndex2d)
[683]40c
[1533]41      zx_tmp_0d=R_incl
[683]42      CALL histwrite(nid_ctesGCM,"R_incl",itau_w,
[1533]43     .               zx_tmp_0d,np,ndex2d)
[683]44c
[1533]45      zx_tmp_0d=solaire
[683]46      CALL histwrite(nid_ctesGCM,"solaire",itau_w,
[1533]47     .               zx_tmp_0d,np,ndex2d)
[683]48c
[1533]49      zx_tmp_0d=co2_ppm
[683]50      CALL histwrite(nid_ctesGCM,"co2_ppm",itau_w,
[1533]51     .               zx_tmp_0d,np,ndex2d)
[683]52c
[1533]53      zx_tmp_0d=CH4_ppb
[683]54      CALL histwrite(nid_ctesGCM,"CH4_ppb",itau_w,
[1533]55     .               zx_tmp_0d,np,ndex2d)
[683]56c
[1533]57      zx_tmp_0d=N2O_ppb
[683]58      CALL histwrite(nid_ctesGCM,"N2O_ppb",itau_w,
[1533]59     .               zx_tmp_0d,np,ndex2d)
[683]60c
[1533]61      zx_tmp_0d=CFC11_ppt
[683]62      CALL histwrite(nid_ctesGCM,"CFC11_ppt",itau_w,
[1533]63     .               zx_tmp_0d,np,ndex2d)
[683]64c
[1533]65      zx_tmp_0d=CFC12_ppt
[683]66      CALL histwrite(nid_ctesGCM,"CFC12_ppt",itau_w,
[1533]67     .               zx_tmp_0d,np,ndex2d)
[683]68c
[1533]69c=================================================================
70c moyennes globales
[683]71c
[1533]72      CALL histwrite(nid_ctesGCM,"bils",itau_w,
73     .               gbils,np,ndex2d)
74      CALL histwrite(nid_ctesGCM,"evap",itau_w,
75     .               gevap,np,ndex2d)
76      CALL histwrite(nid_ctesGCM,"evap_land",itau_w,
77     .               gevapt,np,ndex2d)
78      CALL histwrite(nid_ctesGCM,"flat",itau_w,
79     .               glat,np,ndex2d)
80      CALL histwrite(nid_ctesGCM,"nettop0",itau_w,
81     .               gnet0,np,ndex2d)
82      CALL histwrite(nid_ctesGCM,"nettop",itau_w,
83     .               gnet,np,ndex2d)
84      CALL histwrite(nid_ctesGCM,"precip",itau_w,
85     .               grain,np,ndex2d)
86      CALL histwrite(nid_ctesGCM,"tsol",itau_w,
87     .               gtsol,np,ndex2d)
88      CALL histwrite(nid_ctesGCM,"t2m",itau_w,
89     .               gt2m,np,ndex2d)
90      CALL histwrite(nid_ctesGCM,"prw",itau_w,
91     .               gprw,np,ndex2d)
[683]92c=================================================================
93c
94      if (ok_sync) then
95        call histsync(nid_ctesGCM)
96      endif
97c
[1533]98      endif !(is_mpi_root) then     
99c$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.