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

Last change on this file since 2047 was 1572, checked in by emillour, 8 years ago

All GCMs:
Further adaptations to keep up with changes in LMDZ5 concerning
physics/dynamics separation (up to rev r2500 of LMDZ5)

  • arch:
  • remove ifort debug option '-check all', replace it with '-check bounds,format,output_conversion,pointers,uninit' (i.e. get it to stop complaining about copying into temporary arrays)
  • dyn3d_common:
  • comconst_mod.F90 : add ngroup
  • dyn3d:
  • gcm.F90 : minor bug fix (arguments to a call_abort())
  • leapfrog.F90 : recompute geopotential for bilan_dyn outputs
  • conf_gcm.F90 : read "ngroup" from run.def
  • groupe.F , groupeun.F : ngroup no longer a local parameter
  • dyn3d_par:
  • conf_gcm.F90 : read "ngroup" from run.def
  • groupe_p.F , groupeun_p.F : ngroup no longer a local parameter
  • misc:
  • regr1_step_av_m.F90 : removed (not used)
  • phy_common:
  • mod_phys_lmdz_mpi_transfert.F90 , mod_phys_lmdz_mpi_data.F90 : change is_north_pole and is_south_pole to is_north_pole_dyn and is_south_pole_dyn
  • mod_phys_lmdz_omp_data.F90 : introduce is_nort_pole_phy and is_south_pole_phy
  • dynphy_lonlat:
  • mod_interface_dyn_phys.F90 : use is_north_pole_dyn and is_south_pole_dyn
  • calfis_p.F : use is_north_pole_dyn and is_south_pole_dyn
  • phyvenus:
  • physiq_mod , write_hist*.h : use is_north_pole_phy and is_south_pole_phy to correctly compute mesh area at poles to send to hist*nc files.
  • phytitan:
  • physiq_mod , write_hist*.h : use is_north_pole_phy and is_south_pole_phy to correctly compute mesh area at poles to send to hist*nc files.

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