Ignore:
Timestamp:
2019-06-11T14:58:51+02:00 (18 months ago)
Author:
agn
Message:

complete merge of diahth.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/src/OCE/DIA/diahth.F90

    r11093 r11094  
    8989     REAL(wp), DIMENSION(jpi,jpj) ::   zhtc3   !: heat content of first 300 m                    [W] 
    9090 
    91 <<<<<<< variant A 
    9291     IF (iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1")) THEN 
    9392        ! ------------------------------------------------------------- ! 
     
    9998        zmaxdzT(:,:) = 0._wp 
    10099        IF( nla10 > 1 ) THEN 
    101 >>>>>>> variant B 
    102      IF(iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1").OR.iom_use("mld_dt02") & 
    103           & .OR.iom_use("topthdep").OR.iom_use("mldr10_3").OR.iom_use("pycndep").OR.iom_use("tinv").OR.iom_use("depti")) THEN 
    104         ! initialization 
    105         ztinv  (:,:) = 0._wp   
    106         zdepinv(:,:) = 0._wp   
    107         zmaxdzT(:,:) = 0._wp   
    108         DO jj = 1, jpj 
    109            DO ji = 1, jpi 
    110               zztmp = gdepw_n(ji,jj,mbkt(ji,jj)+1)  
    111               zhth     (ji,jj) = zztmp 
    112               zabs2   (ji,jj) = zztmp 
    113               ztm2    (ji,jj) = zztmp 
    114               zrho10_3(ji,jj) = zztmp 
    115               zpycn   (ji,jj) = zztmp 
    116            END DO 
    117         END DO 
    118         IF( nla10 > 1 ) THEN  
    119 ======= end 
    120100           DO jj = 1, jpj 
    121101              DO ji = 1, jpi 
     
    123103                 zrho0_3(ji,jj) = zztmp 
    124104                 zrho0_1(ji,jj) = zztmp 
    125 <<<<<<< variant A 
    126105                 zhth(ji,jj) = zztmp 
    127 >>>>>>> variant B 
    128 ======= end 
    129               END DO 
    130            END DO 
    131 <<<<<<< variant A 
     106              END DO 
     107           END DO 
    132108        ELSE IF (iom_use("mlddzt")) THEN 
    133109           DO jj = 1, jpj 
     
    139115        ELSE 
    140116           zhth(:,:) = 0._wp 
    141 >>>>>>> variant B 
    142 ======= end 
     117 
    143118        ENDIF 
    144 <<<<<<< variant A 
    145 >>>>>>> variant B 
    146      ENDIF 
    147 ======= end 
    148  
    149 <<<<<<< variant A 
    150 >>>>>>> variant B 
    151      IF (iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1")) THEN 
    152         ! ------------------------------------------------------------- ! 
    153         ! thermocline depth: strongest vertical gradient of temperature ! 
    154         ! turbocline depth (mixing layer depth): avt = zavt5            ! 
    155         ! MLD: rho = rho(1) + zrho3                                     ! 
    156         ! MLD: rho = rho(1) + zrho1                                     ! 
    157         ! ------------------------------------------------------------- ! 
    158 ======= end 
     119 
    159120        DO jk = jpkm1, 2, -1   ! loop from bottom to 2 
    160121           DO jj = 1, jpj 
     
    165126                 zzdep = zzdep * tmask(ji,jj,1) 
    166127 
    167 <<<<<<< variant A 
    168128                 IF( zztmp > zmaxdzT(ji,jj) ) THEN 
    169 >>>>>>> variant B 
    170                  IF( zztmp > zmaxdzT(ji,jj) ) THEN                         
    171 ======= end 
    172129                    zmaxdzT(ji,jj) = zztmp   ;   zhth    (ji,jj) = zzdep                ! max and depth of dT/dz 
    173130                 ENDIF 
    174131 
    175 <<<<<<< variant A 
    176132                 IF( nla10 > 1 ) THEN 
    177 >>>>>>> variant B 
    178                  IF( nla10 > 1 ) THEN  
    179 ======= end 
    180133                    zztmp = rhop(ji,jj,jk) - rhop(ji,jj,1)                             ! delta rho(1) 
    181134                    IF( zztmp > zrho3 )          zrho0_3(ji,jj) = zzdep                ! > 0.03 
     
    188141 
    189142        IF (iom_use("mlddzt")) CALL iom_put( "mlddzt", zhth*tmask(:,:,1) )            ! depth of the thermocline 
    190 <<<<<<< variant A 
    191143        IF( nla10 > 1 ) THEN 
    192 >>>>>>> variant B 
    193         IF( nla10 > 1 ) THEN  
    194 ======= end 
    195144           IF (iom_use("mldr0_3")) CALL iom_put( "mldr0_3", zrho0_3*tmask(:,:,1) )   ! MLD delta rho(surf) = 0.03 
    196145           IF (iom_use("mldr0_1")) CALL iom_put( "mldr0_1", zrho0_1*tmask(:,:,1) )   ! MLD delta rho(surf) = 0.01 
     
    360309     ENDIF 
    361310     ! 
    362       IF( ln_timing )   CALL timing_stop('dia_hth') 
    363      ! 
    364311   END SUBROUTINE dia_hth 
    365312 
Note: See TracChangeset for help on using the changeset viewer.