/[lmdze]/trunk/libf/phylmd/phystokenc.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/phystokenc.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC revision 31 by guez, Thu Apr 1 14:59:19 2010 UTC
# Line 10  c Line 10  c
10       I                   pcoefh,yu1,yv1,ftsol,pctsrf,       I                   pcoefh,yu1,yv1,ftsol,pctsrf,
11       I                   frac_impa,frac_nucl,       I                   frac_impa,frac_nucl,
12       I                   pphis,paire,dtime,itap)       I                   pphis,paire,dtime,itap)
13        USE ioipsl        USE histwrite_m
14          use histcom
15        use dimens_m        use dimens_m
16        use indicesol        use indicesol
17        use dimphy        use dimphy
# Line 115  c   Couche limite: Line 116  c   Couche limite:
116  c======================================================================  c======================================================================
117    
118        ok_sync = .true.        ok_sync = .true.
         print*,'Dans phystokenc.F'  
       print*,'iadvtr= ',iadvtr  
       print*,'istphy= ',istphy  
       print*,'istdyn= ',istdyn  
119    
120        IF (iadvtr.eq.0) THEN        IF (iadvtr.eq.0) THEN
           
121          CALL initphysto('phystoke',          CALL initphysto('phystoke',
122       . rlon,rlat,dtime, dtime*istphy,dtime*istphy,nqmx,physid)       . rlon,rlat,dtime, dtime*istphy,dtime*istphy,nqmx,physid)
           
         write(*,*) 'apres initphysto ds phystokenc'  
   
           
123        ENDIF        ENDIF
124  c  c
125        ndex2d = 0        ndex2d = 0
126        ndex3d = 0        ndex3d = 0
127        i=itap        i=itap
128        CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)        CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
129        CALL histwrite(physid,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)        CALL histwrite(physid,"phis",i,zx_tmp_2d)
130  c  c
131        i=itap        i=itap
132        CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)        CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
133        CALL histwrite(physid,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)        CALL histwrite(physid,"aire",i,zx_tmp_2d)
134    
135        iadvtr=iadvtr+1        iadvtr=iadvtr+1
136  c  c
# Line 245  c Line 237  c
237    
238  ccccc  ccccc
239           CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d)           CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d)
240           CALL histwrite(physid,"t",itap,zx_tmp_3d,           CALL histwrite(physid,"t",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
241    
242           CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d)           CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d)
243        CALL histwrite(physid,"mfu",itap,zx_tmp_3d,        CALL histwrite(physid,"mfu",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
244          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfd, zx_tmp_3d)          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfd, zx_tmp_3d)
245        CALL histwrite(physid,"mfd",itap,zx_tmp_3d,        CALL histwrite(physid,"mfd",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
246          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_u, zx_tmp_3d)          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_u, zx_tmp_3d)
247        CALL histwrite(physid,"en_u",itap,zx_tmp_3d,        CALL histwrite(physid,"en_u",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
248          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_u, zx_tmp_3d)          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_u, zx_tmp_3d)
249        CALL histwrite(physid,"de_u",itap,zx_tmp_3d,        CALL histwrite(physid,"de_u",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
250          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_d, zx_tmp_3d)          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_d, zx_tmp_3d)
251        CALL histwrite(physid,"en_d",itap,zx_tmp_3d,        CALL histwrite(physid,"en_d",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
252          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_d, zx_tmp_3d)                CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_d, zx_tmp_3d)      
253        CALL histwrite(physid,"de_d",itap,zx_tmp_3d,            CALL histwrite(physid,"de_d",itap,zx_tmp_3d)
      .                                   iim*(jjm+1)*klev,ndex3d)  
254          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, coefh, zx_tmp_3d)                  CALL gr_fi_ecrit(klev,klon,iim,jjm+1, coefh, zx_tmp_3d)        
255        CALL histwrite(physid,"coefh",itap,zx_tmp_3d,            CALL histwrite(physid,"coefh",itap,zx_tmp_3d)    
      .                                   iim*(jjm+1)*klev,ndex3d)        
