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

Last change on this file since 1892 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
RevLine 
[1862]1!
2! calcul moyennes globales
3!
[1538]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)
[1862]24!
25!$OMP MASTER
[1538]26      if (is_mpi_root) then     
[1862]27!
[1538]28      ndex2d = 0
[1877]29      itau_w = itau_phy + itap + int(start_time * day_step / iphysiq)
[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.