| 1 | subroutine wrgrads(if,nl,field,name,titlevar) |
|---|
| 2 | implicit none |
|---|
| 3 | |
|---|
| 4 | c Declarations |
|---|
| 5 | c if indice du fichier |
|---|
| 6 | c nl nombre de couches |
|---|
| 7 | c field champ |
|---|
| 8 | c name petit nom |
|---|
| 9 | c titlevar Titre |
|---|
| 10 | |
|---|
| 11 | #include "gradsdef.h" |
|---|
| 12 | |
|---|
| 13 | c arguments |
|---|
| 14 | integer if,nl |
|---|
| 15 | real field(imx*jmx*lmx) |
|---|
| 16 | character*10 name,file |
|---|
| 17 | character*10 titlevar |
|---|
| 18 | |
|---|
| 19 | c local |
|---|
| 20 | |
|---|
| 21 | integer im,jm,lm,i,j,l,lnblnk,iv,iii,iji,iif,ijf |
|---|
| 22 | |
|---|
| 23 | logical writectl |
|---|
| 24 | |
|---|
| 25 | |
|---|
| 26 | writectl=.false. |
|---|
| 27 | |
|---|
| 28 | c print*,if,iid(if),jid(if),ifd(if),jfd(if) |
|---|
| 29 | iii=iid(if) |
|---|
| 30 | iji=jid(if) |
|---|
| 31 | iif=ifd(if) |
|---|
| 32 | ijf=jfd(if) |
|---|
| 33 | im=iif-iii+1 |
|---|
| 34 | jm=ijf-iji+1 |
|---|
| 35 | lm=lmd(if) |
|---|
| 36 | |
|---|
| 37 | c print*,'im,jm,lm,name,firsttime(if)' |
|---|
| 38 | c print*,im,jm,lm,name,firsttime(if) |
|---|
| 39 | |
|---|
| 40 | if(firsttime(if)) then |
|---|
| 41 | if(name.eq.var(1,if)) then |
|---|
| 42 | firsttime(if)=.false. |
|---|
| 43 | ivar(if)=1 |
|---|
| 44 | print*,'fin de l initialiation de l ecriture du fichier' |
|---|
| 45 | print*,file |
|---|
| 46 | print*,'fichier no: ',if |
|---|
| 47 | print*,'unit ',unit(if) |
|---|
| 48 | print*,'nvar ',nvar(if) |
|---|
| 49 | print*,'vars ',(var(iv,if),iv=1,nvar(if)) |
|---|
| 50 | else |
|---|
| 51 | ivar(if)=ivar(if)+1 |
|---|
| 52 | nvar(if)=ivar(if) |
|---|
| 53 | var(ivar(if),if)=name |
|---|
| 54 | tvar(ivar(if),if)=titlevar(1:lnblnk(titlevar)) |
|---|
| 55 | nld(ivar(if),if)=nl |
|---|
| 56 | print*,'initialisation ecriture de ',var(ivar(if),if) |
|---|
| 57 | print*,'if ivar(if) nld ',if,ivar(if),nld(ivar(if),if) |
|---|
| 58 | endif |
|---|
| 59 | writectl=.true. |
|---|
| 60 | itime(if)=1 |
|---|
| 61 | else |
|---|
| 62 | ivar(if)=mod(ivar(if),nvar(if))+1 |
|---|
| 63 | if (ivar(if).eq.nvar(if)) then |
|---|
| 64 | writectl=.true. |
|---|
| 65 | itime(if)=itime(if)+1 |
|---|
| 66 | endif |
|---|
| 67 | |
|---|
| 68 | if(var(ivar(if),if).ne.name) then |
|---|
| 69 | print*,'Il faut stoker la meme succession de champs a chaque' |
|---|
| 70 | print*,'pas de temps' |
|---|
| 71 | print*,'fichier no: ',if |
|---|
| 72 | print*,'unit ',unit(if) |
|---|
| 73 | print*,'nvar ',nvar(if) |
|---|
| 74 | print*,'vars ',(var(iv,if),iv=1,nvar(if)) |
|---|
| 75 | |
|---|
| 76 | stop |
|---|
| 77 | endif |
|---|
| 78 | endif |
|---|
| 79 | |
|---|
| 80 | c print*,'ivar(if),nvar(if),var(ivar(if),if),writectl' |
|---|
| 81 | c print*,ivar(if),nvar(if),var(ivar(if),if),writectl |
|---|
| 82 | do l=1,nl |
|---|
| 83 | irec(if)=irec(if)+1 |
|---|
| 84 | c print*,'Ecrit rec=',irec(if),iii,iif,iji,ijf, |
|---|
| 85 | c s (l-1)*imd(if)*jmd(if)+(iji-1)*imd(if)+iii |
|---|
| 86 | c s ,(l-1)*imd(if)*jmd(if)+(ijf-1)*imd(if)+iif |
|---|
| 87 | write(unit(if)+1,rec=irec(if)) |
|---|
| 88 | s ((field((l-1)*imd(if)*jmd(if)+(j-1)*imd(if)+i) |
|---|
| 89 | s ,i=iii,iif),j=iji,ijf) |
|---|
| 90 | enddo |
|---|
| 91 | if (writectl) then |
|---|
| 92 | |
|---|
| 93 | file=fichier(if) |
|---|
| 94 | c WARNING! on reecrase le fichier .ctl a chaque ecriture |
|---|
| 95 | open(unit(if),file=file(1:lnblnk(file))//'.ctl', |
|---|
| 96 | & form='formatted',status='unknown') |
|---|
| 97 | write(unit(if),'(a5,1x,a40)') |
|---|
| 98 | & 'DSET ','^'//file(1:lnblnk(file))//'.dat' |
|---|
| 99 | |
|---|
| 100 | write(unit(if),'(a12)') 'UNDEF 1.0E30' |
|---|
| 101 | write(unit(if),'(a5,1x,a40)') 'TITLE ',title(if) |
|---|
| 102 | call formcoord(unit(if),im,xd(iii,if),1.,.false.,'XDEF') |
|---|
| 103 | call formcoord(unit(if),jm,yd(iji,if),1.,.true.,'YDEF') |
|---|
| 104 | call formcoord(unit(if),lm,zd(1,if),1.,.false.,'ZDEF') |
|---|
| 105 | write(unit(if),'(a4,i10,a30)') |
|---|
| 106 | & 'TDEF ',itime(if),' LINEAR 02JAN1987 1MO ' |
|---|
| 107 | write(unit(if),'(a4,2x,i5)') 'VARS',nvar(if) |
|---|
| 108 | do iv=1,nvar(if) |
|---|
| 109 | c print*,'if,var(iv,if),nld(iv,if),nld(iv,if)-1/nld(iv,if)' |
|---|
| 110 | c print*,if,var(iv,if),nld(iv,if),nld(iv,if)-1/nld(iv,if) |
|---|
| 111 | write(unit(if),1000) var(iv,if),nld(iv,if)-1/nld(iv,if) |
|---|
| 112 | & ,99,tvar(iv,if) |
|---|
| 113 | enddo |
|---|
| 114 | write(unit(if),'(a7)') 'ENDVARS' |
|---|
| 115 | c |
|---|
| 116 | 1000 format(a5,3x,i4,i3,1x,a39) |
|---|
| 117 | |
|---|
| 118 | close(unit(if)) |
|---|
| 119 | |
|---|
| 120 | endif ! writectl |
|---|
| 121 | |
|---|
| 122 | return |
|---|
| 123 | |
|---|
| 124 | END |
|---|
| 125 | |
|---|
