Changeset 13767
- Timestamp:
- 2020-11-10T16:42:30+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13747_ENHANCE-04_dford_OBSOP_BGC/src/OCE/OBS/diaobs.F90
r13758 r13767 57 57 PUBLIC calc_date ! Compute the date of a timestep 58 58 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 83 84 INTEGER, DIMENSION(imaxavtypes) :: nn_profdavtypes ! Profile data types representing a daily average 84 85 INTEGER :: nproftypes ! Number of profile obs types … … 133 134 & cn_sicfbfiles, & ! Seaice concentration input filenames 134 135 & cn_velfbfiles, & ! Velocity profile input filenames 135 & cn_sstbiasfiles ! SST bias input filenames136 & cn_sstbiasfiles ! SST bias input filenames 136 137 CHARACTER(LEN=128) :: & 137 138 & cn_altbiasfile ! Altimeter bias input filename … … 139 140 & clproffiles, & ! Profile filenames 140 141 & clsurffiles ! Surface filenames 142 CHARACTER(len=8), DIMENSION(:), ALLOCATABLE :: & 143 & clvars ! Expected variable names 141 144 ! 142 145 LOGICAL :: ln_t3d ! Logical switch for temperature profiles … … 452 455 IF ( TRIM(cobstypesprof(jtype)) == 'prof' ) THEN 453 456 nvarsprof(jtype) = 2 454 nextrprof(jtype) = 1 455 IF ( ln_output_clim ) ltype_clim = .TRUE. 457 nextrprof(jtype) = 1 456 458 ALLOCATE( llvar (nvarsprof(jtype)) ) 457 459 ALLOCATE( clvars(nvarsprof(jtype)) ) … … 556 558 & clsurffiles(jtype,1:ifilessurf(jtype)), & 557 559 & 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), & 559 561 & clvars ) 560 562 ! … … 623 625 INTEGER :: ji, jj ! Loop counters 624 626 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: & 625 & zprofvar ,! Model values for variables in a prof ob627 & zprofvar ! Model values for variables in a prof ob 626 628 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: & 627 629 & zprofmask ! Mask associated with zprofvar … … 631 633 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: & 632 634 & zglam, & ! Model longitudes for prof variables 633 & zgphi ,! Model latitudes for prof variables635 & zgphi ! Model latitudes for prof variables 634 636 635 637 !----------------------------------------------------------------------- … … 666 668 SELECT CASE ( TRIM(cobstypesprof(jtype)) ) 667 669 CASE('prof') 668 zprofvar 1(:,:,:) = ts(:,:,:,jp_tem,Kmm)669 zprofvar 2(:,:,:) = ts(:,:,:,jp_sal,Kmm)670 zprofvar(:,:,:,1) = ts(:,:,:,jp_tem,Kmm) 671 zprofvar(:,:,:,2) = ts(:,:,:,jp_sal,Kmm) 670 672 CASE('vel') 671 673 zprofvar(:,:,:,1) = uu(:,:,:,Kmm)
Note: See TracChangeset
for help on using the changeset viewer.