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 15487 for NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_bias.F90 – NEMO

Ignore:
Timestamp:
2021-11-09T16:27:42+01:00 (3 years ago)
Author:
dford
Message:

Remove some comments and add some spaces.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_bias.F90

    r15187 r15487  
    115115      ENDIF 
    116116       
    117       IF(lwp)WRITE(numout,*)  
    118       IF(lwp)WRITE(numout,*) 'obs_app_bias : ' 
    119       IF(lwp)WRITE(numout,*) '----------------- ' 
     117      IF (lwp) WRITE(numout,*)  
     118      IF (lwp) WRITE(numout,*) 'obs_app_bias : ' 
     119      IF (lwp) WRITE(numout,*) '----------------- ' 
    120120      IF ( ll_extvar ) THEN 
    121          IF(lwp)WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cextvars(kvar)) 
     121         IF (lwp) WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cextvars(kvar)) 
    122122      ELSE 
    123          IF(lwp)WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cvars(kvar)) 
     123         IF (lwp) WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cvars(kvar)) 
    124124      ENDIF 
    125125 
    126126      ! Open and read the files 
    127       z_obsbias(:,:,:)=0.0_wp 
     127      z_obsbias(:,:,:) = 0.0_wp 
    128128      DO jtype = 1, knumtypes 
    129129      
    130          numobsbias=0 
    131          IF(lwp)WRITE(numout,*) 'Opening ',cl_bias_files(jtype) 
     130         numobsbias = 0 
     131         IF (lwp) WRITE(numout,*) 'Opening ', cl_bias_files(jtype) 
    132132         CALL iom_open( cl_bias_files(jtype), numobsbias, ldstop=.FALSE. )        
    133133         IF (numobsbias > 0) THEN 
     
    138138               !so have to use NETCDF 
    139139               !routines directly - should be upgraded in the future 
    140                iret=NF90_OPEN(TRIM(cl_bias_files(jtype)), NF90_NOWRITE, incfile) 
     140               iret = NF90_OPEN(TRIM(cl_bias_files(jtype)), NF90_NOWRITE, incfile) 
    141141               IF ( .NOT. ll_extvar ) THEN 
    142142                  iret=NF90_GET_ATT( incfile, NF90_GLOBAL, TRIM(obsdata%cvars(kvar))//"_source", & 
     
    144144                  ibiastypes(jtype) = ifile_source 
    145145               ENDIF 
    146                iret=NF90_CLOSE(incfile) 
     146               iret = NF90_CLOSE(incfile) 
    147147               IF ( iret /= 0  ) CALL ctl_stop( & 
    148148                  'obs_app_bias : Cannot read bias type from file '// & 
     
    196196            & zmask_tmp(2,2,inumtype), & 
    197197            & zbias( 2,2,inumtype ) ) 
    198          jt=1 
     198         jt = 1 
    199199         DO jobs = 1, obsdata%nsurf  
    200200            IF ( obsdata%ntyp(jobs) == ibiastypes(jtype) ) THEN 
     
    211211               &           igrdi_tmp(:,:,:), igrdj_tmp(:,:,:), & 
    212212               &           z_obsbias(:,:,jtype), zbias(:,:,:) ) 
    213          jt=1 
     213         jt = 1 
    214214         DO jobs = 1, obsdata%nsurf 
    215215            IF ( ( obsdata%ntyp(jobs) == ibiastypes(jtype) ) .OR. & 
     
    230230                  obsdata%robs(jobs,kvar) = obsdata%robs(jobs,kvar) - zext(1) 
    231231               ENDIF 
    232                jt=jt+1 
     232               jt = jt + 1 
    233233            ENDIF 
    234234         END DO  
     
    254254         WRITE(numout,*) "Bias correction applied successfully" 
    255255         IF ( .NOT. ll_extvar ) THEN 
    256             WRITE(numout,*) "Obs types: ",ibiastypes(:), & 
    257                                  " Have all been bias corrected" 
     256            WRITE(numout,*) "Obs types: ", ibiastypes(:), & 
     257               &            " Have all been bias corrected" 
    258258         ENDIF 
    259259      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.