Changeset 3545


Ignore:
Timestamp:
2012-11-14T10:04:08+01:00 (8 years ago)
Author:
gm
Message:

branch: dev_r3385_NOCS04_HAMF; #665. bug correction in LIM3

Location:
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limdia.F90

    r3525 r3545  
    238238            ij0 = 136   ;   ij1 = 136              ! Fram strait : 83 = area export 
    239239            ii0 = 134   ;   ii1 = 138              !               84 = volume export 
    240             DO jj = mj0(ij0):mj1(ij1) 
    241                DO ji = mi0(ii0):mi1(ii1) 
     240            DO jj = mj0(ij0),mj1(ij1) 
     241               DO ji = mi0(ii0),mi1(ii1) 
    242242                  vinfor(83) = vinfor(83) - v_ice(ji,jj) * e1t(ji,jj)*at_i(ji,jj)*rdt_ice * 1.e-12_wp 
    243243                  vinfor(84) = vinfor(84) - v_ice(ji,jj) * e1t(ji,jj)*vt_i(ji,jj)*rdt_ice * 1.e-12_wp 
    244244               END DO 
    245245            END DO 
     246         END SELECT 
    246247!!gm   just above, this is NOT the correct way of evaluating the transport ! 
    247248!!gm        mass of snow is missing and v_ice should be the mean between jj and jj+1 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

    r3524 r3545  
    103103      INTEGER  ::   iflt, ial , iadv , ifral, ifrdv   !   -      - 
    104104      REAL(wp) ::   zinda, zemp, zemp_snow, zfmm      ! local scalars 
     105      REAL(wp) ::   zemp_snw                          !   -      - 
    105106      REAL(wp) ::   zfcm1 , zfcm2                     !   -      - 
    106107      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb, zalbp     ! 2D/3D workspace 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

    r3524 r3545  
    593593         ! 
    594594         ! Salt flux 
    595          sfx_thd_1d(ji) = sfx_thd_1d(ji) +        zihgnew  * zsfx_melt(ji)                & 
    596             &                            - (1.0 - zihgnew) * zfmass_i (ji) * sm_i_b(ji) ) * r1_rdtice 
     595         sfx_thd_1d(ji) = sfx_thd_1d(ji) +        zihgnew  * zsfx_melt(ji)               & 
     596            &                            - (1.0 - zihgnew) * zfmass_i (ji) * sm_i_b(ji) * r1_rdtice 
    597597         ! 
    598598         ! Heat flux 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r3517 r3545  
    178178         d_oa_i_thd (:,:,:)   = 0._wp   ;   d_oa_i_trp (:,:,:)   = 0._wp 
    179179         ! 
    180          fseqv    (:,:) = 0._wp 
    181          fsbri    (:,:) = 0._wp     ;   fsalt_res(:,:) = 0._wp 
    182          fsalt_rpo(:,:) = 0._wp 
    183          fhmec    (:,:) = 0._wp     ;   fhbri    (:,:) = 0._wp 
    184          fmmec    (:,:) = 0._wp     ;   fheat_res(:,:) = 0._wp 
    185          fheat_rpo(:,:) = 0._wp     ;   focea2D  (:,:) = 0._wp 
    186          fsup2D   (:,:) = 0._wp 
     180         sfx    (:,:) = 0._wp 
     181         sfx_bri(:,:) = 0._wp   ;   sfx_mec  (:,:) = 0._wp   ;   sfx_res  (:,:) = 0._wp 
     182         fhbri  (:,:) = 0._wp   ;   fheat_mec(:,:) = 0._wp   ;   fheat_res(:,:) = 0._wp 
     183         fhmec  (:,:) = 0._wp   ;    
     184         fmmec  (:,:) = 0._wp      
     185         focea2D(:,:) = 0._wp 
     186         fsup2D (:,:) = 0._wp 
    187187         !  
    188188         diag_sni_gr(:,:) = 0._wp   ;   diag_lat_gr(:,:) = 0._wp 
     
    428428               WRITE(numout,*) ' fdtcn     : ', fdtcn(ji,jj)  
    429429               WRITE(numout,*) ' fhmec     : ', fhmec(ji,jj)  
    430                WRITE(numout,*) ' fheat_rpo : ', fheat_rpo(ji,jj)  
     430               WRITE(numout,*) ' fheat_mec : ', fheat_mec(ji,jj)  
    431431               WRITE(numout,*) ' fheat_res : ', fheat_res(ji,jj)  
    432432               WRITE(numout,*) ' fhbri     : ', fhbri(ji,jj)  
     
    605605        WRITE(numout,*) ' - Heat / FW fluxes ' 
    606606        WRITE(numout,*) '   ~~~~~~~~~~~~~~~~ ' 
    607 !       WRITE(numout,*) ' fsbri      : ', fsbri(ki,kj) 
    608 !       WRITE(numout,*) ' fseqv      : ', fseqv(ki,kj) 
     607!       WRITE(numout,*) ' sfx_bri    : ', sfx_bri  (ki,kj) 
     608!       WRITE(numout,*) ' sfx        : ', sfx      (ki,kj) 
    609609!       WRITE(numout,*) ' fsalt_res  : ', fsalt_res(ki,kj) 
    610         WRITE(numout,*) ' fmmec      : ', fmmec(ki,kj) 
    611         WRITE(numout,*) ' fhmec      : ', fhmec(ki,kj) 
    612         WRITE(numout,*) ' fhbri      : ', fhbri(ki,kj) 
    613         WRITE(numout,*) ' fheat_rpo  : ', fheat_rpo(ki,kj) 
     610        WRITE(numout,*) ' fmmec      : ', fmmec    (ki,kj) 
     611        WRITE(numout,*) ' fhmec      : ', fhmec    (ki,kj) 
     612        WRITE(numout,*) ' fhbri      : ', fhbri    (ki,kj) 
     613        WRITE(numout,*) ' fheat_mec  : ', fheat_mec(ki,kj) 
    614614        WRITE(numout,*)  
    615615        WRITE(numout,*) ' sst        : ', sst_m(ki,kj)   
     
    620620        WRITE(numout,*) ' utau_ice   : ', utau_ice(ki,kj)  
    621621        WRITE(numout,*) ' vtau_ice   : ', vtau_ice(ki,kj) 
    622         WRITE(numout,*) ' utau       : ', utau(ki,kj)  
    623         WRITE(numout,*) ' vtau       : ', vtau(ki,kj) 
    624         WRITE(numout,*) ' oc. vel. u : ', u_oce(ki,kj) 
    625         WRITE(numout,*) ' oc. vel. v : ', v_oce(ki,kj) 
     622        WRITE(numout,*) ' utau       : ', utau    (ki,kj)  
     623        WRITE(numout,*) ' vtau       : ', vtau    (ki,kj) 
     624        WRITE(numout,*) ' oc. vel. u : ', u_oce   (ki,kj) 
     625        WRITE(numout,*) ' oc. vel. v : ', v_oce   (ki,kj) 
    626626     ENDIF 
    627627 
     
    639639        WRITE(numout,*) 
    640640        WRITE(numout,*) ' - Heat fluxes at bottom interface ***' 
    641         WRITE(numout,*) ' qsr        : ', qsr(ki,kj) 
    642         WRITE(numout,*) ' qns        : ', qns(ki,kj) 
     641        WRITE(numout,*) ' qsr       : ', qsr(ki,kj) 
     642        WRITE(numout,*) ' qns       : ', qns(ki,kj) 
    643643        WRITE(numout,*) 
    644644        WRITE(numout,*) ' - Salt fluxes at bottom interface ***' 
    645         WRITE(numout,*) ' sfx        : ', sfx (ki,kj) 
    646         WRITE(numout,*) ' emp        : ', emp(ki,kj) 
    647         WRITE(numout,*) ' fsbri      : ', fsbri(ki,kj) 
    648         WRITE(numout,*) ' fseqv      : ', fseqv(ki,kj) 
    649         WRITE(numout,*) ' fsalt_res  : ', fsalt_res(ki,kj) 
    650         WRITE(numout,*) ' fsalt_rpo  : ', fsalt_rpo(ki,kj) 
     645        WRITE(numout,*) ' emp       : ', emp    (ki,kj) 
     646        WRITE(numout,*) ' sfx_bri   : ', sfx_bri(ki,kj) 
     647        WRITE(numout,*) ' sfx       : ', sfx    (ki,kj) 
     648        WRITE(numout,*) ' sfx_res   : ', sfx_res(ki,kj) 
     649        WRITE(numout,*) ' sfx_mec   : ', sfx_mec(ki,kj) 
    651650        WRITE(numout,*) ' - Heat fluxes at bottom interface ***' 
    652         WRITE(numout,*) ' fheat_res  : ', fheat_res(ki,kj) 
     651        WRITE(numout,*) ' fheat_res : ', fheat_res(ki,kj) 
    653652        WRITE(numout,*) 
    654653        WRITE(numout,*) ' - Momentum fluxes ' 
Note: See TracChangeset for help on using the changeset viewer.