source: LMDZ5/trunk/libf/phylmd/write_histday_seri.h @ 4062

Last change on this file since 4062 was 2421, checked in by Ehouarn Millour, 9 years ago

Stop sending pducov from dynamics to physics; it is not used in the physics anymore (and what was done was only valid when runnning in serial).
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.9 KB
RevLine 
[1862]1!
2! $Header$
3!
[776]4      IF (is_sequential) THEN
5     
[644]6      IF (type_run.EQ."AMIP") THEN
[1862]7!
[644]8      ndex2d = 0
[2343]9      itau_w = itau_phy + itap + start_time * day_step_phy
[1862]10!
11! Champs 2D:
12!
[644]13      pi = ACOS(-1.)
14      pir = 4.0*ATAN(1.0) / 180.0
[1862]15!
[644]16      DO i=1, klon
17       zx_tmp_fi2d(i)=(topsw(i)-toplw(i))
18      ENDDO
[1862]19!
[644]20      ok_msk=.FALSE.
21      msk(1:klon)=pctsrf(1:klon,is_ter)
[2351]22      CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area,  &
[1862]23           ok_msk, msk, moyglo)
[644]24      zx_tmp_fi2d(1:klon)=moyglo
[1862]25!
[2344]26      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
[1862]27      CALL histwrite(nid_day_seri,"bilTOA",itau_w, &
[2344]28                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
[1862]29!
[644]30      ok_msk=.FALSE.
[2351]31      CALL moyglo_pondaire(klon, bils, cell_area,  &
[1862]32           ok_msk, msk, moyglo)
[644]33      zx_tmp_fi2d(1:klon)=moyglo
[1862]34!
[2344]35      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
[1862]36      CALL histwrite(nid_day_seri,"bils",itau_w, &
[2344]37                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
[1862]38!
[644]39      DO k=1, klev
40      DO i=1, klon
[1862]41!IM 080904    zx_tmp_fi3d(i,k)=u(i,k)**2+v(i,k)**2
[644]42       zx_tmp_fi3d(i,k)=(u(i,k)**2+v(i,k)**2)/2.
43      ENDDO
44      ENDDO
[1862]45!
46      CALL moyglo_pondaima(klon, klev, zx_tmp_fi3d,  &
[2351]47           cell_area, paprs, moyglo)
[644]48      zx_tmp_fi2d(1:klon)=moyglo
[1862]49!
[2344]50      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
[1862]51      CALL histwrite(nid_day_seri,"ecin",itau_w, &
[2344]52                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d) 
[1862]53!
54!
55!
[2351]56      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,cell_area,zx_tmp_2d)
[644]57      airetot=0.
58      DO i=1, klon
[2351]59       airetot=airetot+cell_area(i)
[644]60      ENDDO
[1862]61!     IF(itap.EQ.1) PRINT*,'airetotphy=',airetot
62!
[644]63      airetot=0.
[2344]64      DO j=1, nbp_lat
65       DO i=1, nbp_lon
[644]66        airetot=airetot+zx_tmp_2d(i,j)
67       ENDDO
68      ENDDO
[1862]69!
70!     IF(itap.EQ.1) PRINT*,'airetotij=',airetot,
71!    $ '4piR2',4.*pi*RA*RA
72!
[644]73      zx_tmp_fi2d(1:klon)=aam/airetot
[2344]74      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
[1862]75      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d, &
[2344]76                     nbp_lon*nbp_lat,ndex2d)
[1862]77!
[644]78      zx_tmp_fi2d(1:klon)=torsfc/airetot
[2344]79      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
[1862]80      CALL histwrite(nid_day_seri,"torsfc",itau_w,zx_tmp_2d, &
[2344]81                     nbp_lon*nbp_lat,ndex2d)
[1862]82!
83!IM 151004 END
84!
85      CALL moyglo_pondmass(klon, klev, t_seri, &
[2351]86           cell_area, paprs, moyglo)
[644]87      zx_tmp_fi2d(1:klon)=moyglo
[1862]88!
[2344]89      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
[1862]90      CALL histwrite(nid_day_seri,"tamv",itau_w, &
[2344]91                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
[1862]92!
[644]93      ok_msk=.FALSE.
[2351]94      CALL moyglo_pondaire(klon, paprs(:,1), cell_area,  &
[1862]95           ok_msk, msk, moyglo)
[644]96      zx_tmp_fi2d(1:klon)=moyglo
[1862]97!
[2344]98      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
[1862]99      CALL histwrite(nid_day_seri,"psol",itau_w, &
[2344]100                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
[1862]101!
[644]102      ok_msk=.FALSE.
[2351]103      CALL moyglo_pondaire(klon, evap, cell_area,  &
[1862]104           ok_msk, msk, moyglo)
[644]105      zx_tmp_fi2d(1:klon)=moyglo
[1862]106!
[2344]107      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
[1862]108      CALL histwrite(nid_day_seri,"evap",itau_w, &
[2344]109                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
[1862]110!
111!     DO i=1, klon
112!      zx_tmp_fi2d(i)=SnowFrac(i,is_ter)
113!     ENDDO
114!
115!     ok_msk=.TRUE.
116!     msk(1:klon)=pctsrf(1:klon,is_ter)
[2351]117!     CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area, 
[1862]118!    .                     ok_msk, msk, moyglo)
119!     zx_tmp_fi2d(1:klon)=moyglo
120!
[2344]121!     CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
[1862]122!     CALL histwrite(nid_day_seri,"SnowFrac",
[2344]123!    .               itau_w,zx_tmp_2d,nbp_lon*nbp_lat,ndex2d) 
[1862]124!
125!     DO i=1, klon
126!IM 080904    zx_tmp_fi2d(i)=zsnow_mass(i)/330.*rowl
127!      zx_tmp_fi2d(i)=zsnow_mass(i)
128!     ENDDO
129!
130!IM 140904   ok_msk=.FALSE.
131!     ok_msk=.TRUE.
132!     msk(1:klon)=pctsrf(1:klon,is_ter)
[2351]133!     CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area, 
[1862]134!    .     ok_msk, msk, moyglo)
135!     zx_tmp_fi2d(1:klon)=moyglo
136!
[2344]137!     CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
[1862]138!     CALL histwrite(nid_day_seri,"snow_depth",itau_w,
[2344]139!    .               zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
[1862]140!
[644]141      DO i=1, klon
142       zx_tmp_fi2d(i)=ftsol(i,is_oce)
143      ENDDO
[1862]144!
[644]145      ok_msk=.TRUE.
146      msk(1:klon)=pctsrf(1:klon,is_oce)
[2351]147      CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area,  &
[1862]148           ok_msk, msk, moyglo)
[644]149      zx_tmp_fi2d(1:klon)=moyglo
[1862]150!
[2344]151      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d, zx_tmp_2d)
[1862]152      CALL histwrite(nid_day_seri,"tsol_"//clnsurf(is_oce), &
[2344]153                     itau_w,zx_tmp_2d,nbp_lon*nbp_lat,ndex2d) 
[1862]154!
155!=================================================================
156!=================================================================
157!=================================================================
158!
[644]159      if (ok_sync) then
160        call histsync(nid_day_seri)
161      endif
[1862]162!
[644]163      ENDIF !fin test sur type_run.EQ."AMIP"
[766]164     
165      ENDIF  ! mono_cpu
Note: See TracBrowser for help on using the repository browser.