source: trunk/LMDZ.VENUS/libf/phyvenus/write_histmth.h @ 1543

Last change on this file since 1543 was 1543, checked in by emillour, 9 years ago

All models: Further adaptations to keep up with changes in LMDZ5 concerning
physics/dynamics separation:

  • dyn3d:
  • adapted gcm.F so that all physics initializations are now done in iniphysiq.
  • dyn3dpar:
  • adapted gcm.F so that all physics initializations are now done in iniphysiq.
  • updated calfis_p.F to follow up with changes.
  • copied over updated "bands.F90" from LMDZ5.
  • dynphy_lonlat:
  • calfis_p.F90, mod_interface_dyn_phys.F90, follow up of changes in phy_common/mod_* routines
  • phy_common:
  • added "geometry_mod.F90" to store information about the grid (replaces phy*/comgeomphy.F90) and give variables friendlier names: rlond => longitude , rlatd => latitude, airephy => cell_area, cuphy => dx , cvphy => dy
  • added "physics_distribution_mod.F90"
  • updated "mod_grid_phy_lmdz.F90", "mod_phys_lmdz_mpi_data.F90", "mod_phys_lmdz_para.F90", "mod_phys_lmdz_mpi_transfert.F90", "mod_grid_phy_lmdz.F90", "mod_phys_lmdz_omp_data.F90", "mod_phys_lmdz_omp_transfert.F90", "write_field_phy.F90" and "ioipsl_getin_p_mod.F90" to LMDZ5 versions.
  • phy[venus/titan/mars/std]:
  • removed "init_phys_lmdz.F90", "comgeomphy.F90"; adapted routines to use geometry_mod (longitude, latitude, cell_area, etc.)

