Ignore:
Timestamp:
2019-03-29T13:54:25+01:00 (3 years ago)
Author:
smueller
Message:

Replacement of the module variable used to store information about all available tidal components (variable Wave in module tide_mod) by an array used to store information about the selected components only (variable tide_components in module tide_mod), replacement of the corresponding initialisation subroutine, as well as related adjustments in various modules (bdytides, diaharm, sbctide, and tide_mod) and in one include file (tide.h90) (ticket #2194)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/src/OCE/BDY/bdytides.F90

    r10800 r10811  
    115115                    WRITE(numout,*) '             Tidal components: '  
    116116               DO itide = 1, nb_harmo 
    117                   WRITE(numout,*)  '                 ', Wave(ntide(itide))%cname_tide  
     117                  WRITE(numout,*)  '                 ', tide_components(itide)%cname_tide  
    118118               END DO 
    119119            ENDIF  
     
    156156               igrd = 1                       ! Everything is at T-points here 
    157157               DO itide = 1, nb_harmo 
    158                   CALL iom_get( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_z1', ztr(:,:) ) 
    159                   CALL iom_get( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_z2', zti(:,:) )  
     158                  CALL iom_get( inum, jpdom_autoglo, TRIM(tide_components(itide)%cname_tide)//'_z1', ztr(:,:) ) 
     159                  CALL iom_get( inum, jpdom_autoglo, TRIM(tide_components(itide)%cname_tide)//'_z2', zti(:,:) )  
    160160                  DO ib = 1, ilen0(igrd) 
    161161                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    172172               igrd = 2                       ! Everything is at U-points here 
    173173               DO itide = 1, nb_harmo 
    174                   CALL iom_get  ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_u1', ztr(:,:) ) 
    175                   CALL iom_get  ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_u2', zti(:,:) ) 
     174                  CALL iom_get  ( inum, jpdom_autoglo, TRIM(tide_components(itide)%cname_tide)//'_u1', ztr(:,:) ) 
     175                  CALL iom_get  ( inum, jpdom_autoglo, TRIM(tide_components(itide)%cname_tide)//'_u2', zti(:,:) ) 
    176176                  DO ib = 1, ilen0(igrd) 
    177177                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    188188               igrd = 3                       ! Everything is at V-points here 
    189189               DO itide = 1, nb_harmo 
    190                   CALL iom_get  ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_v1', ztr(:,:) ) 
    191                   CALL iom_get  ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_v2', zti(:,:) ) 
     190                  CALL iom_get  ( inum, jpdom_autoglo, TRIM(tide_components(itide)%cname_tide)//'_v1', ztr(:,:) ) 
     191                  CALL iom_get  ( inum, jpdom_autoglo, TRIM(tide_components(itide)%cname_tide)//'_v2', zti(:,:) ) 
    192192                  DO ib = 1, ilen0(igrd) 
    193193                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    217217               DO itide = 1, nb_harmo 
    218218                  !                                                              ! SSH fields 
    219                   clfile = TRIM(filtide)//TRIM(Wave(ntide(itide))%cname_tide)//'_grid_T.nc' 
     219                  clfile = TRIM(filtide)//TRIM(tide_components(itide)%cname_tide)//'_grid_T.nc' 
    220220                  CALL iom_open( clfile, inum ) 
    221221                  CALL fld_map( inum, 'z1' , dta_read(1:ilen0(1),1:1,1:1) , 1,  ibmap_ptr(1) ) 
     
    225225                  CALL iom_close( inum ) 
    226226                  !                                                              ! U fields 
    227                   clfile = TRIM(filtide)//TRIM(Wave(ntide(itide))%cname_tide)//'_grid_U.nc' 
     227                  clfile = TRIM(filtide)//TRIM(tide_components(itide)%cname_tide)//'_grid_U.nc' 
    228228                  CALL iom_open( clfile, inum ) 
    229229                  CALL fld_map( inum, 'u1' , dta_read(1:ilen0(2),1:1,1:1) , 1, ibmap_ptr(2) ) 
     
    233233                  CALL iom_close( inum ) 
    234234                  !                                                              ! V fields 
    235                   clfile = TRIM(filtide)//TRIM(Wave(ntide(itide))%cname_tide)//'_grid_V.nc' 
     235                  clfile = TRIM(filtide)//TRIM(tide_components(itide)%cname_tide)//'_grid_V.nc' 
    236236                  CALL iom_open( clfile, inum ) 
    237237                  CALL fld_map( inum, 'v1' , dta_read(1:ilen0(3),1:1,1:1) , 1, ibmap_ptr(3) ) 
Note: See TracChangeset for help on using the changeset viewer.