Ignore:
Timestamp:
2019-05-22T22:44:28+02:00 (17 months ago)
Author:
rrenshaw
Message:

code fixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis2/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r10974 r11038  
    7575 
    7676  !! * Module variables 
    77   INTEGER :: nn_dct        ! Frequency of computation 
    78   INTEGER :: nn_dctwri     ! Frequency of output 
    79   INTEGER :: nn_secdebug   ! Number of the section to debug 
    80   INTEGER :: nn_dct_h      ! Frequency of computation for NOOS hourly files 
    81   INTEGER :: nn_dctwri_h   ! Frequency of output for NOOS hourly files 
     77  INTEGER :: nn_dct      = 1     ! Frequency of computation 
     78  INTEGER :: nn_dctwri   = 1     ! Frequency of output 
     79  INTEGER :: nn_secdebug = 0     ! Number of the section to debug 
     80  INTEGER :: nn_dct_h    = 1     ! Frequency of computation for NOOS hourly files 
     81  INTEGER :: nn_dctwri_h = 1     ! Frequency of output for NOOS hourly files 
    8282    
    8383  INTEGER, PARAMETER :: nb_class_max  = 12   ! maximum number of classes, i.e. depth levels or density classes 
     
    209209        IF( ln_NOOS ) THEN 
    210210           WRITE(numout,*) "       Calculate NOOS hourly output: ln_dct_calc_noos_hr = ",ln_dct_calc_noos_hr 
    211            WRITE(numout,*) "       Calculate NOOS 25 hour mean output: ln_dct_calc_noos_hr = ",ln_dct_calc_noos_25h 
     211           WRITE(numout,*) "       Calculate NOOS 25 hour mean output: ln_dct_calc_noos_25h = ",ln_dct_calc_noos_25h 
    212212           WRITE(numout,*) "       Use IOM Output: ln_dct_iom_cont = ",ln_dct_iom_cont 
    213213           WRITE(numout,*) "       Output in ASCII (True) or Binary (False): ln_dct_ascii = ",ln_dct_ascii 
     
    254254           ENDIF 
    255255        ELSE 
    256            CALL ctl_opn( numdct_vol,  'volume_transport', 'NEW', 'FORMATTED', 'SEQUENTIAL', -1, numout,  .FALSE. ) 
    257            CALL ctl_opn( numdct_heat, 'heat_transport'  , 'NEW', 'FORMATTED', 'SEQUENTIAL', -1, numout,  .FALSE. ) 
    258            CALL ctl_opn( numdct_salt, 'salt_transport'  , 'NEW', 'FORMATTED', 'SEQUENTIAL', -1, numout,  .FALSE. ) 
     256           CALL ctl_opn( numdct_vol,  'volume_transport', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout,  .FALSE. ) 
     257           CALL ctl_opn( numdct_heat, 'heat_transport'  , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout,  .FALSE. ) 
     258           CALL ctl_opn( numdct_salt, 'salt_transport'  , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout,  .FALSE. ) 
    259259        ENDIF 
    260260     ENDIF 
     
    559559     CALL wrk_alloc( nb_point_max, directemp ) 
    560560 
     561! --------------------------------------------------------------------------------------------- 
     562! This section seems to be needed, else later reads find no sections in section_ijglobal.diadct 
    561563     !open input file 
    562564     !--------------- 
     
    575577     READ(107) isec 
    576578     CLOSE(107) 
     579! --------------------------------------------------------------------------------------------- 
    577580      
    578581     CALL ctl_opn( numdct_in, 'section_ijglobal.diadct', 'OLD', 'UNFORMATTED', 'SEQUENTIAL', -1, numout, .TRUE. ) 
Note: See TracChangeset for help on using the changeset viewer.