Ignore:
Timestamp:
2011-09-26T10:54:50+02:00 (9 years ago)
Author:
cbricaud
Message:

cleanning, minor modifications

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/diadct_sections.f90

    r2849 r2858  
    77     !! 
    88     !! 
    9      !! History: 2011: cbricaud Mercator-Ocean 
     9     !! History: 09/2011: Clement Bricaud ( Mercator-Ocean ) 
    1010     !! 
    1111     !!============================================================================== 
     
    3434     PRINT*,'CREATION OF SECTIONS FOR NEMO diadct.F90 ROUTINE' 
    3535     PRINT*,'~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    36   
    37      !----------------------! 
    38      !0. Read arguments     ! 
    39      !----------------------! 
    40      PRINT*,'              ' 
    41      PRINT*,'READ ARGUMENTS' 
    42      PRINT*,'--------------' 
    43   
    44      !check number of arguments and display usage message if wrong 
    45      narg=iargc() 
    46      PRINT*,'narg= ',narg 
    47      IF ( narg /= 4 ) THEN 
    48           PRINT *,' Usage : generate_sections jpidta jpjdta jpizoom jpjzoom ' 
    49           STOP 
    50      ENDIF 
    51  
    52      ! read on line arguments 
    53      CALL getarg(1,cdum) ; READ(cdum,*) jpidta 
    54      CALL getarg(2,cdum) ; READ(cdum,*) jpjdta 
    55      CALL getarg(3,cdum) ; READ(cdum,*) jpizoom  
    56      CALL getarg(4,cdum) ; READ(cdum,*) jpjzoom 
    57    
    58      PRINT*,'jpidta  jpjdta =',jpidta,jpjdta 
    59      PRINT*,'jpizoom jpjzoom=',jpizoom,jpjzoom 
    60   
    61      !------------------!  
    62      !0. INITIALISATION !       
    63      !------------------!      
    64      PRINT*,'              ' 
    65      PRINT*,'DOMAIN SIZE' 
    66      PRINT*,'--------------' 
    67  
    68      !Domain size 
    69      jpiglo = jpidta-jpizoom+1  ; jpjglo = jpjdta-jpjzoom+1 
    70      jpi    = jpiglo            ; jpj    = jpjglo 
    71      nlci   = jpiglo            ; nlcj   = jpjglo  
    72      nlei   = jpiglo            ; nlej   = jpjglo 
    73  
    74      PRINT*,'jpiglo jpjglo = ',jpiglo,jpjglo 
    75      PRINT*,'jpi    jpj    = ',jpi   ,jpj  
    76      PRINT*,'nlci  nlcj    = ',nlci,nlcj 
    7736  
    7837     !-------------------! 
     
    10867     CALL read_coord_mesh 
    10968  
     69     PRINT*,'domain sizes: ' 
     70     PRINT*,'jpi    jpj    = ',jpi   ,jpj  
    11071     PRINT*,'domain boundaries: ' 
    11172     PRINT*,' 1   1   ',glamt(1,1),gphit(1,1) 
     
    13596     ENDDO 
    13697 
    137      !----------------------! 
    138      !5.ecriture du fichier ! 
    139      !----------------------! 
     98     !--------------------------------! 
     99     !5.Write section_ijglobal.diadct ! 
     100     !--------------------------------! 
    140101     CALL write_sections 
    141102 
Note: See TracChangeset for help on using the changeset viewer.