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

Changeset 15085


Ignore:
Timestamp:
2021-07-05T22:07:24+02:00 (3 years ago)
Author:
clem
Message:

ISF further cleaning and debugging(?) for isfcpl but I think it requires more attention

Location:
NEMO/trunk/src/OCE/ISF
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/OCE/ISF/isfcav.F90

    r15084 r15085  
    109109            zRc(ji,jj) = MAX(rn2(ji,jj,ikt+1), 1.e-20_wp) / MAX( zdku*zdku + zdkv*zdkv, 1.e-20_wp ) 
    110110         END_2D 
    111          CALL lbc_lnk( 'isfmlt', zRc, 'T', 1._wp ) 
     111         CALL lbc_lnk( 'isfcav', zRc, 'T', 1._wp ) 
    112112      ENDIF 
    113113      ! 

  • NEMO/trunk/src/OCE/ISF/isfcpl.F90

    r15062 r15085  
    201201            ENDIF 
    202202         END_2D 
     203         CALL lbc_lnk( 'isfcpl', ssh(:,:,Kmm), 'T', 1.0_wp, zssmask_b(:,:), 'T', 1.0_wp ) 
    203204         ! 
    204205         zssh(:,:) = ssh(:,:,Kmm) 
    205206         zssmask0(:,:) = zssmask_b(:,:) 
    206207         ! 
    207          CALL lbc_lnk( 'iscplrst', zssh, 'T', 1.0_wp, zssmask0, 'T', 1.0_wp ) 
    208208         ! 
    209209      END DO 
    … …  
    359359         END DO 
    360360         ! 
     361         CALL lbc_lnk( 'isfcpl', ts(:,:,:,jp_tem,Kmm), 'T', 1.0_wp, ts(:,:,:,jp_sal,Kmm), 'T', 1.0_wp, ztmask1, 'T', 1.0_wp) 
     362         ! 
    361363         ! update temperature and salinity and mask 
    362364         zts0(:,:,:,:)  = ts(:,:,:,:,Kmm) 
    363365         ztmask0(:,:,:) = ztmask1(:,:,:) 
    364366         ! 
    365          CALL lbc_lnk( 'iscplrst', zts0(:,:,:,jp_tem), 'T', 1.0_wp, zts0(:,:,:,jp_sal), 'T', 1.0_wp, ztmask0, 'T', 1.0_wp) 
    366367         ! 
    367368      END DO  ! nn_drown 
    … …  
    437438               &    - e1v(ji  ,jj-1) * e3v(ji  ,jj-1,jk,Kmm) * vv(ji  ,jj-1,jk,Kmm)  ) & 
    438439               &               * tmask(ji,jj,jk) 
     440            ! 
     441            ! 1.3: get 3d volume flux difference (before - after cpl) (>0 out) 
     442            !      correction to add is _b - _n 
     443            risfcpl_vol(ji,jj,jk) = zqvolb(ji,jj,jk) - zqvoln(ji,jj,jk) 
    439444         END_2D 
    440          ! 
    441          ! 1.3: get 3d volume flux difference (before - after cpl) (>0 out) 
    442          !      correction to add is _b - _n 
    443          risfcpl_vol(:,:,jk) = zqvolb(:,:,jk) - zqvoln(:,:,jk) 
    444445      END DO 
    445446      ! 
    … …  
    455456      END_2D 
    456457      ! 
    457       CALL lbc_lnk( 'iscpl', risfcpl_vol, 'T', 1.0_wp ) 
     458      CALL lbc_lnk( 'isfcpl', risfcpl_vol, 'T', 1.0_wp ) 
    458459      ! 
    459460      ! 3.0: set total correction (div, tr(:,:,:,:,Krhs), ssh) 
    … …  
    691692      ! 
    692693      ! add lbclnk 
    693       CALL lbc_lnk( 'iscplrst', risfcpl_cons_tsc(:,:,:,jp_tem), 'T', 1.0_wp, risfcpl_cons_tsc(:,:,:,jp_sal), 'T', 1.0_wp, & 
    694          &                      risfcpl_cons_vol(:,:,:)       , 'T', 1.0_wp) 
     694      CALL lbc_lnk( 'isfcpl', risfcpl_cons_tsc(:,:,:,jp_tem), 'T', 1.0_wp, risfcpl_cons_tsc(:,:,:,jp_sal), 'T', 1.0_wp, & 
     695         &                    risfcpl_cons_vol(:,:,:)       , 'T', 1.0_wp) 
    695696      ! 
    696697      ! ssh correction (for dynspg_ts) 
Note: See TracChangeset for help on using the changeset viewer.