source: LMDZ5/branches/LF-private/libf/phylmd/write_paramLMDZ_phy.h @ 5442

Last change on this file since 5442 was 1877, checked in by Laurent Fairhead, 11 years ago

Ménage sur le code pour éliminer les calculs spécifiques ISCCP hors COSP
Un certain nombre de variables non utilisées dans physiq.F90 ont aussi
été supprimée. Environ 840 lignes supprimées du code physiq.F90

Code cleanup to eliminate specific references to ISCCP outside the COSP
library. Unused variables have been cleaned up from the physiq.F90 routine
as well.

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