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

Changeset 15445


Ignore:
Timestamp:
2021-10-26T11:01:34+02:00 (10 months ago)
Author:
hadjt
Message:

TRA/tradwl.F90

Bug fix... used tmask(0,:,:) rather than tmask(1,:,:) this assumed everywhere was land. The radiation then penetrated too deep, and so the surface water as too cold.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/TRA/tradwl.F90

    r15318 r15445  
    8080      REAL(wp) ::   zc0, zc1, zc2, zc3   !    -         - 
    8181      !JT 
    82       REAL(wp) ::   hbatt 
     82      REAL(wp), DIMENSION(jpi,jpj) ::   hbatt, qsr_tradwl 
    8383      !JT 
    84       REAL(wp), DIMENSION(jpi,jpj)     ::   zekb, zekg, zekr            ! 2D workspace 
    85       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   ze0, ze1 , ze2, ze3, zea    ! 3D workspace 
    8684      !!---------------------------------------------------------------------- 
    8785      !! HERE GO VARIABLES USED IN POLCOMS CLEAN UP LATER 
     
    120118! 
    121119!--------------------------------------------------------------------------- 
    122       cp=3986.0d0 
     120    !cp=3986.0d0 
    123121 
    124122    DO jj = 2, jpj 
    125123         DO ji = fs_2, fs_jpim1 
    126            qsr(ji,jj)  = qsr(ji,jj)  * (r1_rau0_rcp) 
     124           qsr_tradwl(ji,jj)  = qsr(ji,jj)  * (r1_rau0_rcp) 
    127125         ENDDO       !ji 
    128126    ENDDO            !jj 
    129127!-------------------------------------------------------------------------------- 
    130       if ( first ) then 
     128  
     129 
     130   if ( first ) then 
     131    do jj=2,jpjm1 
     132      do ji = fs_2, fs_jpim1  
     133          IF ( tmask(ji,jj,1) .EQ. 1) THEN ! if land 
     134            hbatt(ji,jj) = sum( e3t_n(ji,jj,:)*tmask(ji,jj,:) ) 
     135        else 
     136            hbatt(ji,jj)= 0. 
     137        endif 
     138      enddo ! ji 
     139    enddo ! jj 
     140 
     141   !CALL iom_put('hbatt_tradwl', hbatt(:,:) ) 
     142 
    131143        rlambda2(:,:) = 0.0 
    132144        first=.false. 
    133145        if ( ln_vary_lambda ) then 
    134 !        do j=1,jesub    ! Original Polcoms style Loop 
    135 !          do i=1,iesub  ! Original Polcoms style Loop 
    136146 
    137147        do jj=2,jpjm1 
    138148          do ji = fs_2, fs_jpim1   ! vector opt.  
    139               if (tmask(ji,jj,0) == 1)  then    
    140  
    141 !             if(ipexb(i,j).ne. 0) then  (Mask, use Tmask instead) 
    142  
    143  
    144               !JT 
    145               !hbatt = gdept_n(ji,jj, k_bot(ji,jj) )  
    146               hbatt = sum( e3t_n(ji,jj,:)*tmask(ji,jj,:) ) 
    147  
    148               rlambda2(ji,jj)=-0.033*log(hbatt)+0.2583    ! JIAs formula 
    149               !JT 
    150  
    151  
    152               !JT rlambda2(ji,jj)=-0.033*log(hbatt(ji,jj))+0.2583    ! JIAs formula 
     149              !IF ( tmask(ji,jj,1) .EQ. 1) THEN ! if land 
     150 
     151 
     152              rlambda2(ji,jj)=-0.033*log(hbatt(ji,jj))+0.2583    ! JIAs formula 
    153153              rlambda2(ji,jj)=max(0.05,rlambda2(ji,jj))     ! limit in deep water 
    154154              rlambda2(ji,jj)=min(0.25,rlambda2(ji,jj))     ! Catch the infinities, from very shallow water/land. 10cm = 0.25 
    155155 
    156               !WRITE(*,300) 'JT tradwl:',jj,ji,njmpp,jpjglo,nimpp,jpiglo,narea, hbatt, rlambda2(ji,jj) 
    157 !300 FORMAT(A14,1X,I4,1X,I4,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,f9.3,1X,f9.2) 
    158  
    159  
    160 !              WRITE(*,300) 'JT tradwl:',jj,ji,njmpp+jj,nimpp+ji,njmpp,nimpp,narea, hbatt, rlambda2(ji,jj) 
    161                !domain size jpjglo,,jpiglo 
    162                !lower lhs of each sub-domain = nimpp,njmpp 
    163                ! index on the global domain??? add or subtract one?? = njmpp+jj,nimpp+ji 
    164 !300 FORMAT(A14,1X,I4,1X,I4,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,f9.3,1X,f9.2) 
    165  
    166  
    167 !              if (kt == 1) WRITE(*,300) 'JT tradwl:',njmpp+jj,nimpp+ji, hbatt, rlambda2(ji,jj) 
    168                !domain size jpjglo,,jpiglo 
    169                !lower lhs of each sub-domain = nimpp,njmpp 
    170                ! index on the global domain??? add or subtract one?? = njmpp+jj,nimpp+ji 
    171 !300 FORMAT(A14,1X,I4,1X,I4,1X,f9.3,1X,f9.2) 
    172  
    173  
    174             else 
    175                 rlambda2(ji,jj)= 0.25 
    176             endif 
     156            !else 
     157            !    rlambda2(ji,jj)= 0.25 
     158            !endif 
    177159          enddo ! ji 
    178160        enddo ! jj 
     
    183165      endif ! If first 
    184166 
    185 !      do j=1,jesub      ! Original Polcoms Style Loop 
    186 !        do i=1,iesub    ! Original Polcoms Style Loop 
     167      ! CALL iom_put('rlambda2_tradwl', rlambda2(:,:) ) 
     168 
    187169      DO jk=2,jpk 
    188170         DO jj=2,jpjm1 
    189171            DO ji = fs_2, fs_jpim1   ! vector opt. 
    190172 
    191               if (tmask(ji,jj,0) == 1)  then    
     173              IF ( tmask(ji,jj,1) .EQ. 1) THEN ! if land 
     174 
    192175    !-------------------------------------------------------------------- 
    193176    ! Calculate change in temperature 
     
    196179    !        rad0 = hfl_in(i,j)   ! change hfl_in to qsr I assume 
    197180 
    198                     rad0 = qsr(ji,jj) 
     181                    rad0 = qsr_tradwl(ji,jj) 
    199182                    rD = rLambda2(ji,jj)  +rLambda      !  Transmissivity to be used here 
    200     !       if rlambda 0 then rlambda2 not zer and vica versa  
     183                          !       if rlambda 0 then rlambda2 not zer and vica versa  
    201184 
    202185                    z2=gdepw_0(ji,jj,jk-1)    ! grid box is from z=z1 to z=z2 
     
    209192                    dtmp(jk)=1.0/(e3t_0(ji,jj,jk))*(Rad2-Rad1) !change in temperature 
    210193                    tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) + dtmp(jk) 
    211                 endif 
     194                endif ! if land 
    212195            enddo  ! ji 
    213196         enddo  ! jj 
     
    249232      REAL(wp) ::   zc2  , zc3  , zchl    !    -         - 
    250233      REAL(wp) ::   zsi0r, zsi1r, zcoef   !    -         - 
    251       REAL(wp), DIMENSION(jpi,jpj)     ::   zekb, zekg, zekr              ! 2D workspace 
    252       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   ze0 , ze1 , ze2 , ze3 , zea   ! 3D workspace 
    253234      !! 
    254235      CHARACTER(len=100) ::   cn_dir   ! Root directory for location of ssr files 
Note: See TracChangeset for help on using the changeset viewer.