| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119 |
- # 0 "<stdin>"
- # 0 "<built-in>"
- # 0 "<command-line>"
- # 1 "/usr/include/stdc-predef.h" 1 3 4
- # 17 "/usr/include/stdc-predef.h" 3 4
- # 2 "<command-line>" 2
- # 1 "<stdin>"
- # 10 "<stdin>"
- module m_put_mod_fld
- ! KAL -- This routine puts one of the fields to the restart file, specified
- ! KAL -- by name, vertical level and time level.
- ! KAL -- Its a bit dangerous to use -- indx must be updated correctly (max one
- ! KAL -- increment per call of this routine), otherwise there wil be a
- ! KAL -- inconsistency between .a and .b files
- contains
- !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- ! KAL - This is for the new file type
- subroutine put_mod_fld(memfile,fld,iens,cfld,vlevel,tlevel,indx,nx,ny)
- use mod_raw_io
- implicit none
- integer, intent(in) :: nx,ny
- integer, intent(in) :: iens ! Ensemble member to read
- real, dimension(nx,ny), intent(in) :: fld ! output fld
- character(len=*), intent(in) :: memfile! base name of input files
- character(len=8), intent(in) :: cfld ! name of fld
- integer, intent(in) :: tlevel ! time level
- integer, intent(in) :: vlevel ! vertical level
- integer, intent(in) :: indx ! index into file
- real*4:: amin, amax,spval
- real*4:: writefldr4(nx,ny)
- integer , parameter :: nop=123
- integer :: ios
- writefldr4 = real(fld, 4)
- ! Write fld into index of .a file -- Direct access file
- call writeraw(writefldr4, &! feltet som skal skrivast
- amin,amax, &! min og max (returnerast fra writeraw)
- nx,ny, &! dimensjon
- .false.,spval, &! feltet vil ikkje faa "no-value" verdiar
- memfile//'.a', &! forelopig filnavn...
- indx) ! indexen i .a fila
- ! Skriv header -- .b fil -- Hold tunga rett i munnen her og utanfor rutina,
- ! ellers blir det inkonsistens mellom .a og .b filer - dette vil fangast opp
- ! av postprosessering. Dette er tungvint men vanskelig aa omgaa fordi
- ! .b-fila er sekvensiell mens .a fila er direct access.
- if (indx==1) then
- ! forste indeks - vi overskriv evt gamle filer
- open(nop,file=memfile//'.b',status='replace')
- else
- ! Ellers legg vi til
- open(nop,file=memfile//'.b',status='old', position='append')
- end if
- ! Skriv i vei !
- write(nop,4100,iostat=ios) cfld,vlevel,tlevel,amin,amax
- close(nop)
- 4100 format(a,': layer,tlevel,range = ',i3,i3,2x,1p2e16.7)
- ! format (a8,23x,i3,i3,2x,2e16.7)
- end subroutine
- end module m_put_mod_fld
|
