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 9016 for branches/UKMO – NEMO

Changeset 9016 for branches/UKMO


Ignore:
Timestamp:
2017-12-13T15:30:00+01:00 (6 years ago)
Author:
dford
Message:

Change naming convention and fix a couple of bugs - code now compiles and runs.

Location:
branches/UKMO/dev_r5518_obs_oper_update_bgc3d/NEMOGCM
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_obs_oper_update_bgc3d/NEMOGCM/CONFIG/SHARED/namelist_ref

    r8695 r9016  
    12011201   ln_vel3d   = .false.             ! Logical switch for velocity observations 
    12021202   ln_sss     = .false.             ! Logical swithc for SSS observations 
    1203    ln_logchl  = .false.             ! Logical switch for log(Chl) obs 
    1204    ln_spm     = .false.             ! Logical switch for SPM obs 
    1205    ln_fco2    = .false.              
    1206    ln_pco2    = .false. 
     1203   ln_slchltot = .false.            ! Logical switch for surface total              log10(chlorophyll) obs 
     1204   ln_slchldia = .false.            ! Logical switch for surface diatom             log10(chlorophyll) obs 
     1205   ln_slchlnon = .false.            ! Logical switch for surface non-diatom         log10(chlorophyll) obs 
     1206   ln_slchldin = .false.            ! Logical switch for surface dinoflagellate     log10(chlorophyll) obs 
     1207   ln_slchlmic = .false.            ! Logical switch for surface microphytoplankton log10(chlorophyll) obs 
     1208   ln_slchlnan = .false.            ! Logical switch for surface nanophytoplankton  log10(chlorophyll) obs 
     1209   ln_slchlpic = .false.            ! Logical switch for surface picophytoplankton  log10(chlorophyll) obs 
     1210   ln_schltot  = .false.            ! Logical switch for surface total              chlorophyll        obs 
     1211   ln_sspm     = .false.            ! Logical switch for surface suspended particulate matter obs 
     1212   ln_sfco2    = .false.            ! Logical switch for surface fugacity         of carbon dioxide obs 
     1213   ln_spco2    = .false.            ! Logical switch for surface partial pressure of carbon dioxide obs 
    12071214   ln_altbias = .false.             ! Logical switch for altimeter bias correction 
    12081215   ln_sstbias = .false.             ! Logical switch for SST bias correction 

  • branches/UKMO/dev_r5518_obs_oper_update_bgc3d/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r8695 r9016  
    9898      & profdataqc           !: Profile data after quality control 
    9999 
    100    CHARACTER(len=25), PUBLIC, DIMENSION(:), ALLOCATABLE :: & 
     100   CHARACTER(len=8), PUBLIC, DIMENSION(:), ALLOCATABLE :: & 
    101101      & cobstypesprof, &     !: Profile obs types 
    102102      & cobstypessurf        !: Surface obs types 
    … …  
    147147 
    148148      CHARACTER(len=128), DIMENSION(jpmaxnfiles) :: & 
    149          & cn_profbfiles,                 & ! T/S profile input filenames 
    150          & cn_sstfbfiles,                 & ! Sea surface temperature input filenames 
    151          & cn_slafbfiles,                 & ! Sea level anomaly input filenames 
    152          & cn_sicfbfiles,                 & ! Seaice concentration input filenames 
    153          & cn_velfbfiles,                 & ! Velocity profile input filenames 
    154          & cn_sssfbfiles,                 & ! Sea surface salinity input filenames 
    155          & cn_surf_logchl_totalfbfiles,  & ! Surface total              log10(chlorophyll) input filenames 
    156          & cn_surf_logchl_diatfbfiles,    & ! Surface diatom             log10(chlorophyll) input filenames 
    157          & cn_surf_logchl_nondiatfbfiles, & ! Surface non-diatom         log10(chlorophyll) input filenames 
    158          & cn_surf_logchl_dinofbfiles,    & ! Surface dinoflagellate     log10(chlorophyll) input filenames 
    159          & cn_surf_logchl_microfbfiles,  & ! Surface microphytoplankton log10(chlorophyll) input filenames 
    160          & cn_surf_logchl_nanofbfiles,    & ! Surface nanophytoplankton  log10(chlorophyll) input filenames 
    161          & cn_surf_logchl_picofbfiles,    & ! Surface picophytoplankton  log10(chlorophyll) input filenames 
    162          & cn_surf_chl_totalfbfiles,      & ! Surface total              chlorophyll        input filenames 
    163          & cn_surf_spmfbfiles,            & ! Surface suspended particulate matter input filenames 
    164          & cn_surf_fco2fbfiles,           & ! Surface fugacity         of carbon dioxide input filenames 
    165          & cn_surf_pco2fbfiles,           & ! Surface partial pressure of carbon dioxide input filenames 
    166          & cn_sstbiasfiles                  ! SST bias input filenames 
     149         & cn_profbfiles,      & ! T/S profile input filenames 
     150         & cn_sstfbfiles,      & ! Sea surface temperature input filenames 
     151         & cn_slafbfiles,      & ! Sea level anomaly input filenames 
     152         & cn_sicfbfiles,      & ! Seaice concentration input filenames 
     153         & cn_velfbfiles,      & ! Velocity profile input filenames 
     154         & cn_sssfbfiles,      & ! Sea surface salinity input filenames 
     155         & cn_slchltotfbfiles, & ! Surface total              log10(chlorophyll) input filenames 
     156         & cn_slchldiafbfiles, & ! Surface diatom             log10(chlorophyll) input filenames 
     157         & cn_slchlnonfbfiles, & ! Surface non-diatom         log10(chlorophyll) input filenames 
     158         & cn_slchldinfbfiles, & ! Surface dinoflagellate     log10(chlorophyll) input filenames 
     159         & cn_slchlmicfbfiles, & ! Surface microphytoplankton log10(chlorophyll) input filenames 
     160         & cn_slchlnanfbfiles, & ! Surface nanophytoplankton  log10(chlorophyll) input filenames 
     161         & cn_slchlpicfbfiles, & ! Surface picophytoplankton  log10(chlorophyll) input filenames 
     162         & cn_schltotfbfiles,  & ! Surface total              chlorophyll        input filenames 
     163         & cn_sspmfbfiles,     & ! Surface suspended particulate matter input filenames 
     164         & cn_sfco2fbfiles,    & ! Surface fugacity         of carbon dioxide input filenames 
     165         & cn_spco2fbfiles,    & ! Surface partial pressure of carbon dioxide input filenames 
     166         & cn_sstbiasfiles       ! SST bias input filenames 
    167167 
    168168      CHARACTER(LEN=128) :: & 
    … …  
    170170 
    171171 
    172       LOGICAL :: ln_t3d                 ! Logical switch for temperature profiles 
    173       LOGICAL :: ln_s3d                 ! Logical switch for salinity profiles 
    174       LOGICAL :: ln_sla                 ! Logical switch for sea level anomalies  
    175       LOGICAL :: ln_sst                 ! Logical switch for sea surface temperature 
    176       LOGICAL :: ln_sic                 ! Logical switch for sea ice concentration 
    177       LOGICAL :: ln_sss                 ! Logical switch for sea surface salinity obs 
    178       LOGICAL :: ln_vel3d               ! Logical switch for velocity (u,v) obs 
    179       LOGICAL :: ln_surf_logchl_total   ! Logical switch for surface total              log10(chlorophyll) obs 
    180       LOGICAL :: ln_surf_logchl_diat    ! Logical switch for surface diatom             log10(chlorophyll) obs 
    181       LOGICAL :: ln_surf_logchl_nondiat ! Logical switch for surface non-diatom         log10(chlorophyll) obs 
    182       LOGICAL :: ln_surf_logchl_dino    ! Logical switch for surface dinoflagellate     log10(chlorophyll) obs 
    183       LOGICAL :: ln_surf_logchl_micro   ! Logical switch for surface microphytoplankton log10(chlorophyll) obs 
    184       LOGICAL :: ln_surf_logchl_nano    ! Logical switch for surface nanophytoplankton  log10(chlorophyll) obs 
    185       LOGICAL :: ln_surf_logchl_pico    ! Logical switch for surface picophytoplankton  log10(chlorophyll) obs 
    186       LOGICAL :: ln_surf_chl_total      ! Logical switch for surface total              chlorophyll        obs 
    187       LOGICAL :: ln_surf_spm            ! Logical switch for surface suspended particulate matter obs 
    188       LOGICAL :: ln_surf_fco2           ! Logical switch for surface fugacity         of carbon dioxide obs 
    189       LOGICAL :: ln_surf_pco2           ! Logical switch for surface partial pressure of carbon dioxide obs 
    190       LOGICAL :: ln_nea                 ! Logical switch to remove obs near land 
    191       LOGICAL :: ln_altbias             ! Logical switch for altimeter bias 
    192       LOGICAL :: ln_sstbias             ! Logical switch for bias correction of SST 
    193       LOGICAL :: ln_ignmis              ! Logical switch for ignoring missing files 
    194       LOGICAL :: ln_s_at_t              ! Logical switch to compute model S at T obs 
    195       LOGICAL :: ln_bound_reject        ! Logical switch for rejecting obs near the boundary 
     172      LOGICAL :: ln_t3d          ! Logical switch for temperature profiles 
     173      LOGICAL :: ln_s3d          ! Logical switch for salinity profiles 
     174      LOGICAL :: ln_sla          ! Logical switch for sea level anomalies  
     175      LOGICAL :: ln_sst          ! Logical switch for sea surface temperature 
     176      LOGICAL :: ln_sic          ! Logical switch for sea ice concentration 
     177      LOGICAL :: ln_sss          ! Logical switch for sea surface salinity obs 
     178      LOGICAL :: ln_vel3d        ! Logical switch for velocity (u,v) obs 
     179      LOGICAL :: ln_slchltot     ! Logical switch for surface total              log10(chlorophyll) obs 
     180      LOGICAL :: ln_slchldia     ! Logical switch for surface diatom             log10(chlorophyll) obs 
     181      LOGICAL :: ln_slchlnon    ! Logical switch for surface non-diatom         log10(chlorophyll) obs 
     182      LOGICAL :: ln_slchldin     ! Logical switch for surface dinoflagellate     log10(chlorophyll) obs 
     183      LOGICAL :: ln_slchlmic     ! Logical switch for surface microphytoplankton log10(chlorophyll) obs 
     184      LOGICAL :: ln_slchlnan     ! Logical switch for surface nanophytoplankton  log10(chlorophyll) obs 
     185      LOGICAL :: ln_slchlpic     ! Logical switch for surface picophytoplankton  log10(chlorophyll) obs 
     186      LOGICAL :: ln_schltot      ! Logical switch for surface total              chlorophyll        obs 
     187      LOGICAL :: ln_sspm         ! Logical switch for surface suspended particulate matter obs 
     188      LOGICAL :: ln_sfco2        ! Logical switch for surface fugacity         of carbon dioxide obs 
     189      LOGICAL :: ln_spco2        ! Logical switch for surface partial pressure of carbon dioxide obs 
     190      LOGICAL :: ln_nea          ! Logical switch to remove obs near land 
     191      LOGICAL :: ln_altbias      ! Logical switch for altimeter bias 
     192      LOGICAL :: ln_sstbias      ! Logical switch for bias correction of SST 
     193      LOGICAL :: ln_ignmis       ! Logical switch for ignoring missing files 
     194      LOGICAL :: ln_s_at_t       ! Logical switch to compute model S at T obs 
     195      LOGICAL :: ln_bound_reject ! Logical switch for rejecting obs near the boundary 
    196196 
    197197      REAL(dp) :: rn_dobsini     ! Obs window start date YYYYMMDD.HHMMSS 
    … …  
    207207      LOGICAL :: llvar1          ! Logical for profile variable 1 
    208208      LOGICAL :: llvar2          ! Logical for profile variable 1 
    209       LOGICAL :: ltype_fpindegs ! Local version of ln_*_fp_indegs 
     209      LOGICAL :: ltype_fp_indegs ! Local version of ln_*_fp_indegs 
    210210      LOGICAL :: ltype_night     ! Local version of ln_sstnight (false for other variables) 
    211211 
    … …  
    223223      NAMELIST/namobs/ln_diaobs, ln_t3d, ln_s3d, ln_sla,              & 
    224224         &            ln_sst, ln_sic, ln_sss, ln_vel3d,               & 
    225          &            ln_surf_logchl_total,   ln_surf_logchl_diat,    & 
    226          &            ln_surf_logchl_nondiat, ln_surf_logchl_dino,    & 
    227          &            ln_surf_logchl_micro,   ln_surf_logchl_nano,    & 
    228          &            ln_surf_logchl_pico,    ln_surf_chl_total,      & 
    229          &            ln_surf_spm, ln_surf_fco2, ln_surf_pco2,        & 
     225         &            ln_slchltot, ln_slchldia, ln_slchlnon,          & 
     226         &            ln_slchldin, ln_slchlmic, ln_slchlnan,          & 
     227         &            ln_slchlpic, ln_schltot,                        & 
     228         &            ln_sspm,     ln_sfco2,    ln_spco2,             & 
    230229         &            ln_altbias, ln_sstbias, ln_nea,                 & 
    231230         &            ln_grid_global, ln_grid_search_lookup,          & 
    … …  
    237236         &            cn_sstfbfiles, cn_sicfbfiles,                   & 
    238237         &            cn_velfbfiles, cn_sssfbfiles,                   & 
    239          &            cn_surf_logchl_totalfbfiles,                    & 
    240          &            cn_surf_logchl_diatfbfiles,                     & 
    241          &            cn_surf_logchl_nondiatfbfiles,                  & 
    242          &            cn_surf_logchl_dinofbfiles,                     & 
    243          &            cn_surf_logchl_microfbfiles,                    & 
    244          &            cn_surf_logchl_nanofbfiles,                     & 
    245          &            cn_surf_logchl_picofbfiles,                     & 
    246          &            cn_surf_chl_totalfbfiles, cn_surf_spmfbfiles,   & 
    247          &            cn_surf_fco2fbfiles, cn_surf_pco2fbfiles,       & 
     238         &            cn_slchltotfbfiles, cn_slchldiafbfiles,         & 
     239         &            cn_slchlnonfbfiles, cn_slchldinfbfiles,         & 
     240         &            cn_slchlmicfbfiles, cn_slchlnanfbfiles,         & 
     241         &            cn_slchlpicfbfiles,                             & 
     242         &            cn_schltotfbfiles, cn_sspmfbfiles,              & 
     243         &            cn_sfco2fbfiles, cn_spco2fbfiles,               & 
    248244         &            cn_sstbiasfiles, cn_altbiasfile,                & 
    249245         &            cn_gridsearchfile, rn_gridsearchres,            & 
    … …  
    272268 
    273269      ! Some namelist arrays need initialising 
    274       cn_profbfiles(:)                 = '' 
    275       cn_slafbfiles(:)                 = '' 
    276       cn_sstfbfiles(:)                 = '' 
    277       cn_sicfbfiles(:)                 = '' 
    278       cn_velfbfiles(:)                 = '' 
    279       cn_sssfbfiles(:)                 = '' 
    280       cn_surf_logchl_totalfbfiles(:)  = '' 
    281       cn_surf_logchl_diatfbfiles(:)    = '' 
    282       cn_surf_logchl_nondiatfbfiles(:) = '' 
    283       cn_surf_logchl_dinofbfiles(:)    = '' 
    284       cn_surf_logchl_microfbfiles(:)  = '' 
    285       cn_surf_logchl_nanofbfiles(:)    = '' 
    286       cn_surf_logchl_picofbfiles(:)    = '' 
    287       cn_surf_chl_totalfbfiles(:)      = '' 
    288       cn_surf_spmfbfiles(:)            = '' 
    289       cn_surf_fco2fbfiles(:)           = '' 
    290       cn_surf_pco2fbfiles(:)           = '' 
    291       cn_sstbiasfiles(:)               = '' 
    292       nn_profdavtypes(:)               = -1 
     270      cn_profbfiles(:)      = '' 
     271      cn_slafbfiles(:)      = '' 
     272      cn_sstfbfiles(:)      = '' 
     273      cn_sicfbfiles(:)      = '' 
     274      cn_velfbfiles(:)      = '' 
     275      cn_sssfbfiles(:)      = '' 
     276      cn_slchltotfbfiles(:) = '' 
     277      cn_slchldiafbfiles(:) = '' 
     278      cn_slchlnonfbfiles(:) = '' 
     279      cn_slchldinfbfiles(:) = '' 
     280      cn_slchlmicfbfiles(:) = '' 
     281      cn_slchlnanfbfiles(:) = '' 
     282      cn_slchlpicfbfiles(:) = '' 
     283      cn_schltotfbfiles(:)  = '' 
     284      cn_sspmfbfiles(:)     = '' 
     285      cn_sfco2fbfiles(:)    = '' 
     286      cn_spco2fbfiles(:)    = '' 
     287      cn_sstbiasfiles(:)    = '' 
     288      nn_profdavtypes(:)    = -1 
    293289 
    294290      CALL ini_date( rn_dobsini ) 
    … …  
    328324         WRITE(numout,*) '             Logical switch for velocity observations               ln_vel3d = ', ln_vel3d 
    329325         WRITE(numout,*) '             Logical switch for SSS observations                      ln_sss = ', ln_sss 
    330          WRITE(numout,*) '             Logical switch for surf_logchl_total obs   ln_surf_logchl_total = ', ln_surf_logchl_total 
    331          WRITE(numout,*) '             Logical switch for surf_logchl_diat obs     ln_surf_logchl_diat = ', ln_surf_logchl_diat 
    332          WRITE(numout,*) '             Logical switch for surf_logchl_nondiat   ln_surf_logchl_nondiat = ', ln_surf_logchl_nondiat 
    333          WRITE(numout,*) '             Logical switch for surf_logchl_dino obs     ln_surf_logchl_dino = ', ln_surf_logchl_dino 
    334          WRITE(numout,*) '             Logical switch for surf_logchl_micro obs   ln_surf_logchl_micro = ', ln_surf_logchl_micro 
    335          WRITE(numout,*) '             Logical switch for surf_logchl_nano obs     ln_surf_logchl_nano = ', ln_surf_logchl_nano 
    336          WRITE(numout,*) '             Logical switch for surf_logchl_pico obs     ln_surf_logchl_pico = ', ln_surf_logchl_pico 
    337          WRITE(numout,*) '             Logical switch for surf_chl_total obs         ln_surf_chl_total = ', ln_surf_chl_total 
    338          WRITE(numout,*) '             Logical switch for surface SPM observations         ln_surf_spm = ', ln_surf_spm 
    339          WRITE(numout,*) '             Logical switch for surface fCO2 observations       ln_surf_fco2 = ', ln_surf_fco2 
    340          WRITE(numout,*) '             Logical switch for surface pCO2 observations       ln_surf_pco2 = ', ln_surf_pco2 
     326         WRITE(numout,*) '             Logical switch for surface total logchl obs         ln_slchltot = ', ln_slchltot 
     327         WRITE(numout,*) '             Logical switch for surface diatom logchl obs        ln_slchldia = ', ln_slchldia 
     328         WRITE(numout,*) '             Logical switch for surface non-diatom logchl obs    ln_slchlnon = ', ln_slchlnon 
     329         WRITE(numout,*) '             Logical switch for surface dino logchl obs          ln_slchldin = ', ln_slchldin 
     330         WRITE(numout,*) '             Logical switch for surface micro logchl obs         ln_slchlmic = ', ln_slchlmic 
     331         WRITE(numout,*) '             Logical switch for surface nano logchl obs          ln_slchlnan = ', ln_slchlnan 
     332         WRITE(numout,*) '             Logical switch for surface pico logchl obs          ln_slchlpic = ', ln_slchlpic 
     333         WRITE(numout,*) '             Logical switch for surface total chl obs             ln_schltot = ', ln_schltot 
     334         WRITE(numout,*) '             Logical switch for surface SPM observations             ln_sspm = ', ln_sspm 
     335         WRITE(numout,*) '             Logical switch for surface fCO2 observations           ln_sfco2 = ', ln_sfco2 
     336         WRITE(numout,*) '             Logical switch for surface pCO2 observations           ln_spco2 = ', ln_spco2 
    341337         WRITE(numout,*) '             Global distribution of observations              ln_grid_global = ', ln_grid_global 
    342338         WRITE(numout,*) '             Logical switch for obs grid search lookup ln_grid_search_lookup = ', ln_grid_search_lookup 
    … …  
    380376 
    381377      nproftypes = COUNT( (/ln_t3d .OR. ln_s3d, ln_vel3d /) ) 
    382       nsurftypes = COUNT( (/ln_sla, ln_sst, ln_sic, ln_sss,              & 
    383          &                  ln_surf_logchl_total, ln_surf_logchl_diat,   & 
    384          &                  ln_surf_logchl_nondiat, ln_surf_logchl_dino, & 
    385          &                  ln_surf_logchl_micro, ln_surf_logchl_nano,   & 
    386          &                  ln_surf_logchl_pico, ln_surf_chl_total,      & 
    387          &                  ln_surf_spm, ln_surf_fco2, ln_surf_pco2 /) ) 
     378      nsurftypes = COUNT( (/ln_sla, ln_sst, ln_sic, ln_sss,                     & 
     379         &                  ln_slchltot, ln_slchldia, ln_slchlnon, ln_slchldin, & 
     380         &                  ln_slchlmic, ln_slchlnan, ln_slchlpic, ln_schltot,  & 
     381         &                  ln_sspm, ln_sfco2, ln_spco2 /) ) 
    388382 
    389383      IF ( nproftypes == 0 .AND. nsurftypes == 0 ) THEN 
    … …  
    453447         ENDIF 
    454448 
    455          IF (ln_surf_logchl_total) THEN 
     449         IF (ln_slchltot) THEN 
    456450            jtype = jtype + 1 
    457             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_total', & 
    458                &                   cn_surf_logchl_totalfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    459          ENDIF 
    460  
    461          IF (ln_surf_logchl_diat) THEN 
     451            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchltot', & 
     452               &                   cn_slchltotfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     453         ENDIF 
     454 
     455         IF (ln_slchldia) THEN 
    462456            jtype = jtype + 1 
    463             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_diat', & 
    464                &                   cn_surf_logchl_diatfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    465          ENDIF 
    466  
    467          IF (ln_surf_logchl_nondiat) THEN 
     457            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchldia', & 
     458               &                   cn_slchldiafbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     459         ENDIF 
     460 
     461         IF (ln_slchlnon) THEN 
    468462            jtype = jtype + 1 
    469             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_nondiat', & 
    470                &                   cn_surf_logchl_nondiatfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    471          ENDIF 
    472  
    473          IF (ln_surf_logchl_dino) THEN 
     463            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchlnon', & 
     464               &                   cn_slchlnonfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     465         ENDIF 
     466 
     467         IF (ln_slchldin) THEN 
    474468            jtype = jtype + 1 
    475             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_dino', & 
    476                &                   cn_surf_logchl_dinofbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    477          ENDIF 
    478  
    479          IF (ln_surf_logchl_micro) THEN 
     469            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchldin', & 
     470               &                   cn_slchldinfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     471         ENDIF 
     472 
     473         IF (ln_slchlmic) THEN 
    480474            jtype = jtype + 1 
    481             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_micro', & 
    482                &                   cn_surf_logchl_microfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    483          ENDIF 
    484  
    485          IF (ln_surf_logchl_nano) THEN 
     475            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchlmic', & 
     476               &                   cn_slchlmicfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     477         ENDIF 
     478 
     479         IF (ln_slchlnan) THEN 
    486480            jtype = jtype + 1 
    487             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_nano', & 
    488                &                   cn_surf_logchl_nanofbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    489          ENDIF 
    490  
    491          IF (ln_surf_logchl_pico) THEN 
     481            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchlnan', & 
     482               &                   cn_slchlnanfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     483         ENDIF 
     484 
     485         IF (ln_slchlpic) THEN 
    492486            jtype = jtype + 1 
    493             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_logchl_pico', & 
    494                &                   cn_surf_logchl_picofbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    495          ENDIF 
    496  
    497          IF (ln_surf_chl_total) THEN 
     487            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'slchlpic', & 
     488               &                   cn_slchlpicfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     489         ENDIF 
     490 
     491         IF (ln_schltot) THEN 
    498492            jtype = jtype + 1 
    499             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_chl_total', & 
    500                &                   cn_surf_chl_totalfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    501          ENDIF 
    502  
    503          IF (ln_surf_spm) THEN 
     493            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'schltot', & 
     494               &                   cn_schltotfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     495         ENDIF 
     496 
     497         IF (ln_sspm) THEN 
    504498            jtype = jtype + 1 
    505             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_spm', & 
    506                &                   cn_surf_spmfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    507          ENDIF 
    508  
    509          IF (ln_surf_fco2) THEN 
     499            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'sspm', & 
     500               &                   cn_sspmfbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     501         ENDIF 
     502 
     503         IF (ln_sfco2) THEN 
    510504            jtype = jtype + 1 
    511             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_fco2', & 
    512                &                   cn_surf_fco2fbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    513          ENDIF 
    514  
    515          IF (ln_surf_pco2) THEN 
     505            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'sfco2', & 
     506               &                   cn_sfco2fbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     507         ENDIF 
     508 
     509         IF (ln_spco2) THEN 
    516510            jtype = jtype + 1 
    517             CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'surf_pco2', & 
    518                &                   cn_surf_pco2fbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
     511            CALL obs_settypefiles( nsurftypes, jpmaxnfiles, jtype, 'spco2', & 
     512               &                   cn_spco2fbfiles, ifilessurf, cobstypessurf, clsurffiles ) 
    519513         ENDIF 
    520514 
    521515         DO jtype = 1, nsurftypes 
    522516 
    523             IF ( TRIM(n2dintsurf(jtype)) == 'sla' ) THEN 
     517            IF ( TRIM(cobstypessurf(jtype)) == 'sla' ) THEN 
    524518               IF ( nn_2dint_sla == -1 ) THEN 
    525519                  n2dint_type  = nn_2dint_default 
    … …  
    529523               rtype_avglamscl = rn_sla_avglamscl 
    530524               rtype_avgphiscl = rn_sla_avgphiscl 
    531                ltype_fpindegs  = ln_sla_fpindegs 
     525               ltype_fp_indegs = ln_sla_fp_indegs 
    532526               ltype_night     = .FALSE. 
    533             ELSE IF ( TRIM(n2dintsurf(jtype)) == 'sst' ) THEN 
     527            ELSE IF ( TRIM(cobstypessurf(jtype)) == 'sst' ) THEN 
    534528               IF ( nn_2dint_sst == -1 ) THEN 
    535529                  n2dint_type  = nn_2dint_default 
    … …  
    539533               rtype_avglamscl = rn_sst_avglamscl 
    540534               rtype_avgphiscl = rn_sst_avgphiscl 
    541                ltype_fpindegs  = ln_sst_fpindegs 
     535               ltype_fp_indegs = ln_sst_fp_indegs 
    542536               ltype_night     = ln_sstnight 
    543             ELSE IF ( TRIM(n2dintsurf(jtype)) == 'sic' ) THEN 
     537            ELSE IF ( TRIM(cobstypessurf(jtype)) == 'sic' ) THEN 
    544538               IF ( nn_2dint_sic == -1 ) THEN 
    545539                  n2dint_type  = nn_2dint_default 
    … …  
    549543               rtype_avglamscl = rn_sic_avglamscl 
    550544               rtype_avgphiscl = rn_sic_avgphiscl 
    551                ltype_fpindegs  = ln_sic_fpindegs 
     545               ltype_fp_indegs = ln_sic_fp_indegs 
    552546               ltype_night     = .FALSE. 
    553             ELSE IF ( TRIM(n2dintsurf(jtype)) == 'sss' ) THEN 
     547            ELSE IF ( TRIM(cobstypessurf(jtype)) == 'sss' ) THEN 
    554548               IF ( nn_2dint_sss == -1 ) THEN 
    555549                  n2dint_type  = nn_2dint_default 
    … …  
    559553               rtype_avglamscl = rn_sss_avglamscl 
    560554               rtype_avgphiscl = rn_sss_avgphiscl 
    561                ltype_fpindegs  = ln_sss_fpindegs 
     555               ltype_fp_indegs = ln_sss_fp_indegs 
    562556               ltype_night     = .FALSE. 
    563557            ELSE 
    … …  
    565559               rtype_avglamscl = rn_default_avglamscl 
    566560               rtype_avgphiscl = rn_default_avgphiscl 
    567                ltype_fpindegs  = ln_default_fpindegs 
     561               ltype_fp_indegs = ln_default_fp_indegs 
    568562               ltype_night     = .FALSE. 
    569563            ENDIF 
    570564             
    571             CALL obs_setinterpopts( nsurftypes, jtype, TRIM(n2dintsurf(jtype)), & 
    572                &                    nn_2dint_default, n2dint_type,              & 
    573                &                    rtype_avglamscl, rtype_avgphiscl,           & 
    574                &                    ltype_fpindegs, ltype_night,                & 
    575                &                    n2dintsurf, ravglamscl, ravgphiscl,         & 
     565            CALL obs_setinterpopts( nsurftypes, jtype, TRIM(cobstypessurf(jtype)), & 
     566               &                    nn_2dint_default, n2dint_type,                 & 
     567               &                    rtype_avglamscl, rtype_avgphiscl,              & 
     568               &                    ltype_fp_indegs, ltype_night,                  & 
     569               &                    n2dintsurf, ravglamscl, ravgphiscl,            & 
    576570               &                    lfpindegs, llnightav ) 
    577571 
    … …  
    923917               ENDIF 
    924918 
    925             CASE('surf_logchl_total') 
     919            CASE('slchltot') 
    926920#if defined key_hadocc 
    927921               ! Surface chlorophyll from HadOCC 
    … …  
    935929                  &            trn(:,:,1,jp_fabm_chl3) + trn(:,:,1,jp_fabm_chl4) 
    936930#else 
    937                CALL ctl_stop( ' Trying to run surf_logchl_total observation operator', & 
     931               CALL ctl_stop( ' Trying to run slchltot observation operator', & 
    938932                  &           ' but no biogeochemical model appears to have been defined' ) 
    939933#endif 
    … …  
    948942               END WHERE 
    949943 
    950             CASE('surf_logchl_diat') 
     944            CASE('slchldia') 
    951945#if defined key_hadocc 
    952                CALL ctl_stop( ' Trying to run surf_logchl_diat observation operator', & 
     946               CALL ctl_stop( ' Trying to run slchldia observation operator', & 
    953947                  &           ' but HadOCC does not explicitly simulate diatoms' ) 
    954948#elif defined key_medusa && defined key_foam_medusa 
    … …  
    959953               zsurfvar(:,:) = trn(:,:,1,jp_fabm_chl1) 
    960954#else 
    961                CALL ctl_stop( ' Trying to run surf_logchl_diat observation operator', & 
     955               CALL ctl_stop( ' Trying to run slchldia observation operator', & 
    962956                  &           ' but no biogeochemical model appears to have been defined' ) 
    963957#endif 
    … …  
    972966               END WHERE 
    973967 
    974             CASE('surf_logchl_nondiat') 
     968            CASE('slchlnon') 
    975969#if defined key_hadocc 
    976                CALL ctl_stop( ' Trying to run surf_logchl_nondiat observation operator', & 
     970               CALL ctl_stop( ' Trying to run slchlnon observation operator', & 
    977971                  &           ' but HadOCC does not explicitly simulate non-diatoms' ) 
    978972#elif defined key_medusa && defined key_foam_medusa 
    … …  
    984978                  &            trn(:,:,1,jp_fabm_chl3) + trn(:,:,1,jp_fabm_chl4) 
    985979#else 
    986                CALL ctl_stop( ' Trying to run surf_logchl_nondiat observation operator', & 
     980               CALL ctl_stop( ' Trying to run slchlnon observation operator', & 
    987981                  &           ' but no biogeochemical model appears to have been defined' ) 
    988982#endif 
    … …  
    997991               END WHERE 
    998992 
    999             CASE('surf_logchl_dino') 
     993            CASE('slchldin') 
    1000994#if defined key_hadocc 
    1001                CALL ctl_stop( ' Trying to run surf_logchl_dino observation operator', & 
     995               CALL ctl_stop( ' Trying to run slchldin observation operator', & 
    1002996                  &           ' but HadOCC does not explicitly simulate dinoflagellates' ) 
    1003997#elif defined key_medusa && defined key_foam_medusa 
    1004                CALL ctl_stop( ' Trying to run surf_logchl_dino observation operator', & 
     998               CALL ctl_stop( ' Trying to run slchldin observation operator', & 
    1005999                  &           ' but MEDUSA does not explicitly simulate dinoflagellates' ) 
    10061000#elif defined key_fabm 
    … …  
    10081002               zsurfvar(:,:) = trn(:,:,1,jp_fabm_chl4) 
    10091003#else 
    1010                CALL ctl_stop( ' Trying to run surf_logchl_dino observation operator', & 
     1004               CALL ctl_stop( ' Trying to run slchldin observation operator', & 
    10111005                  &           ' but no biogeochemical model appears to have been defined' ) 
    10121006#endif 
    … …  
    10211015               END WHERE 
    10221016 
    1023             CASE('surf_logchl_micro') 
     1017            CASE('slchlmic') 
    10241018#if defined key_hadocc 
    1025                CALL ctl_stop( ' Trying to run surf_logchl_micro observation operator', & 
     1019               CALL ctl_stop( ' Trying to run slchlmic observation operator', & 
    10261020                  &           ' but HadOCC does not explicitly simulate microphytoplankton' ) 
    10271021#elif defined key_medusa && defined key_foam_medusa 
    1028                CALL ctl_stop( ' Trying to run surf_logchl_micro observation operator', & 
     1022               CALL ctl_stop( ' Trying to run slchlmic observation operator', & 
    10291023                  &           ' but MEDUSA does not explicitly simulate microphytoplankton' ) 
    10301024#elif defined key_fabm 
    … …  
    10321026               zsurfvar(:,:) = trn(:,:,1,jp_fabm_chl1) + trn(:,:,1,jp_fabm_chl4) 
    10331027#else 
    1034                CALL ctl_stop( ' Trying to run surf_logchl_micro observation operator', & 
     1028               CALL ctl_stop( ' Trying to run slchlmic observation operator', & 
    10351029                  &           ' but no biogeochemical model appears to have been defined' ) 
    10361030#endif 
    … …  
    10451039               END WHERE 
    10461040 
    1047             CASE('surf_logchl_nano') 
     1041            CASE('slchlnan') 
    10481042#if defined key_hadocc 
    1049                CALL ctl_stop( ' Trying to run surf_logchl_nano observation operator', & 
     1043               CALL ctl_stop( ' Trying to run slchlnan observation operator', & 
    10501044                  &           ' but HadOCC does not explicitly simulate nanophytoplankton' ) 
    10511045#elif defined key_medusa && defined key_foam_medusa 
    1052                CALL ctl_stop( ' Trying to run surf_logchl_nano observation operator', & 
     1046               CALL ctl_stop( ' Trying to run slchlnan observation operator', & 
    10531047                  &           ' but MEDUSA does not explicitly simulate nanophytoplankton' ) 
    10541048#elif defined key_fabm 
    … …  
    10561050               zsurfvar(:,:) = trn(:,:,1,jp_fabm_chl2) 
    10571051#else 
    1058                CALL ctl_stop( ' Trying to run surf_logchl_nano observation operator', & 
     1052               CALL ctl_stop( ' Trying to run slchlnan observation operator', & 
    10591053                  &           ' but no biogeochemical model appears to have been defined' ) 
    10601054#endif 
    … …  
    10691063               END WHERE 
    10701064 
    1071             CASE('surf_logchl_pico') 
     1065            CASE('slchlpic') 
    10721066#if defined key_hadocc 
    1073                CALL ctl_stop( ' Trying to run surf_logchl_pico observation operator', & 
     1067               CALL ctl_stop( ' Trying to run slchlpic observation operator', & 
    10741068                  &           ' but HadOCC does not explicitly simulate picophytoplankton' ) 
    10751069#elif defined key_medusa && defined key_foam_medusa 
    1076                CALL ctl_stop( ' Trying to run surf_logchl_pico observation operator', & 
     1070               CALL ctl_stop( ' Trying to run slchlpic observation operator', & 
    10771071                  &           ' but MEDUSA does not explicitly simulate picophytoplankton' ) 
    10781072#elif defined key_fabm 
    … …  
    10801074               zsurfvar(:,:) = trn(:,:,1,jp_fabm_chl3) 
    10811075#else 
    1082                CALL ctl_stop( ' Trying to run surf_logchl_pico observation operator', & 
     1076               CALL ctl_stop( ' Trying to run slchlpic observation operator', & 
    10831077                  &           ' but no biogeochemical model appears to have been defined' ) 
    10841078#endif 
    … …  
    10931087               END WHERE 
    10941088 
    1095             CASE('surf_chl_total') 
     1089            CASE('schltot') 
    10961090#if defined key_hadocc 
    10971091               ! Surface chlorophyll from HadOCC 
    … …  
    11051099                  &            trn(:,:,1,jp_fabm_chl3) + trn(:,:,1,jp_fabm_chl4) 
    11061100#else 
    1107                CALL ctl_stop( ' Trying to run surf_chl_total observation operator', & 
     1101               CALL ctl_stop( ' Trying to run schltot observation operator', & 
    11081102                  &           ' but no biogeochemical model appears to have been defined' ) 
    11091103#endif 
    11101104               zsurfmask(:,:) = tmask(:,:,1)         ! create a special mask to exclude certain things 
    11111105 
    1112             CASE('surf_spm') 
     1106            CASE('sspm') 
    11131107#if defined key_spm 
    11141108               zsurfvar(:,:) = 0.0 
    … …  
    11171111               END DO 
    11181112#else 
    1119                CALL ctl_stop( ' Trying to run spm observation operator', & 
     1113               CALL ctl_stop( ' Trying to run sspm observation operator', & 
    11201114                  &           ' but no spm model appears to have been defined' ) 
    11211115#endif 
    11221116 
    1123             CASE('surf_fco2') 
     1117            CASE('sfco2') 
    11241118#if defined key_hadocc 
    11251119               zsurfvar(:,:) = HADOCC_FCO2(:,:)    ! fCO2 from HadOCC 
    … …  
    11571151                  &            (82.0578 * (tsn(:,:,1,jp_tem)+rt0))) 
    11581152#else 
    1159                CALL ctl_stop( ' Trying to run fco2 observation operator', & 
     1153               CALL ctl_stop( ' Trying to run sfco2 observation operator', & 
    11601154                  &           ' but no biogeochemical model appears to have been defined' ) 
    11611155#endif 
    11621156 
    1163             CASE('surf_pco2') 
     1157            CASE('spco2') 
    11641158#if defined key_hadocc 
    11651159               zsurfvar(:,:) = HADOCC_PCO2(:,:)    ! pCO2 from HadOCC 
    … …  
    11771171               zsurfvar(:,:) = pco2_3d(:,:,1) 
    11781172#else 
    1179                CALL ctl_stop( ' Trying to run pCO2 observation operator', & 
     1173               CALL ctl_stop( ' Trying to run spco2 observation operator', & 
    11801174                  &           ' but no biogeochemical model appears to have been defined' ) 
    11811175#endif 
    … …  
    14991493          &                   ifiles      ! Out appended number of files for this type 
    15001494 
    1501        CHARACTER(len=25), INTENT(IN) :: ctypein  
     1495       CHARACTER(len=8), INTENT(IN) :: ctypein  
    15021496       CHARACTER(len=128), DIMENSION(jpmaxnfiles), INTENT(IN) :: & 
    15031497          &                   cfilestype  ! In list of files for this obs type 
    1504        CHARACTER(len=25), DIMENSION(ntypes), INTENT(INOUT) :: & 
     1498       CHARACTER(len=8), DIMENSION(ntypes), INTENT(INOUT) :: & 
    15051499          &                   cobstypes   ! Out appended list of obs types 
    15061500       CHARACTER(len=128), DIMENSION(ntypes, jpmaxnfiles), INTENT(INOUT) :: & 
    … …  
    15481542       LOGICAL, INTENT(IN)  :: lfp_indegs_type    !T=> footprint in degrees, F=> in metres 
    15491543       LOGICAL, INTENT(IN)  :: lavnight_type      !T=> obs represent night time average 
    1550        CHARACTER(len=25), INTENT(IN) :: ctypein  
     1544       CHARACTER(len=8), INTENT(IN) :: ctypein  
    15511545 
    15521546       INTEGER, DIMENSION(ntypes), INTENT(INOUT) :: & 
    … …  
    15591553       lavnight(jtype) = lavnight_type 
    15601554 
    1561        IF ( (n2dint_type >= 1) .AND. (n2dint_type <= 6) ) THEN 
     1555       IF ( (n2dint_type >= 0) .AND. (n2dint_type <= 6) ) THEN 
    15621556          n2dint(jtype) = n2dint_type 
    15631557       ELSE IF ( n2dint_type == -1 ) THEN 

  • branches/UKMO/dev_r5518_obs_oper_update_bgc3d/NEMOGCM/NEMO/OPA_SRC/OBS/obs_write.F90

    r8693 r9016  
    8383      TYPE(obfbdata) :: fbdata 
    8484      CHARACTER(LEN=40) :: clfname 
    85       CHARACTER(LEN=27) :: clfiletype 
     85      CHARACTER(LEN=10) :: clfiletype 
    8686      INTEGER :: ilevel 
    8787      INTEGER :: jvar 
    … …  
    320320      TYPE(obfbdata) :: fbdata 
    321321      CHARACTER(LEN=40) :: clfname         ! netCDF filename 
    322       CHARACTER(LEN=27) :: clfiletype 
     322      CHARACTER(LEN=10) :: clfiletype 
    323323      CHARACTER(LEN=12), PARAMETER :: cpname = 'obs_wri_surf' 
    324324      CHARACTER(LEN=ilenlong) :: cllongname  ! Long name of variable 
    … …  
    399399         clgrid     = 'T' 
    400400 
    401       CASE('SURF_LOGCHL_TOTAL','LOGCHL','LogChl','logchl') 
    402  
    403          clfiletype = 'surf_logchl_totalfb' 
     401      CASE('SLCHLTOT','LOGCHL','LogChl','logchl') 
     402 
     403         clfiletype = 'slchltotfb' 
    404404         cllongname = 'Surface total log10(chlorophyll)' 
    405405         clunits    = 'log10(mg/m3)' 
    406406         clgrid     = 'T' 
    407407 
    408       CASE('SURF_LOGCHL_DIAT') 
    409  
    410          clfiletype = 'surf_logchl_diatfb' 
     408      CASE('SLCHLDIA') 
     409 
     410         clfiletype = 'slchldiafb' 
    411411         cllongname = 'Surface diatom log10(chlorophyll)' 
    412412         clunits    = 'log10(mg/m3)' 
    413413         clgrid     = 'T' 
    414414 
    415       CASE('SURF_LOGCHL_NONDIAT') 
    416  
    417          clfiletype = 'surf_logchl_nondiatfb' 
     415      CASE('SLCHLNON') 
     416 
     417         clfiletype = 'slchlnonfb' 
    418418         cllongname = 'Surface non-diatom log10(chlorophyll)' 
    419419         clunits    = 'log10(mg/m3)' 
    420420         clgrid     = 'T' 
    421421 
    422       CASE('SURF_LOGCHL_DINO') 
    423  
    424          clfiletype = 'surf_logchl_dinofb' 
     422      CASE('SLCHLDIN') 
     423 
     424         clfiletype = 'slchldinfb' 
    425425         cllongname = 'Surface dinoflagellate log10(chlorophyll)' 
    426426         clunits    = 'log10(mg/m3)' 
    427427         clgrid     = 'T' 
    428428 
    429       CASE('SURF_LOGCHL_MICRO') 
    430  
    431          clfiletype = 'surf_logchl_microfb' 
     429      CASE('SLCHLMIC') 
     430 
     431         clfiletype = 'slchlmicfb' 
    432432         cllongname = 'Surface microphytoplankton log10(chlorophyll)' 
    433433         clunits    = 'log10(mg/m3)' 
    434434         clgrid     = 'T' 
    435435 
    436       CASE('SURF_LOGCHL_NANO') 
    437  
    438          clfiletype = 'surf_logchl_nanofb' 
     436      CASE('SLCHLNAN') 
     437 
     438         clfiletype = 'slchlnanfb' 
    439439         cllongname = 'Surface nanophytoplankton log10(chlorophyll)' 
    440440         clunits    = 'log10(mg/m3)' 
    441441         clgrid     = 'T' 
    442442 
    443       CASE('SURF_LOGCHL_PICO') 
    444  
    445          clfiletype = 'surf_logchl_picofb' 
     443      CASE('SLCHLPIC') 
     444 
     445         clfiletype = 'slchlpicfb' 
    446446         cllongname = 'Surface picophytoplankton log10(chlorophyll)' 
    447447         clunits    = 'log10(mg/m3)' 
    448448         clgrid     = 'T' 
    449449 
    450       CASE('SURF_CHL_TOTAL') 
    451  
    452          clfiletype = 'surf_chl_totalfb' 
     450      CASE('SCHLTOT') 
     451 
     452         clfiletype = 'schltotfb' 
    453453         cllongname = 'Surface total chlorophyll' 
    454454         clunits    = 'mg/m3' 
    455455         clgrid     = 'T' 
    456456 
    457       CASE('SURF_SPM') 
    458  
    459          clfiletype = 'surf_spmfb' 
     457      CASE('SSPM') 
     458 
     459         clfiletype = 'sspmfb' 
    460460         cllongname = 'Surface suspended particulate matter' 
    461461         clunits    = 'g/m3' 
    462462         clgrid     = 'T' 
    463463 
    464       CASE('SURF_FCO2','FCO2','fCO2','fco2') 
    465  
    466          clfiletype = 'surf_fco2fb' 
     464      CASE('SFCO2','FCO2','fCO2','fco2') 
     465 
     466         clfiletype = 'sfco2fb' 
    467467         cllongname = 'Surface fugacity of carbon dioxide' 
    468468         clunits    = 'uatm' 
    469469         clgrid     = 'T' 
    470470 
    471       CASE('SURF_PCO2') 
    472  
    473          clfiletype = 'surf_pco2fb' 
     471      CASE('SPCO2') 
     472 
     473         clfiletype = 'spco2fb' 
    474474         cllongname = 'Surface partial pressure of carbon dioxide' 
    475475         clunits    = 'uatm' 
Note: See TracChangeset for help on using the changeset viewer.