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

Last change on this file since 1543 was 1543, checked in by emillour, 10 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: 4.2 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/write_histday.h,v 1.2 2004/06/01 09:27:10 lmdzadmin Exp $
3!
4      IF (ok_journe) THEN
5
6         itau_w = itau_phy + itap
7
8c-------------------------------------------------------
9      IF(lev_histday.GE.1) THEN
10
11ccccccccccccc 2D fields, basics
12
13      call histwrite_phy(nid_day,.false.,"phis",itau_w,pphis)
14      call histwrite_phy(nid_day,.false.,"aire",itau_w,cell_area)
15      call histwrite_phy(nid_day,.false.,"tsol",itau_w,ftsol)
16      call histwrite_phy(nid_day,.false.,"psol",itau_w,paprs(:,1))
17c     call histwrite_phy(nid_day,.false.,"ue",itau_w,ue)
18c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
19c     call histwrite_phy(nid_day,.false.,"ve",itau_w,-1.*ve)
20c     call histwrite_phy(nid_day,.false.,"cdragh",itau_w,cdragh)
21c     call histwrite_phy(nid_day,.false.,"cdragm",itau_w,cdragm)
22
23      ENDIF !lev_histday.GE.1
24
25c-------------------------------------------------------
26      IF(lev_histday.GE.2) THEN
27
28ccccccccccccc 3D fields, basics
29
30      call histwrite_phy(nid_day,.false.,"temp",itau_w,t_seri)
31      call histwrite_phy(nid_day,.false.,"pres",itau_w,pplay)
32      call histwrite_phy(nid_day,.false.,"geop",itau_w,zphi)
33      call histwrite_phy(nid_day,.false.,"vitu",itau_w,u_seri)
34c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
35      call histwrite_phy(nid_day,.false.,"vitv",itau_w,-1.*v_seri)
36      call histwrite_phy(nid_day,.false.,"vitw",itau_w,omega)
37c en (m/s)/s     
38      call histwrite_phy(nid_day,.false.,"dudyn",itau_w,d_u_dyn)
39c en (m/s)/s     
40      call histwrite_phy(nid_day,.false.,"duvdf",itau_w,d_u_vdf)
41c     call histwrite_phy(nid_day,.false.,"mang",itau_w,mang)
42c     call histwrite_phy(nid_day,.false.,"Kz",itau_w,ycoefh)
43
44c plusieurs traceurs
45       IF (iflag_trac.eq.1) THEN
46         DO iq=1,nqmax
47       call histwrite_phy(nid_day,.false.,tname(iq),itau_w,qx(:,:,iq))
48         ENDDO
49       ENDIF
50
51      call histwrite_phy(nid_day,.false.,"tops",itau_w,topsw)
52
53      ENDIF !lev_histday.GE.2
54
55c-------------------------------------------------------
56      IF(lev_histday.GE.3) THEN
57
58cccccccccccccccccc  Radiative transfer
59
60c 2D
61
62      call histwrite_phy(nid_day,.false.,"topl",itau_w,toplw)
63      call histwrite_phy(nid_day,.false.,"sols",itau_w,solsw)
64      call histwrite_phy(nid_day,.false.,"soll",itau_w,sollw)
65
66c 3D
67
68      call histwrite_phy(nid_day,.false.,"SWnet",itau_w,swnet)
69      call histwrite_phy(nid_day,.false.,"LWnet",itau_w,lwnet)
70c     call histwrite_phy(nid_day,.false.,"fluxvdf",itau_w,fluxt)
71c     call histwrite_phy(nid_day,.false.,"fluxdyn",itau_w,flux_dyn)
72c     call histwrite_phy(nid_day,.false.,"fluxajs",itau_w,flux_ajs)
73c     call histwrite_phy(nid_day,.false.,"fluxec",itau_w,flux_ec)
74
75      ENDIF !lev_histday.GE.3
76
77c-------------------------------------------------------
78      IF(lev_histday.GE.4) THEN
79
80c en K/s     
81      call histwrite_phy(nid_day,.false.,"dtdyn",itau_w,d_t_dyn)
82c en K/s     
83c     call histwrite_phy(nid_day,.false.,"dtphy",itau_w,d_t)
84c en K/s     
85      call histwrite_phy(nid_day,.false.,"dtvdf",itau_w,d_t_vdf)
86c en K/s     
87      call histwrite_phy(nid_day,.false.,"dtajs",itau_w,d_t_ajs)
88c en K/s     
89      call histwrite_phy(nid_day,.false.,"dtswr",itau_w,dtsw)
90c en K/s     
91      call histwrite_phy(nid_day,.false.,"dtlwr",itau_w,dtlw)
92c en K/s     
93c     call histwrite_phy(nid_day,.false.,"dtec",itau_w,d_t_ec)
94c en (m/s)/s     
95      call histwrite_phy(nid_day,.false.,"duajs",itau_w,d_u_ajs)
96c en (m/s)/s     
97      call histwrite_phy(nid_day,.false.,"dugwo",itau_w,d_u_oro)
98c en (m/s)/s     
99      call histwrite_phy(nid_day,.false.,"dugwno",itau_w,d_u_hin)
100c en (m/s)/s     
101c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
102c     call histwrite_phy(nid_day,.false.,"dvvdf",itau_w,-1.*d_v_vdf)
103
104      ENDIF !lev_histday.GE.4
105
106c-------------------------------------------------------
107      IF(lev_histday.GE.5) THEN
108
109c     call histwrite_phy(nid_day,.false.,"taux_",itau_w,fluxu)
110c     call histwrite_phy(nid_day,.false.,"tauy_",itau_w,fluxv)
111c     call histwrite_phy(nid_day,.false.,"cdrm",itau_w,cdragm)
112c     call histwrite_phy(nid_day,.false.,"cdrh",itau_w,cdragh)
113
114      ENDIF !lev_histday.GE.5
115c-------------------------------------------------------
116
117      if (ok_sync) then
118        call histsync(nid_day)
119      endif
120
121      ENDIF
Note: See TracBrowser for help on using the repository browser.