source: LMDZ.3.3/trunk/libf/dyn3d/read_dstoke.F @ 2793

c $Header$
            subroutine read_dstoke(irec,zdtvr,ziadvtrac,ziadvtrac2)


            IMPLICIT NONE

#include "netcdf.inc"
#include "dimensions.h"
#include "paramet.h"
#include "comgeom.h"
#include "comvert.h"

        
        integer ncidfd
        save ncidfd
        integer variddt,varididvt,varididvp
        save variddt,varididvt,varididvp
        real*4 dtv(1,1,1),adv2(1,1,1),adv(1,1,1)
        real zdtvr,ziadvtrac,ziadvtrac2
        integer status,irec
        real rcode
        INTEGER epais(2), debut(2)
        
        if (irec .eq. 0) then
        
        ncidfd=NCOPN('defstoke.nc',NCNOWRIT,rcode)

        variddt=NCVID(ncidfd,'dtvr',rcode)
        print*,'ncidfd,variddt',ncidfd,variddt

        varididvt=NCVID(ncidfd,'istdyn',rcode)
        print*,'ncidfd,varididvt',ncidfd,varididvt
        
        varididvp=NCVID(ncidfd,'istphy',rcode)
        print*,'ncidfd,varididvp',ncidfd,varididvp

c lecture de zdtvr et ziadvtrac
        
        epais(1) = 1
        epais(2) = 1
        debut(1) = 1
        debut(2) = 1

#ifdef NC_DOUBLE
        status=NF_GET_VARA_DOUBLE(ncidfd,variddt,debut,epais,dtv)
#else
        status=NF_GET_VARA_REAL(ncidfd,variddt,debut,epais,dtv)
#endif
        zdtvr=dtv(1,1,1)

#ifdef NC_DOUBLE
        status=NF_GET_VARA_DOUBLE(ncidfd,varididvt,debut,epais,adv)
#else
        status=NF_GET_VARA_REAL(ncidfd,varididvt,debut,epais,adv)
#endif
        ziadvtrac= adv(1,1,1)
        
#ifdef NC_DOUBLE
        status=NF_GET_VARA_DOUBLE(ncidfd,varididvp,debut,epais,adv2)
#else
        status=NF_GET_VARA_REAL(ncidfd,varididvp,debut,epais,adv2)
#endif
        ziadvtrac2= adv2(1,1,1)

        write(*,*) 'ds read_dstoke zdtvr = ', zdtvr
        write(*,*) 'ds read_dstoke ziadvtrac = ', ziadvtrac
        write(*,*) 'ds read_dstoke ziadvtrac2 = ', ziadvtrac2
        
c           status=NF_CLOSE(ncidfd)

        else
            stop'Pas bon irec ne 0'
        endif
        
        return

        end