source: LMDZ5/trunk/libf/phylmd/write_paramLMDZ_phy.h @ 2401

Last change on this file since 2401 was 2351, checked in by Ehouarn Millour, 9 years ago

More on physics/dynamics separation and cleanup:

  • Set things up so that all physics-related initializations are done via iniphysiq.
  • Created a "geometry_mod.F90" module in phy_common to store information on the loacl grid (i.e. replaces comgeomphy) and moreover give these variables more obvious names (e.g.: rlond => longitude, rlatd => latitude, airephy => cell_area).
  • removed obsolete comgeomphy.h and comgeomphy.F90

EM
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 3.4 KB
RevLine 
[1862]1!
2! calcul moyennes globales
3!
[2351]4       zx_tmp_fi2d=bils*cell_area
5       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gbils)
6       zx_tmp_fi2d=evap*cell_area
7       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gevap)
8       zx_tmp_fi2d(:)=fevap(:, is_ter)*cell_area(:)
9       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gevapt)
10       zx_tmp_fi2d=zxfluxlat*cell_area
11       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,glat)
12       zx_tmp_fi2d=(topsw0-toplw0)*cell_area
13       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gnet0)
14       zx_tmp_fi2d=(topsw-toplw)*cell_area
15       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gnet)
16       zx_tmp_fi2d=(rain_fall+snow_fall)*cell_area
17       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,grain)
18       zx_tmp_fi2d=zxtsol*cell_area
19       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gtsol)
20       zx_tmp_fi2d=zt2m*cell_area
21       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gt2m)
22       zx_tmp_fi2d=prw*cell_area
23       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gprw)
[1862]24!
25!$OMP MASTER
[1538]26      if (is_mpi_root) then     
[1862]27!
[1538]28      ndex2d = 0
[2343]29      itau_w = itau_phy + itap + int(start_time * day_step_phy)
[1862]30!
31! Variables globales
32!
[1538]33      zx_tmp_0d=R_ecc
[1862]34      CALL histwrite(nid_ctesGCM,"R_ecc",itau_w, &
35                     zx_tmp_0d,np,ndex2d)
36!
[1538]37      zx_tmp_0d=R_peri
[1862]38      CALL histwrite(nid_ctesGCM,"R_peri",itau_w, &
39                     zx_tmp_0d,np,ndex2d)
40!
[1538]41      zx_tmp_0d=R_incl
[1862]42      CALL histwrite(nid_ctesGCM,"R_incl",itau_w, &
43                     zx_tmp_0d,np,ndex2d)
44!
[1538]45      zx_tmp_0d=solaire
[1862]46      CALL histwrite(nid_ctesGCM,"solaire",itau_w, &
47                     zx_tmp_0d,np,ndex2d)
48!
[1538]49      zx_tmp_0d=co2_ppm
[1862]50      CALL histwrite(nid_ctesGCM,"co2_ppm",itau_w, &
51                     zx_tmp_0d,np,ndex2d)
52!
[1538]53      zx_tmp_0d=CH4_ppb
[1862]54      CALL histwrite(nid_ctesGCM,"CH4_ppb",itau_w, &
55                     zx_tmp_0d,np,ndex2d)
56!
[1538]57      zx_tmp_0d=N2O_ppb
[1862]58      CALL histwrite(nid_ctesGCM,"N2O_ppb",itau_w, &
59                     zx_tmp_0d,np,ndex2d)
60!
[1538]61      zx_tmp_0d=CFC11_ppt
[1862]62      CALL histwrite(nid_ctesGCM,"CFC11_ppt",itau_w, &
63                     zx_tmp_0d,np,ndex2d)
64!
[1538]65      zx_tmp_0d=CFC12_ppt
[1862]66      CALL histwrite(nid_ctesGCM,"CFC12_ppt",itau_w, &
67                     zx_tmp_0d,np,ndex2d)
68!
69!=================================================================
70! moyennes globales
71!
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)
92!=================================================================
93!
[683]94      if (ok_sync) then
95        call histsync(nid_ctesGCM)
96      endif
[1862]97!
[1538]98      endif !(is_mpi_root) then     
[1862]99!$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.