/[MITgcm]/MITgcm/pkg/mnc/mnc_cw_readwrite.template
ViewVC logotype

Diff of /MITgcm/pkg/mnc/mnc_cw_readwrite.template

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.32 by edhill, Sat Sep 10 18:30:07 2005 UTC revision 1.33 by edhill, Mon Sep 19 02:24:40 2005 UTC
# Line 776  C           There is only one face Line 776  C           There is only one face
776            ENDIF            ENDIF
777    
778  C         Try to open the PER-FACE file  C         Try to open the PER-FACE file
779    C         WRITE(*,*) 'trying: "', path_fname, '"'
780            err = NF_OPEN(path_fname, NF_NOWRITE, fid)            err = NF_OPEN(path_fname, NF_NOWRITE, fid)
781            IF ( err .EQ. NF_NOERR ) THEN            IF ( err .EQ. NF_NOERR ) THEN
782              f_or_t = 0              f_or_t = 1
783            ELSE            ELSE
784    
785  C           Create the PER-TILE file name  C           Create the PER-TILE file name
# Line 788  C           Create the PER-TILE file nam Line 789  C           Create the PER-TILE file nam
789              WRITE(path_fname,'(2a,a1,a,a3)')              WRITE(path_fname,'(2a,a1,a,a3)')
790       &           mnc_indir_str(1:npath), fname(1:ntot), 't',       &           mnc_indir_str(1:npath), fname(1:ntot), 't',
791       &           tmpnm(1:k), '.nc'       &           tmpnm(1:k), '.nc'
792    C           WRITE(*,*) 'trying: "', path_fname, '"'
793              err = NF_OPEN(path_fname, NF_NOWRITE, fid)              err = NF_OPEN(path_fname, NF_NOWRITE, fid)
794              IF ( err .EQ. NF_NOERR ) THEN              IF ( err .EQ. NF_NOERR ) THEN
795                f_or_t = 1                f_or_t = 0
796              ELSE              ELSE
797                k = ILNBLNK(path_fname)                k = ILNBLNK(path_fname)
798                write(msgbuf,'(4a)')                write(msgbuf,'(4a)')
# Line 811  C           Create the PER-TILE file nam Line 813  C           Create the PER-TILE file nam
813            err = NF_INQ_VARID(fid, vtype, idv)            err = NF_INQ_VARID(fid, vtype, idv)
814            CALL MNC_HANDLE_ERR(err, msgbuf, myThid)            CALL MNC_HANDLE_ERR(err, msgbuf, myThid)
815    
816              k = ILNBLNK(path_fname)
817              fname(1:k) = path_fname(1:k)
818              nfname = k
819    
820  #endif  #endif
821    
822            IF ( f_or_t .EQ. 1 ) THEN            IF ( f_or_t .EQ. 1 ) THEN
             write(msgbuf,'(4a)')  
      &           'MNC_CW_RX_R: per-face reads are not yet ',  
      &           'implemented -- so pester Ed to finish them'  
             CALL print_error(msgbuf, mythid)  
             STOP 'ABNORMAL END: S/R MNC_CW_RX_W'  
823    
824    C           write(msgbuf,'(2a)')
825    C           &           'MNC_CW_RX_R: per-face reads are not yet ',
826    C           &           'implemented -- so pester Ed to finish them'
827    C           CALL print_error(msgbuf, mythid)
828    C           STOP 'ABNORMAL END: S/R MNC_CW_RX_W'
829                
830  C           Get the X,Y PER-FACE offsets  C           Get the X,Y PER-FACE offsets
831              CALL MNC_CW_GET_XYFO(lbi,lbj, ixoff,iyoff, myThid)              CALL MNC_CW_GET_XYFO(lbi,lbj, ixoff,iyoff, myThid)
832    
833            ENDIF            ENDIF
834    
835    C         WRITE(*,*) 'f_or_t = ',f_or_t
836    
837  C         Check that the current sNy,sNy values and the in-file values  C         Check that the current sNy,sNy values and the in-file values
838  C         are compatible and WARN (only warn) if not  C         are compatible and WARN (only warn) if not
839  C           f_sNx = -1  C           f_sNx = -1
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22