source: lmdz_wrf/WRFV3/lmdz/write_paramLMDZ_phy.h @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

  • Property svn:executable set to *
File size: 4.0 KB
Line 
1!c
2!c calcul moyennes globales
3!c
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)
24!c
25!$OMP MASTER
26      if (is_mpi_root) then     
27!c
28      ndex2d = 0
29      itau_w = itau_phy + itap + start_time * day_step / iphysiq
30!c
31!c Variables globales
32!c
33      zx_tmp_0d=R_ecc
34      CALL histwrite(nid_ctesGCM,"R_ecc",itau_w,                                     &
35       &               zx_tmp_0d,np,ndex2d)
36!c
37      zx_tmp_0d=R_peri
38      CALL histwrite(nid_ctesGCM,"R_peri",itau_w,                                    &
39       &               zx_tmp_0d,np,ndex2d)
40!c
41      zx_tmp_0d=R_incl
42      CALL histwrite(nid_ctesGCM,"R_incl",itau_w,                                    &
43       &               zx_tmp_0d,np,ndex2d)
44!c
45      zx_tmp_0d=solaire
46      CALL histwrite(nid_ctesGCM,"solaire",itau_w,                                   &
47       &               zx_tmp_0d,np,ndex2d)
48!c
49      zx_tmp_0d=co2_ppm
50      CALL histwrite(nid_ctesGCM,"co2_ppm",itau_w,                                   &
51       &               zx_tmp_0d,np,ndex2d)
52!c
53      zx_tmp_0d=CH4_ppb
54      CALL histwrite(nid_ctesGCM,"CH4_ppb",itau_w,                                   &
55       &               zx_tmp_0d,np,ndex2d)
56!c
57      zx_tmp_0d=N2O_ppb
58      CALL histwrite(nid_ctesGCM,"N2O_ppb",itau_w,                                   &
59       &               zx_tmp_0d,np,ndex2d)
60!c
61      zx_tmp_0d=CFC11_ppt
62      CALL histwrite(nid_ctesGCM,"CFC11_ppt",itau_w,                                 &
63       &               zx_tmp_0d,np,ndex2d)
64!c
65      zx_tmp_0d=CFC12_ppt
66      CALL histwrite(nid_ctesGCM,"CFC12_ppt",itau_w,                                 &
67       &               zx_tmp_0d,np,ndex2d)
68!c
69!c=================================================================
70!c moyennes globales
71!c
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!c=================================================================
93!c
94      if (ok_sync) then
95        call histsync(nid_ctesGCM)
96      endif
97!c
98      endif !(is_mpi_root) then     
99!$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.