EM
File size: 5.3 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/write_histmth.h,v 1.2 2005/03/09 12:30:16 fairhead Exp $
3!
4      IF (ok_mensuel) THEN
5
6         itau_w = itau_phy + itap
7
8c-------------------------------------------------------
9      IF(lev_histmth.GE.1) THEN
10
11ccccccccccccc 2D fields, basics
12
13      call histwrite_phy(nid_mth,.false.,"phis",itau_w,pphis)
14      call histwrite_phy(nid_mth,.false.,"aire",itau_w,cell_area)
15      call histwrite_phy(nid_mth,.false.,"tsol",itau_w,ftsol)
16      call histwrite_phy(nid_mth,.false.,"psol",itau_w,paprs(:,1))
17c     call histwrite_phy(nid_mth,.false.,"ue",itau_w,ue)
18c VENUS: regardee a l envers!!!!!!!!!!!!!!!
19c     call histwrite_phy(nid_mth,.false.,"ve",itau_w,-1.*ve)
20c     call histwrite_phy(nid_mth,.false.,"cdragh",itau_w,cdragh)
21c     call histwrite_phy(nid_mth,.false.,"cdragm",itau_w,cdragm)
22
23      ENDIF !lev_histmth.GE.1
24
25c-------------------------------------------------------
26      IF(lev_histmth.GE.2) THEN
27
28ccccccccccccc 3D fields, basics
29
30      call histwrite_phy(nid_mth,.false.,"temp",itau_w,t_seri)
31      call histwrite_phy(nid_mth,.false.,"pres",itau_w,pplay)
32      call histwrite_phy(nid_mth,.false.,"geop",itau_w,zphi)
33      call histwrite_phy(nid_mth,.false.,"vitu",itau_w,u_seri)
34c VENUS: regardee a l envers!!!!!!!!!!!!!!!
35      call histwrite_phy(nid_mth,.false.,"vitv",itau_w,-1.*v_seri)
36      call histwrite_phy(nid_mth,.false.,"vitw",itau_w,omega)
37c en (m/s)/s     
38      call histwrite_phy(nid_mth,.false.,"dudyn",itau_w,d_u_dyn)
39c en (m/s)/s     
40      call histwrite_phy(nid_mth,.false.,"duvdf",itau_w,d_u_vdf)
41c     call histwrite_phy(nid_mth,.false.,"mang",itau_w,mang)
42c     call histwrite_phy(nid_mth,.false.,"Kz",itau_w,ycoefh)
43      call histwrite_phy(nid_mth,.false.,"mmean",itau_w,mmean)
44c     call histwrite_phy(nid_mth,.false.,"rho",itau_w,rho)
45
46c plusieurs traceurs  !!!outputs in [vmr]
47       IF (iflag_trac.eq.1) THEN
48         DO iq=1,nqmax
49       call histwrite_phy(nid_mth,.false.,tname(iq),itau_w,qx(:,:,iq)
50     &                   *mmean(:,:)/M_tr(iq))
51         ENDDO
52       ENDIF
53
54       IF (callthermos .and. ok_chem) THEN
55       call histwrite_phy(nid_mth,.false.,"d_qmoldifCO2",itau_w,
56     .                 d_q_moldif(:,:,i_co2))
57       call histwrite_phy(nid_mth,.false.,"d_qmoldifO3p",itau_w,
58     .                  d_q_moldif(:,:,i_o))
59       call histwrite_phy(nid_mth,.false.,"d_qmoldifN2",itau_w,
60     .                  d_q_moldif(:,:,i_n2))
61       ENDIF
62
63      call histwrite_phy(nid_mth,.false.,"tops",itau_w,topsw)
64
65      ENDIF !lev_histmth.GE.2
66
67c-------------------------------------------------------
68      IF(lev_histmth.GE.3) THEN
69
70cccccccccccccccccc  Radiative transfer
71
72c 2D
73
74      call histwrite_phy(nid_mth,.false.,"topl",itau_w,toplw)
75      call histwrite_phy(nid_mth,.false.,"sols",itau_w,solsw)
76      call histwrite_phy(nid_mth,.false.,"soll",itau_w,sollw)
77
78c 3D
79
80      call histwrite_phy(nid_mth,.false.,"SWnet",itau_w,swnet)
81      call histwrite_phy(nid_mth,.false.,"LWnet",itau_w,lwnet)
82c     call histwrite_phy(nid_mth,.false.,"fluxvdf",itau_w,fluxt)
83c     call histwrite_phy(nid_mth,.false.,"fluxdyn",itau_w,flux_dyn)
84c     call histwrite_phy(nid_mth,.false.,"fluxajs",itau_w,flux_ajs)
85c     call histwrite_phy(nid_mth,.false.,"fluxec",itau_w,flux_ec)
86
87      ENDIF !lev_histmth.GE.3
88
89c-------------------------------------------------------
90      IF(lev_histmth.GE.4) THEN
91
92c en K/s     
93      call histwrite_phy(nid_mth,.false.,"dtdyn",itau_w,d_t_dyn)
94c en K/s     
95c     call histwrite_phy(nid_mth,.false.,"dtphy",itau_w,d_t)
96c en K/s     
97      call histwrite_phy(nid_mth,.false.,"dtvdf",itau_w,d_t_vdf)
98c en K/s     
99      call histwrite_phy(nid_mth,.false.,"dtajs",itau_w,d_t_ajs)
100c en K/s     
101      call histwrite_phy(nid_mth,.false.,"dtswr",itau_w,dtsw)
102      call histwrite_phy(nid_mth,.false.,"dtswrNLTE",itau_w,d_t_nirco2)
103      call histwrite_phy(nid_mth,.false.,"dtswrDCrisp",itau_w,heat)
104c en K/s     
105      call histwrite_phy(nid_mth,.false.,"dtlwr",itau_w,dtlw)
106c en K/s     
107      call histwrite_phy(nid_mth,.false.,"dtlwrNLTE",itau_w,d_t_nlte)
108      call histwrite_phy(nid_mth,.false.,"dtlwrLTE",itau_w,-1.*cool)
109      call histwrite_phy(nid_mth,.false.,"dteuv",itau_w,d_t_euv)
110      call histwrite_phy(nid_mth,.false.,"dtcond",itau_w,d_t_conduc)
111      call histwrite_phy(nid_mth,.false.,"dumolvis",itau_w,d_u_molvis)
112      call histwrite_phy(nid_mth,.false.,"dvmolvis",itau_w,d_v_molvis)
113
114c     call histwrite_phy(nid_mth,.false.,"dtec",itau_w,d_t_ec)
115c en (m/s)/s     
116      call histwrite_phy(nid_mth,.false.,"duajs",itau_w,d_u_ajs)
117c en (m/s)/s     
118      call histwrite_phy(nid_mth,.false.,"dugwo",itau_w,d_u_oro)
119c en (m/s)/s     
120      call histwrite_phy(nid_mth,.false.,"dugwno",itau_w,d_u_hin)
121c en (m/s)/s     
122
123   
124c     VENUS: regardee a l envers!!!!!!!!!!!!!!!
125c     call histwrite_phy(nid_mth,.false.,"dvvdf",itau_w,-1.*d_v_vdf)
126
127      ENDIF !lev_histmth.GE.4
128
129c-------------------------------------------------------
130      IF(lev_histmth.GE.5) THEN
131
132c     call histwrite_phy(nid_mth,.false.,"taux_",itau_w,fluxu)
133c     call histwrite_phy(nid_mth,.false.,"tauy_",itau_w,fluxv)
134c     call histwrite_phy(nid_mth,.false.,"cdrm",itau_w,cdragm)
135c     call histwrite_phy(nid_mth,.false.,"cdrh",itau_w,cdragh)
136
137      ENDIF !lev_histmth.GE.5
138c-------------------------------------------------------
139
140      if (ok_sync) then
141        call histsync(nid_mth)
142      endif
143
144      ENDIF
145
Note: See TracBrowser for help on using the repository browser.