Changeset 11456


Ignore:
Timestamp:
2019-08-19T18:27:52+02:00 (14 months ago)
Author:
mattmartin
Message:

First wrking version which produces correct seeming results.

Location:
branches/UKMO/dev_r5518_obs_oper_update_addclim/NEMOGCM/NEMO/OPA_SRC/OBS
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_addclim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r11455 r11456  
    840840            IF ( surfdataqc%lclim ) THEN   
    841841               CALL obs_int_h2d( 1, 1, zweig, zclim(:,:,iobs), zclm ) 
    842                IF (lwp) THEN 
    843                   WRITE(numout,*)'zclim: ', iobs, zclim(:,:,iobs), zclm 
    844                ENDIF 
    845842            ENDIF 
    846843 
  • branches/UKMO/dev_r5518_obs_oper_update_addclim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90

    r11455 r11456  
    698698                                 &                inpfiles(jj)%padd(ij,ji,1,jvar) 
    699699                           ENDIF 
     700                           IF ( profdata%lclim ) THEN 
     701                               profdata%var(jvar)%vclm(ivart(jvar)) = fbrmdi 
     702                           ENDIF                           
    700703                           ! Count number of profile var1 data as function of type 
    701704                           itypvar( profdata%ntyp(iprof) + 1, jvar ) = & 
  • branches/UKMO/dev_r5518_obs_oper_update_addclim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_surf.F90

    r11455 r11456  
    463463                  IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) surfdata%rext(iobs,:) = fbrmdi 
    464464               ENDIF 
    465  
     465                
     466               ! Initialise climatology if set 
     467               IF ( surfdata%lclim ) surfdata%rclm(iobs,1) = fbrmdi 
     468                
    466469               ! STD (obs error standard deviation) read from file and passed through obs operator 
    467470               IF ( iadd_std(jj) /= -1 ) THEN 
Note: See TracChangeset for help on using the changeset viewer.