Ignore:
Timestamp:
2015-04-17T12:24:12+02:00 (6 years ago)
Author:
cetlod
Message:

dev_r5204_CNRS_PISCES_dcy: some improvments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

    r5207 r5222  
    219219 
    220220 
     221      IF( lk_cpl ) THEN  
     222         !!! LIM2 version zqsr = qsr_tot(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj,1) ) * ( 1.0 - pfrld(ji,jj) ) 
     223          zfcm1 = qsr_tot(ji,jj) 
     224          DO jl = 1, jpl 
     225             zfcm1 = zfcm1 + ( ftr_ice(ji,jj,jl) - qsr_ice(ji,jj,jl) ) * a_i_b(ji,jj,jl) 
     226          END DO 
     227      ELSE 
     228          !!! LIM2 version zqsr = pfrld(ji,jj) * qsr(ji,jj)  + ( 1.  - pfrld(ji,jj) ) * fstric(ji,jj) 
     229          zfcm1   = pfrld(ji,jj) * qsr(ji,jj) 
     230          DO jl = 1, jpl 
     231             zfcm1   = zfcm1 + a_i_b(ji,jj,jl) * ftr_ice(ji,jj,jl) 
     232          END DO 
     233      ENDIF 
     234 
    221235      !   daily mean qsr when diurnal cycle is applied on physics - for BGC models 
    222236      IF( l_trcdm2dc ) THEN 
    223          DO jj = 1, jpj 
    224             DO ji = 1, jpi 
    225                zfcm1 = pfrld(ji,jj) * qsr_mean(ji,jj) 
    226                DO jl = 1, jpl 
    227                   zfcm1 = zfcm1 + ftr_ice_mean(ji,jj,jl) 
    228                END DO 
    229                qsr_mean(ji,jj) = zfcm1                                       
     237         IF( lk_cpl ) THEN  
     238            DO jj = 1, jpj 
     239               DO ji = 1, jpi 
     240                  zfcm1 = qsr_mean(ji,jj) 
     241                  DO jl = 1, jpl 
     242                    zfcm1 = zfcm1 + ftr_ice_mean(ji,jj,jl) - qsr_ice_mean(ji,jj,jl) * a_i_b(ji,jj,jl) 
     243                 END DO 
     244                 qsr_mean(ji,jj) = zfcm1                                       
     245               ENDDO 
    230246            ENDDO 
    231          ENDDO 
     247         ELSE 
     248            DO jj = 1, jpj 
     249               DO ji = 1, jpi 
     250                  zfcm1 = pfrld(ji,jj) * qsr_mean(ji,jj) 
     251                  DO jl = 1, jpl 
     252                     zfcm1 = zfcm1 + ftr_ice_mean(ji,jj,jl) 
     253                  END DO 
     254                  qsr_mean(ji,jj) = zfcm1                                       
     255               ENDDO 
     256            ENDDO 
     257         ENDIF 
    232258      ENDIF 
    233259 
Note: See TracChangeset for help on using the changeset viewer.