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 – NEMO

Changeset 15487


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

Remove some comments and add some spaces.

Location:
NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/DYN/dynspg_ts.F90

    r15248 r15487  
    5151   USE agrif_oce 
    5252#endif 
    53 ! REMOVING THIS AS PER V3.6, BUT IS THIS THE BEST BRANCH TO DO IT IN? 
    54 !#if defined key_asminc    
    55 !   USE asminc          ! Assimilation increment 
    56 !#endif 
    5753   ! 
    5854   USE in_out_manager  ! I/O manager 
     
    342338      ENDIF 
    343339      ! 
    344 ! REMOVING THIS AS PER V3.6, BUT IS THIS THE BEST BRANCH TO DO IT IN? 
    345 !#if defined key_asminc 
    346 !      !                                   !=  Add the IAU weighted SSH increment  =! 
    347 !      !                                   !  ------------------------------------  ! 
    348 !      IF( lk_asminc .AND. ln_sshinc .AND. ln_asmiau ) THEN 
    349 !         zssh_frc(:,:) = zssh_frc(:,:) - ssh_iau(:,:) 
    350 !      ENDIF 
    351 !#endif 
    352340      !                                   != Fill boundary data arrays for AGRIF 
    353341      !                                   ! ------------------------------------ 
  • 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 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90

    r15285 r15487  
    210210         zdaystp = 1.0 / REAL( kdaystp ) 
    211211         IF ( idayend == 0 ) THEN 
    212             IF (lwp) WRITE(numout,*) 'Calculating prodatqc%vdmean on time-step: ',kt 
     212            IF (lwp) WRITE(numout,*) 'Calculating prodatqc%vdmean on time-step: ', kt 
    213213            CALL FLUSH(numout) 
    214214            DO jk = 1, jpk 
     
    658658         ! Added kt == 0 test to catch restart case 
    659659         IF ( ( imeanend == 1 ) .OR. ( kt == 0 ) ) THEN 
    660             IF (lwp) WRITE(numout,*) 'Reset surfdataqc%vdmean on time-step: ',kt 
     660            IF (lwp) WRITE(numout,*) 'Reset surfdataqc%vdmean on time-step: ', kt 
    661661            DO jj = 1, jpj 
    662662               DO ji = 1, jpi 
     
    667667 
    668668         ! On each time-step, increment the field for computing time mean 
    669          IF (lwp) WRITE(numout,*)'Accumulating surfdataqc%vdmean on time-step: ',kt 
     669         IF (lwp) WRITE(numout,*)'Accumulating surfdataqc%vdmean on time-step: ', kt 
    670670         DO jj = 1, jpj 
    671671            DO ji = 1, jpi 
     
    678678         IF ( imeanend == 0 ) THEN 
    679679            zmeanstp = 1.0 / REAL( kmeanstp ) 
    680             IF (lwp) WRITE(numout,*)'Calculating surfdataqc%vdmean time mean on time-step: ',kt,' with weight: ',zmeanstp 
     680            IF (lwp) WRITE(numout,*) 'Calculating surfdataqc%vdmean time mean on time-step: ', & 
     681               &                     kt, ' with weight: ', zmeanstp 
    681682            DO jj = 1, jpj 
    682683               DO ji = 1, jpi 
     
    732733         zdaystp = 1.0 / REAL( kdaystp ) 
    733734         IF ( idayend == 0 ) THEN 
    734             IF (lwp) WRITE(numout,*) 'Calculating surfdataqc%vdmean on time-step: ',kt 
     735            IF (lwp) WRITE(numout,*) 'Calculating surfdataqc%vdmean on time-step: ', kt 
    735736            DO jj = 1, jpj 
    736737               DO ji = 1, jpi 
     
    817818         IF ( imeanend == 0 ) THEN 
    818819            ALLOCATE( zsurfm(imaxifp,imaxjfp,isurf) ) 
    819             IF (lwp) WRITE(numout,*)' Interpolating the time mean values on time step: ',kt 
     820            IF (lwp) WRITE(numout,*)' Interpolating the time mean values on time step: ', kt 
    820821            CALL obs_int_comm_2d( imaxifp, imaxjfp, isurf, kpi, kpj, & 
    821822               &                  igrdi, igrdj, surfdataqc%vdmean(:,:,kvar), zsurfm ) 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_prep.F90

    r15248 r15487  
    136136 
    137137      ! Set QC cutoff to optional value if provided 
    138       IF ( PRESENT(kqc_cutoff) ) iqc_cutoff=kqc_cutoff 
     138      IF ( PRESENT(kqc_cutoff) ) iqc_cutoff = kqc_cutoff 
    139139 
    140140      ! ----------------------------------------------------------------------- 
     
    212212         DO jvar = 1, surfdataqc%nvar        
    213213            IF ( jvar == 1 ) THEN 
    214                cout1=TRIM(surfdataqc%cvars(1))                   
     214               cout1 = TRIM(surfdataqc%cvars(1))                   
    215215            ELSE 
    216216               WRITE(cout1,'(A,A1,A)') TRIM(cout1), '/', TRIM(surfdataqc%cvars(jvar))             
     
    245245         WRITE(numout,*) ' Number of observations per time step :' 
    246246         WRITE(numout,*) 
    247          WRITE(numout,'(10X,A,10X,A)')'Time step',TRIM(cout1) 
    248          WRITE(numout,'(10X,A,5X,A)')'---------','-----------------' 
     247         WRITE(numout,'(10X,A,10X,A)') 'Time step', TRIM(cout1) 
     248         WRITE(numout,'(10X,A,5X,A)')  '---------', '-----------------' 
    249249         CALL FLUSH(numout) 
    250250      ENDIF 
     
    377377 
    378378      ! Set QC cutoff to optional value if provided 
    379       IF ( PRESENT(kqc_cutoff) ) iqc_cutoff=kqc_cutoff 
     379      IF ( PRESENT(kqc_cutoff) ) iqc_cutoff = kqc_cutoff 
    380380 
    381381      ! ----------------------------------------------------------------------- 
     
    515515               &            ilanvobsmpp(jvar) 
    516516            IF (ld_nea) THEN 
    517                WRITE(numout,*) ' Remaining '//prodatqc%cvars(jvar)//' data near land points (removed) = ',& 
     517               WRITE(numout,*) ' Remaining '//prodatqc%cvars(jvar)//' data near land points (removed) = ', & 
    518518                  &            inlavobsmpp(jvar) 
    519519            ELSE 
    520                WRITE(numout,*) ' Remaining '//prodatqc%cvars(jvar)//' data near land points (kept)    = ',& 
     520               WRITE(numout,*) ' Remaining '//prodatqc%cvars(jvar)//' data near land points (kept)    = ', & 
    521521                  &            inlavobsmpp(jvar) 
    522522            ENDIF 
     
    530530               ENDIF 
    531531            ENDIF 
    532             WRITE(numout,*) ' Remaining '//prodatqc%cvars(jvar)//' data near open boundary (removed) = ',& 
     532            WRITE(numout,*) ' Remaining '//prodatqc%cvars(jvar)//' data near open boundary (removed) = ', & 
    533533                  &            ibdyvobsmpp(jvar) 
    534534            WRITE(numout,*) ' '//prodatqc%cvars(jvar)//' data accepted                             = ', & 
     
    832832      !! * Local declarations 
    833833      INTEGER :: jobs 
    834       INTEGER :: iqc_cutoff=255 
     834      INTEGER :: iqc_cutoff = 255 
    835835 
    836836      !----------------------------------------------------------------------- 
     
    11221122         & ln_zps,        & 
    11231123         & mbkt 
    1124 ! I THINK MBKT IS CORRECT (V3.6 WAS MBATHY) BUT CONFIRM 
     1124 
    11251125      !! * Arguments 
    11261126      INTEGER, INTENT(IN) :: kprofno      ! Number of profiles 
     
    12741274 
    12751275            ! Calculate max T and W depths of 2x2 grid 
    1276             maxdept=zgdept(1,1,NINT(zbathy(1,1,jobs)),jobs) 
    1277             maxdepw=zgdepw(1,1,NINT(zbathy(1,1,jobs))+1,jobs) 
     1276            maxdept = zgdept(1,1,NINT(zbathy(1,1,jobs)),jobs) 
     1277            maxdepw = zgdepw(1,1,NINT(zbathy(1,1,jobs))+1,jobs) 
    12781278            DO jj = 1, 2 
    12791279               DO ji = 1, 2 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_profiles_def.F90

    r15224 r15487  
    783783         DO ji = 1, prof%nprof 
    784784            IF ( lvalid%luse(ji) ) THEN 
    785                inprof=inprof+1 
     785               inprof = inprof + 1 
    786786               DO jvar = 1, prof%nvar 
    787787                  DO jj = prof%npvsta(ji,jvar), prof%npvend(ji,jvar) 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_read_altbias.F90

    r14075 r15487  
    116116      numaltbias=0 
    117117 
    118       IF(lwp)WRITE(numout,*) 'Opening ',bias_file 
     118      IF (lwp) WRITE(numout,*) 'Opening ', bias_file 
    119119 
    120120      CALL iom_open( bias_file, numaltbias, ldstop=.FALSE. ) 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_read_surf.F90

    r15288 r15487  
    325325            ENDIF 
    326326 
    327             IF (lwp) WRITE(numout,*)'Observation file contains ',inpfiles(jj)%nobs,' observations' 
     327            IF (lwp) WRITE(numout,*) 'Observation file contains ', inpfiles(jj)%nobs, ' observations' 
    328328 
    329329            !------------------------------------------------------------------ 
     
    595595                  ityp(itype+1) = ityp(itype+1) + 1 
    596596               ELSE 
    597                   IF(lwp)WRITE(numout,*)'WARNING:Increase jpsurfmaxtype in ',& 
    598                      &                  cpname 
     597                  IF(lwp)WRITE(numout,*) 'WARNING: Increase jpsurfmaxtype in ', & 
     598                     &                   cpname 
    599599               ENDIF 
    600600 
     
    665665         DO jj = 1,8 
    666666            IF ( itypmpp(jj) > 0 ) THEN 
    667                WRITE(numout,'(1X,A4,I4,A3,I10)')'Type ', jj,' = ',itypmpp(jj) 
     667               WRITE(numout,'(1X,A4,I4,A3,I10)') 'Type ', jj, ' = ', itypmpp(jj) 
    668668            ENDIF 
    669669         END DO 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_surf_def.F90

    r15225 r15487  
    580580      DO ji = 1, surf%nsurf 
    581581 
    582          jj=surf%nsind(ji) 
     582         jj = surf%nsind(ji) 
    583583 
    584584         oldsurf%mi(jj,:)  = surf%mi(ji,:) 
     
    605605         DO ji = 1, surf%nsurf 
    606606             
    607             jj=surf%nsind(ji) 
     607            jj = surf%nsind(ji) 
    608608 
    609609            oldsurf%robs(jj,jk)  = surf%robs(ji,jk) 
     
    622622         DO ji = 1, surf%nsurf 
    623623             
    624             jj=surf%nsind(ji) 
     624            jj = surf%nsind(ji) 
    625625 
    626626            oldsurf%rext(jj,jk)  = surf%rext(ji,jk) 
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_write.F90

    r15288 r15487  
    152152      IF(lwp) THEN 
    153153         WRITE(numout,*) 
    154          WRITE(numout,*)'obs_wri_prof :' 
    155          WRITE(numout,*)'~~~~~~~~~~~~~' 
    156          WRITE(numout,*)'Writing '//TRIM(clfiletype)//' feedback file : ',TRIM(clfname) 
     154         WRITE(numout,*) 'obs_wri_prof :' 
     155         WRITE(numout,*) '~~~~~~~~~~~~~' 
     156         WRITE(numout,*) 'Writing '//TRIM(clfiletype)//' feedback file : ', TRIM(clfname) 
    157157      ENDIF 
    158158 
     
    369369      IF(lwp) THEN 
    370370         WRITE(numout,*) 
    371          WRITE(numout,*)'obs_wri_surf :' 
    372          WRITE(numout,*)'~~~~~~~~~~~~~' 
    373          WRITE(numout,*)'Writing '//TRIM(surfdata%cvars(1))//' feedback file : ',TRIM(clfname) 
     371         WRITE(numout,*) 'obs_wri_surf :' 
     372         WRITE(numout,*) '~~~~~~~~~~~~~' 
     373         WRITE(numout,*) 'Writing '//TRIM(surfdata%cvars(1))//' feedback file : ', TRIM(clfname) 
    374374      ENDIF 
    375375 
     
    494494 
    495495      DO jvar = 1, fbdata%nvar 
    496          zsumx=0.0_wp 
    497          zsumx2=0.0_wp 
    498          inumgoodobs=0 
     496         zsumx = 0.0_wp 
     497         zsumx2 = 0.0_wp 
     498         inumgoodobs = 0 
    499499         DO jo = 1, fbdata%nobs 
    500500            DO jk = 1, fbdata%nlev 
     
    503503                  & ( fbdata%padd(jk,jo,1,jvar) < 9999.0 ) ) THEN 
    504504 
    505                   zomb=fbdata%pob(jk, jo, jvar)-fbdata%padd(jk, jo, 1, jvar) 
    506                   zsumx=zsumx+zomb 
    507                   zsumx2=zsumx2+zomb**2 
    508                   inumgoodobs=inumgoodobs+1 
     505                  zomb = fbdata%pob(jk, jo, jvar)-fbdata%padd(jk, jo, 1, jvar) 
     506                  zsumx = zsumx + zomb 
     507                  zsumx2 = zsumx2 + zomb**2 
     508                  inumgoodobs = inumgoodobs + 1 
    509509               ENDIF 
    510510            ENDDO 
     
    516516 
    517517         IF (lwp) THEN 
    518             WRITE(numout,*) 'Type: ',fbdata%cname(jvar),'  Total number of good observations: ',inumgoodobsmpp 
     518            WRITE(numout,*) 'Type: ', fbdata%cname(jvar), '  Total number of good observations: ', inumgoodobsmpp 
    519519            IF ( inumgoodobsmpp > 0 ) THEN 
    520                WRITE(numout,*) 'Overall mean obs minus model of the good observations: ',zsumx/inumgoodobsmpp 
    521                WRITE(numout,*) 'Overall RMS obs minus model of the good observations: ',sqrt( zsumx2/inumgoodobsmpp ) 
     520               WRITE(numout,*) 'Overall mean obs minus model of the good observations: ', zsumx/inumgoodobsmpp 
     521               WRITE(numout,*) 'Overall RMS obs minus model of the good observations: ', sqrt( zsumx2/inumgoodobsmpp ) 
    522522            ENDIF 
    523523            WRITE(numout,*) '' 
Note: See TracChangeset for help on using the changeset viewer.