Ignore:
Timestamp:
Apr 13, 2001, 12:44:53 PM (23 years ago)
Author:
lmdz
Message:

Debogage du guidage et de la version debranchee et abandon de la version
debranchee non-netcdf FH/MAF
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F

    r198 r204  
     1c
    12c $Header
     3c
    24      subroutine read_reanalyse(timestep,u,v,t,masse,ps,mode)
    35
     
    1416      integer nlevnc
    1517cModef 11-2-99      parameter (nlevnc=15)
    16       parameter (nlevnc=15)
     18      parameter (nlevnc=21)
    1719      integer timestep,mode,l
    1820
     
    5860            varidt=NCVID(ncidt,'AIR',rcode)
    5961            print*,'ncidt,varidt',ncidt,varidt
    60             ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)
    61             varidps=NCVID(ncidps,'SP',rcode)
    62             print*,'ncidps,varidps',ncidps,varidps
     62c           ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)
     63c           varidps=NCVID(ncidps,'SP',rcode)
     64c           print*,'ncidps,varidps',ncidps,varidps
    6365      endif
    6466
     
    106108      count(4)=0
    107109c  ps
    108       status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc)
     110c     status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc)
    109111c     call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ')
    110       call correctbid(iim,jjp1,psnc)
     112c     call correctbid(iim,jjp1,psnc)
    111113
    112114c   Transformations
     
    187189c -----------------------------------------------------------------
    188190
    189       do j=1,jjp1
    190          do i=1,iim
    191             ps(i,j)=psnc(i,jjp1+1-j)
    192          enddo
    193          ps(iip1,j)=ps(1,j)
    194       enddo
     191c     do j=1,jjp1
     192c        do i=1,iim
     193c           ps(i,j)=psnc(i,jjp1+1-j)
     194c        enddo
     195c        ps(iip1,j)=ps(1,j)
     196c     enddo
    195197
    196198      CALL pression( ip1jmp1, ap, bp, ps, p )
     
    200202      unskap=1./kappa
    201203      prefkap =  preff  ** kappa
    202       PRINT *,' Pref kappa unskap  ',preff,kappa,unskap
     204c     PRINT *,' Pref kappa unskap  ',preff,kappa,unskap
    203205      DO l = 1, llm
    204206       DO j=1,jjp1
Note: See TracChangeset for help on using the changeset viewer.