New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 13767 – NEMO

Changeset 13767


Ignore:
Timestamp:
2020-11-10T16:42:30+01:00 (4 years ago)
Author:
dford
Message:

Further bug fixes for compilation.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13747_ENHANCE-04_dford_OBSOP_BGC/src/OCE/OBS/diaobs.F90

    r13758 r13767  
    5757   PUBLIC calc_date        ! Compute the date of a timestep 
    5858 
    59    LOGICAL, PUBLIC :: ln_diaobs          !: Logical switch for the obs operator 
    60    LOGICAL         :: ln_sstnight        !  Logical switch for night mean SST obs 
    61    LOGICAL         :: ln_sla_fp_indegs   !  T=> SLA obs footprint size specified in degrees, F=> in metres 
    62    LOGICAL         :: ln_sst_fp_indegs   !  T=> SST obs footprint size specified in degrees, F=> in metres 
    63    LOGICAL         :: ln_sss_fp_indegs   !  T=> SSS obs footprint size specified in degrees, F=> in metres 
    64    LOGICAL         :: ln_sic_fp_indegs   !  T=> sea-ice obs footprint size specified in degrees, F=> in metres 
    65  
    66    REAL(wp) ::   rn_default_avglamscl    ! E/W diameter of SLA observation footprint (metres) 
    67    REAL(wp) ::   rn_default_avgphiscl    ! N/S diameter of SLA observation footprint (metre 
    68    REAL(wp) ::   rn_sla_avglamscl        ! E/W diameter of SLA observation footprint (metres) 
    69    REAL(wp) ::   rn_sla_avgphiscl        ! N/S diameter of SLA observation footprint (metres) 
    70    REAL(wp) ::   rn_sst_avglamscl        ! E/W diameter of SST observation footprint (metres) 
    71    REAL(wp) ::   rn_sst_avgphiscl        ! N/S diameter of SST observation footprint (metres) 
    72    REAL(wp) ::   rn_sss_avglamscl        ! E/W diameter of SSS observation footprint (metres) 
    73    REAL(wp) ::   rn_sss_avgphiscl        ! N/S diameter of SSS observation footprint (metres) 
    74    REAL(wp) ::   rn_sic_avglamscl        ! E/W diameter of sea-ice observation footprint (metres) 
    75    REAL(wp) ::   rn_sic_avgphiscl        ! N/S diameter of sea-ice observation footprint (metres) 
    76  
    77    INTEGER :: nn_1dint                   ! Vertical interpolation method 
    78    INTEGER :: nn_2dint_default           ! Default horizontal interpolation method 
    79    INTEGER :: nn_2dint_sla               ! SLA horizontal interpolation method  
    80    INTEGER :: nn_2dint_sst               ! SST horizontal interpolation method  
    81    INTEGER :: nn_2dint_sss               ! SSS horizontal interpolation method  
    82    INTEGER :: nn_2dint_sic               ! Seaice horizontal interpolation method  
     59   LOGICAL, PUBLIC :: ln_diaobs            !: Logical switch for the obs operator 
     60   LOGICAL         :: ln_sstnight          !  Logical switch for night mean SST obs 
     61   LOGICAL         :: ln_default_fp_indegs !  T=> Default obs footprint size specified in degrees, F=> in metres 
     62   LOGICAL         :: ln_sla_fp_indegs     !  T=> SLA obs footprint size specified in degrees, F=> in metres 
     63   LOGICAL         :: ln_sst_fp_indegs     !  T=> SST obs footprint size specified in degrees, F=> in metres 
     64   LOGICAL         :: ln_sss_fp_indegs     !  T=> SSS obs footprint size specified in degrees, F=> in metres 
     65   LOGICAL         :: ln_sic_fp_indegs     !  T=> sea-ice obs footprint size specified in degrees, F=> in metres 
     66 
     67   REAL(wp) ::   rn_default_avglamscl      ! E/W diameter of SLA observation footprint (metres) 
     68   REAL(wp) ::   rn_default_avgphiscl      ! N/S diameter of SLA observation footprint (metre 
     69   REAL(wp) ::   rn_sla_avglamscl          ! E/W diameter of SLA observation footprint (metres) 
     70   REAL(wp) ::   rn_sla_avgphiscl          ! N/S diameter of SLA observation footprint (metres) 
     71   REAL(wp) ::   rn_sst_avglamscl          ! E/W diameter of SST observation footprint (metres) 
     72   REAL(wp) ::   rn_sst_avgphiscl          ! N/S diameter of SST observation footprint (metres) 
     73   REAL(wp) ::   rn_sss_avglamscl          ! E/W diameter of SSS observation footprint (metres) 
     74   REAL(wp) ::   rn_sss_avgphiscl          ! N/S diameter of SSS observation footprint (metres) 
     75   REAL(wp) ::   rn_sic_avglamscl          ! E/W diameter of sea-ice observation footprint (metres) 
     76   REAL(wp) ::   rn_sic_avgphiscl          ! N/S diameter of sea-ice observation footprint (metres) 
     77 
     78   INTEGER :: nn_1dint                     ! Vertical interpolation method 
     79   INTEGER :: nn_2dint_default             ! Default horizontal interpolation method 
     80   INTEGER :: nn_2dint_sla                 ! SLA horizontal interpolation method  
     81   INTEGER :: nn_2dint_sst                 ! SST horizontal interpolation method  
     82   INTEGER :: nn_2dint_sss                 ! SSS horizontal interpolation method  
     83   INTEGER :: nn_2dint_sic                 ! Seaice horizontal interpolation method  
    8384   INTEGER, DIMENSION(imaxavtypes) ::   nn_profdavtypes   ! Profile data types representing a daily average 
    8485   INTEGER :: nproftypes     ! Number of profile obs types 
     
    133134         & cn_sicfbfiles, &      ! Seaice concentration input filenames 
    134135         & cn_velfbfiles, &      ! Velocity profile input filenames 
    135          & cn_sstbiasfiles      ! SST bias input filenames 
     136         & cn_sstbiasfiles       ! SST bias input filenames 
    136137      CHARACTER(LEN=128) :: & 
    137138         & cn_altbiasfile        ! Altimeter bias input filename 
     
    139140         & clproffiles, &        ! Profile filenames 
    140141         & clsurffiles           ! Surface filenames 
     142      CHARACTER(len=8), DIMENSION(:), ALLOCATABLE :: & 
     143         & clvars                ! Expected variable names 
    141144         ! 
    142145      LOGICAL :: ln_t3d          ! Logical switch for temperature profiles 
     
    452455            IF ( TRIM(cobstypesprof(jtype)) == 'prof' ) THEN 
    453456               nvarsprof(jtype) = 2 
    454                nextrprof(jtype) = 1 
    455                IF ( ln_output_clim ) ltype_clim = .TRUE.               
     457               nextrprof(jtype) = 1              
    456458               ALLOCATE( llvar (nvarsprof(jtype)) ) 
    457459               ALLOCATE( clvars(nvarsprof(jtype)) ) 
     
    556558               &               clsurffiles(jtype,1:ifilessurf(jtype)), & 
    557559               &               nvarssurf(jtype), nextrsurf(jtype), nitend-nit000+2, & 
    558                &               rn_dobsini, rn_dobsend, ln_ignmis, .FALSE., llnightav(jtype), * 
     560               &               rn_dobsini, rn_dobsend, ln_ignmis, .FALSE., llnightav(jtype), & 
    559561               &               clvars ) 
    560562               ! 
     
    623625      INTEGER :: ji, jj            ! Loop counters 
    624626      REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: & 
    625          & zprofvar,               ! Model values for variables in a prof ob 
     627         & zprofvar                ! Model values for variables in a prof ob 
    626628      REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: & 
    627629         & zprofmask               ! Mask associated with zprofvar 
     
    631633      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: & 
    632634         & zglam,    &             ! Model longitudes for prof variables 
    633          & zgphi,                  ! Model latitudes for prof variables 
     635         & zgphi                   ! Model latitudes for prof variables 
    634636 
    635637      !----------------------------------------------------------------------- 
     
    666668            SELECT CASE ( TRIM(cobstypesprof(jtype)) ) 
    667669            CASE('prof') 
    668                zprofvar1(:,:,:) = ts(:,:,:,jp_tem,Kmm) 
    669                zprofvar2(:,:,:) = ts(:,:,:,jp_sal,Kmm) 
     670               zprofvar(:,:,:,1) = ts(:,:,:,jp_tem,Kmm) 
     671               zprofvar(:,:,:,2) = ts(:,:,:,jp_sal,Kmm) 
    670672            CASE('vel') 
    671673               zprofvar(:,:,:,1) = uu(:,:,:,Kmm) 
Note: See TracChangeset for help on using the changeset viewer.