copyright (C) 2001 MSC-RPN COMM %%%MC2%%% *subroutine out_sfilec (stepno) 2,2 implicit none * integer stepno * #include "lcldim.cdk"
#include "cdate.cdk"
#include "partopo.cdk"
#include "yomdyn1.cdk"
#include "grd.cdk"
#include "sor.cdk"
#include "rec.cdk"
#include "levels.cdk"
* integer longueur external longueur character*15 datev,startindx character*512 filen integer i,k,nis,njs,ind_o(gnk+1) real*8 dayfrac,sec_in_day parameter ( sec_in_day=86400.0d0 ) * *------------------------------------------------------------------ * if (out_unfc.gt.0) return * out_unfc = 79 if (blocme.eq.0) then dayfrac = dble(stepno) * dble(grdt) / sec_in_day call incdatsd
(datev,gcrunstrt,dayfrac) write (startindx,'((i7.7),a1,(i7.7))') out_idg,'-',out_jdg filen='3df_'//datev//'_'//startindx c filen= '../../output/casc_n/'//filen(1:longueur(filen)) filen= '../../output/casc/'//filen(1:longueur(filen)) open (out_unfc,file=filen,access='SEQUENTIAL', $ form='UNFORMATTED') call wrgeocasc
( xpx(g_id+hx-1),out_nisg+1,ypx(g_jd+hy-1), $ out_njsg+1,zt,htop,iscal,gnk,nvardyn,out_unfc ) endif * *------------------------------------------------------------------ return end *
subroutine wrgeocasc ( xp,ni,yp,nj,zt,htop,iscal,nk,ntr,unf ) 1 implicit none * integer ni,nj,nk,iscal(2),ntr,unf real zt(nk), htop real*8 xp(ni), yp(nj) * *------------------------------------------------------------------ write (unf) '>>^^',ni,nj,ntr write (unf) xp,yp if (nk.gt.0) then write (unf) 'HTHT',nk,iscal,htop write (unf) zt endif * *------------------------------------------------------------------ return end