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.
Ticket Diff – NEMO

Changes between Initial Version and Version 1 of Ticket #2564


Ignore:
Timestamp:
2020-12-03T13:49:00+01:00 (3 years ago)
Author:
hadcv
Comment:

Legend:

Unmodified
Added
Removed
Modified

  • Ticket #2564

    • Property Owner changed from systeam to hadcv
    • Property Status changed from new to assigned

  • Ticket #2564 – Description

    initial v1  
    5050 
    5151I gather that this is because `sf_tsd(jp_tem)%fnow` is only updated once per record when temporal interpolation is not used. 
     52 
     53**Solution** 
     54 
     55Check to see if the record number has changed (implying new data has been read in): 
     56 
     57 {{{#!diff 
     58         IF( cn_cfg == "orca" .OR. cn_cfg == "ORCA" ) THEN 
     59            IF( nn_cfg == 2 .AND. ln_tsd_dmp ) THEN    ! some hand made alterations 
     60+              irec_n(jp_tem) = sf_tsd(jp_tem)%nrec(1,sf_tsd(jp_tem)%naa)            ! Determine if there is new data (ln_tint = F) 
     61+              irec_n(jp_sal) = sf_tsd(jp_sal)%nrec(1,sf_tsd(jp_sal)%naa)            ! If not, then do not apply the increments 
     62+              IF( kt == nit000 ) irec_b = irec_n 
     63               ! 
     64               ij0 = 101 + nn_hls       ;   ij1 = 109 + nn_hls                       ! Reduced T & S in the Alboran Sea 
     65               ii0 = 141 + nn_hls - 1   ;   ii1 = 155 + nn_hls - 1 
     66               DO jj = mj0(ij0), mj1(ij1) 
     67                  DO ji = mi0(ii0), mi1(ii1) 
     68+                    IF( sf_tsd(jp_tem)%ln_tint .OR. irec_n(jp_tem) /= irec_b(jp_tem) ) THEN 
     69                        sf_tsd(jp_tem)%fnow(ji,jj,13:13) = sf_tsd(jp_tem)%fnow(ji,jj,13:13) - 0.20_wp 
     70                        sf_tsd(jp_tem)%fnow(ji,jj,14:15) = sf_tsd(jp_tem)%fnow(ji,jj,14:15) - 0.35_wp 
     71                        sf_tsd(jp_tem)%fnow(ji,jj,16:25) = sf_tsd(jp_tem)%fnow(ji,jj,16:25) - 0.40_wp 
     72+                       irec_b(jp_tem) = irec_n(jp_tem) 
     73+                    ENDIF 
     74                     ! 
     75+                    IF( sf_tsd(jp_sal)%ln_tint .OR. irec_n(jp_sal) /= irec_b(jp_sal) ) THEN 
     76                        sf_tsd(jp_sal)%fnow(ji,jj,13:13) = sf_tsd(jp_sal)%fnow(ji,jj,13:13) - 0.15_wp 
     77                        sf_tsd(jp_sal)%fnow(ji,jj,14:15) = sf_tsd(jp_sal)%fnow(ji,jj,14:15) - 0.25_wp 
     78                        sf_tsd(jp_sal)%fnow(ji,jj,16:17) = sf_tsd(jp_sal)%fnow(ji,jj,16:17) - 0.30_wp 
     79                        sf_tsd(jp_sal)%fnow(ji,jj,18:25) = sf_tsd(jp_sal)%fnow(ji,jj,18:25) - 0.35_wp 
     80+                       irec_b(jp_sal) = irec_n(jp_sal) 
     81+                    ENDIF 
     82                  END DO 
     83               END DO 
     84               ij0 =  87 + nn_hls       ;   ij1 =  96 + nn_hls                       ! Reduced temperature in Red Sea 
     85               ii0 = 148 + nn_hls - 1   ;   ii1 = 160 + nn_hls - 1 
     86               sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ,  4:10 ) = 7.0_wp 
     87               sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 11:13 ) = 6.5_wp 
     88               sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 14:20 ) = 6.0_wp 
     89            ENDIF 
     90         ENDIF 
     91 }}}