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

original lot of changes to diahth.F90

File:
1 edited

Legend:

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

    r11091 r11093  
    8989     REAL(wp), DIMENSION(jpi,jpj) ::   zhtc3   !: heat content of first 300 m                    [W] 
    9090 
     91<<<<<<< variant A 
    9192     IF (iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1")) THEN 
    9293        ! ------------------------------------------------------------- ! 
     
    9899        zmaxdzT(:,:) = 0._wp 
    99100        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 
    100120           DO jj = 1, jpj 
    101121              DO ji = 1, jpi 
     
    103123                 zrho0_3(ji,jj) = zztmp 
    104124                 zrho0_1(ji,jj) = zztmp 
     125<<<<<<< variant A 
    105126                 zhth(ji,jj) = zztmp 
    106               END DO 
    107            END DO 
     127>>>>>>> variant B 
     128======= end 
     129              END DO 
     130           END DO 
     131<<<<<<< variant A 
    108132        ELSE IF (iom_use("mlddzt")) THEN 
    109133           DO jj = 1, jpj 
     
    115139        ELSE 
    116140           zhth(:,:) = 0._wp 
     141>>>>>>> variant B 
     142======= end 
    117143        ENDIF 
    118  
     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 
    119159        DO jk = jpkm1, 2, -1   ! loop from bottom to 2 
    120160           DO jj = 1, jpj 
     
    125165                 zzdep = zzdep * tmask(ji,jj,1) 
    126166 
     167<<<<<<< variant A 
    127168                 IF( zztmp > zmaxdzT(ji,jj) ) THEN 
     169>>>>>>> variant B 
     170                 IF( zztmp > zmaxdzT(ji,jj) ) THEN                         
     171======= end 
    128172                    zmaxdzT(ji,jj) = zztmp   ;   zhth    (ji,jj) = zzdep                ! max and depth of dT/dz 
    129173                 ENDIF 
    130174 
     175<<<<<<< variant A 
    131176                 IF( nla10 > 1 ) THEN 
     177>>>>>>> variant B 
     178                 IF( nla10 > 1 ) THEN  
     179======= end 
    132180                    zztmp = rhop(ji,jj,jk) - rhop(ji,jj,1)                             ! delta rho(1) 
    133181                    IF( zztmp > zrho3 )          zrho0_3(ji,jj) = zzdep                ! > 0.03 
     
    140188 
    141189        IF (iom_use("mlddzt")) CALL iom_put( "mlddzt", zhth*tmask(:,:,1) )            ! depth of the thermocline 
     190<<<<<<< variant A 
    142191        IF( nla10 > 1 ) THEN 
     192>>>>>>> variant B 
     193        IF( nla10 > 1 ) THEN  
     194======= end 
    143195           IF (iom_use("mldr0_3")) CALL iom_put( "mldr0_3", zrho0_3*tmask(:,:,1) )   ! MLD delta rho(surf) = 0.03 
    144196           IF (iom_use("mldr0_1")) CALL iom_put( "mldr0_1", zrho0_1*tmask(:,:,1) )   ! MLD delta rho(surf) = 0.01 
     
    308360     ENDIF 
    309361     ! 
     362      IF( ln_timing )   CALL timing_stop('dia_hth') 
     363     ! 
    310364   END SUBROUTINE dia_hth 
    311365 
Note: See TracChangeset for help on using the changeset viewer.