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 7306 for branches/2016/dev_r6859_LIM3_meltponds – NEMO

Ignore:
Timestamp:
2016-11-22T14:57:00+01:00 (7 years ago)
Author:
vancop
Message:

First bugs found

Location:
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r7293 r7306  
    281281   LOGICAL , PUBLIC ::   ln_limMP         !: activate melt ponds or not 
    282282   INTEGER , PUBLIC ::   nn_limMP         !: type of melt pond implementation (1=full, 2=radiation only, 3=freshwater only) 
     283   ! END MV MP 2016 
    283284 
    284285   !                                     !!** some other parameters  
     
    527528      ALLOCATE( t_i(jpi,jpj,nlay_i,jpl) , e_i(jpi,jpj,nlay_i,jpl) , s_i(jpi,jpj,nlay_i,jpl) , STAT=ierr(ii) ) 
    528529 
     530      ! MV MP 2016 
     531      ii = ii + 1 
     532      ALLOCATE( a_ip(jpi,jpj,jpl) , v_ip(jpi,jpj,jpl) , a_ip_frac(jpi,jpj,jpl) , & 
     533         &      h_ip(jpi,jpj,jpl) , STAT = ierr(ii) ) 
     534      ii = ii + 1 
     535      ALLOCATE( at_ip(jpi,jpj) , vt_ip(jpi,jpj) , STAT = ierr(ii) ) 
     536      ! END MV MP 2016 
     537 
    529538      ! * Moments for advection 
    530539      ii = ii + 1 
     
    544553         &      syye(jpi,jpj,nlay_i,jpl) , sxye(jpi,jpj,nlay_i,jpl)                            , STAT=ierr(ii) ) 
    545554 
     555      ! MV MP 2016 
     556      ALLOCATE( sxap(jpi,jpj,jpl) , syap(jpi,jpj,jpl) , sxxap(jpi,jpj,jpl) , syyap(jpi,jpj,jpl) , sxyap(jpi,jpj,jpl) ,   & 
     557         &      sxvp(jpi,jpj,jpl) , syvp(jpi,jpj,jpl) , sxxvp(jpi,jpj,jpl) , syyvp(jpi,jpj,jpl) , sxyvp(jpi,jpj,jpl) ,   & 
     558         &      STAT = ierr(ii) ) 
     559      ! END MV MP 2016 
     560 
    546561      ! * Old values of global variables 
    547562      ii = ii + 1 
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90

    r7293 r7306  
    474474      ! MV MP 2016 
    475475      IF ( ln_limMP ) THEN 
    476          sxap  (:,:) = 0._wp    ; sxvp  (:,:) = 0._wp  
    477          syap  (:,:) = 0._wp    ; syvp  (:,:) = 0._wp  
    478          sxxap (:,:) = 0._wp    ; sxxvp (:,:) = 0._wp  
    479          syyap (:,:) = 0._wp    ; syyvp (:,:) = 0._wp  
    480          sxyap (:,:) = 0._wp    ; sxyvp (:,:) = 0._wp 
     476         sxap  (:,:,:) = 0._wp    ; sxvp  (:,:,:) = 0._wp  
     477         syap  (:,:,:) = 0._wp    ; syvp  (:,:,:) = 0._wp  
     478         sxxap (:,:,:) = 0._wp    ; sxxvp (:,:,:) = 0._wp  
     479         syyap (:,:,:) = 0._wp    ; syyvp (:,:,:) = 0._wp  
     480         sxyap (:,:,:) = 0._wp    ; sxyvp (:,:,:) = 0._wp 
    481481      ENDIF 
    482482      ! END MV MP 2016 
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90

    r7293 r7306  
    154154      IF ( ln_limMP ) THEN 
    155155         DO jl = 1, jpl  
     156            WRITE(zchar,'(I2.2)') jl 
    156157            znam = 'a_ip'//'_htc'//zchar 
    157158            z2d(:,:) = a_ip(:,:,jl) 
     
    317318         IF ( ln_limMP ) THEN 
    318319            DO jl = 1, jpl  
     320               WRITE(zchar,'(I2.2)') jl 
    319321               znam = 'sxap'//'_htc'//zchar 
    320322               z2d(:,:) = sxap(:,:,jl) 
     
    439441      IF ( ln_limMP ) THEN 
    440442         DO jl = 1, jpl  
     443            WRITE(zchar,'(I2.2)') jl 
    441444            znam = 'a_ip'//'_htc'//zchar 
    442445            CALL iom_get( numrir, jpdom_autoglo, znam , z2d ) 
     
    447450         END DO 
    448451      ENDIF 
     452      ! END MV MP 2016 
    449453 
    450454      DO jl = 1, jpl  
     
    573577         IF ( ln_limMP ) THEN 
    574578            DO jl = 1, jpl  
     579               WRITE(zchar,'(I2.2)') jl 
    575580               znam = 'sxap'//'_htc'//zchar 
    576581               CALL iom_get( numrir, jpdom_autoglo, znam , z2d ) 
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90

    r7293 r7306  
    9292 
    9393      ! MV MP 2016 
    94       IF ( ln_limMP ) 
    95          at_i_p(:,:) = SUM( a_ip, dim=3 ) 
    96          vt_i_p(:,:) = SUM( v_ip, dim=3 ) 
     94      IF ( ln_limMP ) THEN 
     95         at_ip(:,:) = SUM( a_ip, dim=3 ) 
     96         vt_ip(:,:) = SUM( v_ip, dim=3 ) 
    9797      ENDIF 
    9898      ! END MP 2016 
Note: See TracChangeset for help on using the changeset viewer.