256    
257  c ajou...  c ajou...
258          do k=1,klev          do k=1,klev
# Line 278  c ajou... Line 262  c ajou...
262          enddo          enddo
263    
264        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, fm_therm1, zx_tmp_3d)        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, fm_therm1, zx_tmp_3d)
265        CALL histwrite(physid,"fm_th",itap,zx_tmp_3d,        CALL histwrite(physid,"fm_th",itap,zx_tmp_3d)
      .                                 iim*(jjm+1)*klev,ndex3d)  
266  c  c
267        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, entr_therm, zx_tmp_3d)        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, entr_therm, zx_tmp_3d)
268        CALL histwrite(physid,"en_th",itap,zx_tmp_3d,        CALL histwrite(physid,"en_th",itap,zx_tmp_3d)
      .                                iim*(jjm+1)*klev,ndex3d)  
269  cccc  cccc
270         CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_impa,zx_tmp_3d)         CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_impa,zx_tmp_3d)
271          CALL histwrite(physid,"frac_impa",itap,zx_tmp_3d,          CALL histwrite(physid,"frac_impa",itap,zx_tmp_3d)
      .  iim*(jjm+1)*klev,ndex3d)  
272    
273          CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_nucl,zx_tmp_3d)          CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_nucl,zx_tmp_3d)
274          CALL histwrite(physid,"frac_nucl",itap,zx_tmp_3d,          CALL histwrite(physid,"frac_nucl",itap,zx_tmp_3d)
      .  iim*(jjm+1)*klev,ndex3d)  
275    
276          CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d)          CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d)
277        CALL histwrite(physid,"pyu1",itap,zx_tmp_2d,iim*(jjm+1),        CALL histwrite(physid,"pyu1",itap,zx_tmp_2d)
      .                                                ndex2d)  
278                    
279          CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d)          CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d)
280        CALL histwrite(physid,"pyv1",itap,zx_tmp_2d,iim*(jjm+1)        CALL histwrite(physid,"pyv1",itap,zx_tmp_2d)
      .                                                ,ndex2d)  
281                    
282          CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol1, zx_tmp_2d)          CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol1, zx_tmp_2d)
283        CALL histwrite(physid,"ftsol1",itap,zx_tmp_2d,        CALL histwrite(physid,"ftsol1",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
284           CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol2, zx_tmp_2d)           CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol2, zx_tmp_2d)
285        CALL histwrite(physid,"ftsol2",itap,zx_tmp_2d,        CALL histwrite(physid,"ftsol2",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
286            CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol3, zx_tmp_2d)            CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol3, zx_tmp_2d)
287        CALL histwrite(physid,"ftsol3",itap,zx_tmp_2d,        CALL histwrite(physid,"ftsol3",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
288           CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d)           CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d)
289        CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d,        CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
290    
291          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf1, zx_tmp_2d)          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf1, zx_tmp_2d)
292        CALL histwrite(physid,"psrf1",itap,zx_tmp_2d,          CALL histwrite(physid,"psrf1",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
293          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf2, zx_tmp_2d)          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf2, zx_tmp_2d)
294        CALL histwrite(physid,"psrf2",itap,zx_tmp_2d,        CALL histwrite(physid,"psrf2",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
295          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf3, zx_tmp_2d)          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf3, zx_tmp_2d)
296        CALL histwrite(physid,"psrf3",itap,zx_tmp_2d,        CALL histwrite(physid,"psrf3",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
297          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf4, zx_tmp_2d)          CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf4, zx_tmp_2d)
298        CALL histwrite(physid,"psrf4",itap,zx_tmp_2d,        CALL histwrite(physid,"psrf4",itap,zx_tmp_2d)
      .                                   iim*(jjm+1),ndex2d)  
299    
300        if (ok_sync) call histsync(physid)        if (ok_sync) call histsync(physid)
301  c     if (ok_sync) call histsync  c     if (ok_sync) call histsync
Legend:
Removed from v.12  
changed lines
  Added in v.31
  ViewVC Help
Powered by ViewVC 1.1.21