source: trunk/LMDZ.VENUS/libf/phyvenus/write_histins.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: 6.2 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/write_histins.h,v 1.1.1.1 2004/05/19 12:53:09 lmdzadmin Exp $
3!
4      IF (ok_instan) THEN
5
6         itau_w = itau_phy + itap
7
8c-------------------------------------------------------
9      IF(lev_histins.GE.1) THEN
10
11ccccccccccccc 2D fields, basics
12
13      call histwrite_phy(nid_ins,.false.,"phis",itau_w,pphis)
14      call histwrite_phy(nid_ins,.false.,"aire",itau_w,cell_area)
15      call histwrite_phy(nid_ins,.false.,"tsol",itau_w,ftsol)
16      call histwrite_phy(nid_ins,.false.,"psol",itau_w,paprs(:,1))
17c     call histwrite_phy(nid_ins,.false.,"ue",itau_w,ue)
18c VENUS: regardee a l' envers!!!!!!!!!!!!!!!
19c     call histwrite_phy(nid_ins,.false.,"ve",itau_w,-1.*ve)
20c     call histwrite_phy(nid_ins,.false.,"cdragh",itau_w,cdragh)
21c     call histwrite_phy(nid_ins,.false.,"cdragm",itau_w,cdragm)
22
23      ENDIF !lev_histins.GE.1
24
25c-------------------------------------------------------
26      IF(lev_histins.GE.2) THEN
27
28ccccccccccccc 3D fields, basics
29
30      call histwrite_phy(nid_ins,.false.,"temp",itau_w,t_seri)
31      call histwrite_phy(nid_ins,.false.,"pres",itau_w,pplay)
32      call histwrite_phy(nid_ins,.false.,"geop",itau_w,zphi)
33      call histwrite_phy(nid_ins,.false.,"vitu",itau_w,u_seri)
34c VENUS: regardee a l' envers !!!!!!!!!!!!!!!
35      call histwrite_phy(nid_ins,.false.,"vitv",itau_w,-1.*v_seri)
36      call histwrite_phy(nid_ins,.false.,"vitw",itau_w,omega)
37c en (m/s)/s     
38      call histwrite_phy(nid_ins,.false.,"dudyn",itau_w,d_u_dyn)
39c en (m/s)/s     
40      call histwrite_phy(nid_ins,.false.,"duvdf",itau_w,d_u_vdf)
41c     call histwrite_phy(nid_ins,.false.,"mang",itau_w,mang)
42c     call histwrite_phy(nid_ins,.false.,"Kz",itau_w,ycoefh)
43      call histwrite_phy(nid_ins,.false.,"mmean",itau_w,mmean)
44      call histwrite_phy(nid_ins,.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_ins,.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_ins,.false.,"d_qmoldifCO2",itau_w,
56     .                 d_q_moldif(:,:,i_co2))
57       call histwrite_phy(nid_ins,.false.,"d_qmoldifO3p",itau_w,
58     .                  d_q_moldif(:,:,i_o))
59       call histwrite_phy(nid_ins,.false.,"d_qmoldifN2",itau_w,
60     .                  d_q_moldif(:,:,i_n2))
61       ENDIF
62     
63       call histwrite_phy(nid_ins,.false.,"tops",itau_w,topsw)
64     
65       IF (ok_cloud) THEN
66       
67       IF (nb_mode.GE.1) THEN
68      call histwrite_phy(nid_ins,.false.,"NBRTOTm1",itau_w,
69     & NBRTOT(:,:,1))
70     
71c      call histwrite_phy(nid_ins,.false.,"R_MEDIANm1",itau_w,
72c     & R_MEDIAN(:,:,1))
73     
74c      call histwrite_phy(nid_ins,.false.,"STDDEVm1",itau_w,
75c     & STDDEV(:,:,1))
76       
77       IF (nb_mode.GE.2) THEN
78      call histwrite_phy(nid_ins,.false.,"NBRTOTm2",itau_w,
79     & NBRTOT(:,:,2))
80     
81c      call histwrite_phy(nid_ins,.false.,"R_MEDIANm2",itau_w,
82c     & R_MEDIAN(:,:,2)) 
83     
84c      call histwrite_phy(nid_ins,.false.,"STDDEVm2",itau_w,
85c     & STDDEV(:,:,2))
86         
87       IF (nb_mode.GE.3) THEN
88      call histwrite_phy(nid_ins,.false.,"NBRTOTm3",itau_w,
89     & NBRTOT(:,:,3))
90     
91c      call histwrite_phy(nid_ins,.false.,"R_MEDIANm3",itau_w,
92c     & R_MEDIAN(:,:,3))
93     
94c      call histwrite_phy(nid_ins,.false.,"STDDEVm3",itau_w,
95c     & STDDEV(:,:,3))
96
97       ENDIF
98       ENDIF
99       ENDIF
100             
101      call histwrite_phy(nid_ins,.false.,"WH2SO4",itau_w,WH2SO4)
102     
103      call histwrite_phy(nid_ins,.false.,"rho_droplet",itau_w,
104     & rho_droplet)
105                ENDIF
106       IF (ok_sedim) THEN
107     
108      call histwrite_phy(nid_ins,.false.,"d_tr_sed_H2SO4",itau_w,
109     & d_tr_sed(:,:,1))     
110      call histwrite_phy(nid_ins,.false.,"d_tr_sed_H2O",itau_w,
111     & d_tr_sed(:,:,2))
112     
113      call histwrite_phy(nid_ins,.false.,"F_sedim",itau_w,Fsedim)
114                ENDIF             
115      ENDIF !lev_histins.GE.2
116
117c-------------------------------------------------------
118      IF(lev_histins.GE.3) THEN
119
120cccccccccccccccccc  Radiative transfer
121
122c 2D
123
124      call histwrite_phy(nid_ins,.false.,"topl",itau_w,toplw)
125      call histwrite_phy(nid_ins,.false.,"sols",itau_w,solsw)
126      call histwrite_phy(nid_ins,.false.,"soll",itau_w,sollw)
127
128c 3D
129
130      call histwrite_phy(nid_ins,.false.,"SWnet",itau_w,swnet)
131      call histwrite_phy(nid_ins,.false.,"LWnet",itau_w,lwnet)
132c     call histwrite_phy(nid_ins,.false.,"fluxvdf",itau_w,fluxt)
133c     call histwrite_phy(nid_ins,.false.,"fluxdyn",itau_w,flux_dyn)
134c     call histwrite_phy(nid_ins,.false.,"fluxajs",itau_w,flux_ajs)
135c     call histwrite_phy(nid_ins,.false.,"fluxec",itau_w,flux_ec)
136
137      ENDIF !lev_histins.GE.3
138
139c-------------------------------------------------------
140      IF(lev_histins.GE.4) THEN
141
142c en K/s     
143      call histwrite_phy(nid_ins,.false.,"dtdyn",itau_w,d_t_dyn)
144c en K/s     
145c     call histwrite_phy(nid_ins,.false.,"dtphy",itau_w,d_t)
146c en K/s     
147      call histwrite_phy(nid_ins,.false.,"dtvdf",itau_w,d_t_vdf)
148c en K/s     
149      call histwrite_phy(nid_ins,.false.,"dtajs",itau_w,d_t_ajs)
150c K/day ==> K/s
151      call histwrite_phy(nid_ins,.false.,"dtswr",itau_w,heat/RDAY)
152c K/day ==> K/s     
153      call histwrite_phy(nid_ins,.false.,"dtlwr",itau_w,-1.*cool/RDAY)
154c en K/s     
155c     call histwrite_phy(nid_ins,.false.,"dtec",itau_w,d_t_ec)
156c en (m/s)/s     
157      call histwrite_phy(nid_ins,.false.,"duajs",itau_w,d_u_ajs)
158c en (m/s)/s     
159      call histwrite_phy(nid_ins,.false.,"dugwo",itau_w,d_u_oro)
160c en (m/s)/s     
161      call histwrite_phy(nid_ins,.false.,"dugwno",itau_w,d_u_hin)
162c en (m/s)/s     
163c VENUS: regardee a l envers!!!!!!!!!!!!!!!
164c     call histwrite_phy(nid_ins,.false.,"dvvdf",itau_w,-1.*d_v_vdf)
165
166      ENDIF !lev_histins.GE.4
167
168c-------------------------------------------------------
169      IF(lev_histins.GE.5) THEN
170
171c     call histwrite_phy(nid_ins,.false.,"taux_",itau_w,fluxu)
172c     call histwrite_phy(nid_ins,.false.,"tauy_",itau_w,fluxv)
173c     call histwrite_phy(nid_ins,.false.,"cdrm",itau_w,cdragm)
174c     call histwrite_phy(nid_ins,.false.,"cdrh",itau_w,cdragh)
175
176      ENDIF !lev_histins.GE.5
177c-------------------------------------------------------
178
179      if (ok_sync) then
180        call histsync(nid_ins)
181      endif
182
183      ENDIF
Note: See TracBrowser for help on using the repository browser.