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 2239 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/ASM – NEMO

Ignore:
Timestamp:
2010-10-13T00:03:11+02:00 (14 years ago)
Author:
rblod
Message:

First commit to compile ORCA2_LIM for next nemo v3.3

Location:
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/ASM
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/ASM/asminc.F90

    r2218 r2239  
    4747      & ub, un, ua,    & 
    4848      & vb, vn, va,    & 
    49       & tb, tn, ta,    & 
    50       & sb, sn, sa,    & 
     49      & tsb, tsn, tsa, & 
    5150      & sshb, sshn,    & 
    5251      & rhd, rhop,     & 
    5352      & rotb, rotn,    & 
    5453      & hdivb, hdivn,  & 
    55       & gtu, gsu, gru, & 
    56       & gtv, gsv, grv  
     54      & gtsu, gru,    & 
     55      & gtsv, grv  
    5756   USE divcur, ONLY : &   ! Horizontal divergence and relative vorticity 
    5857      & div_cur 
     
    724723            ! Update the tracer tendencies 
    725724            DO jk = 1, jpkm1 
    726                ta(:,:,jk) = ta(:,:,jk) + t_bkginc(:,:,jk) * zincwgt   
    727                sa(:,:,jk) = sa(:,:,jk) + s_bkginc(:,:,jk) * zincwgt 
     725               tsa(:,:,jk,jp_tem) = tsa(:,:,jk,jp_tem) + t_bkginc(:,:,jk) * zincwgt   
     726               tsa(:,:,jk,jp_sal) = tsa(:,:,jk,jp_sal) + s_bkginc(:,:,jk) * zincwgt 
    728727            END DO 
    729728 
     
    733732                  DO jj = 1, jpj 
    734733                     DO ji= 1, jpi 
    735                         sa(ji,jj,jk) = MAX( sa(ji,jj,jk), salfixmin ) 
     734                        tsa(ji,jj,jk,jp_sal) = MAX( tsa(ji,jj,jk,jp_sal), salfixmin ) 
    736735                     END DO 
    737736                  END DO 
     
    758757 
    759758            ! Initialize the now fields with the background + increment 
    760             tn(:,:,:) = t_bkg(:,:,:) + t_bkginc(:,:,:)    
    761             sn(:,:,:) = s_bkg(:,:,:) + s_bkginc(:,:,:)    
     759            tsn(:,:,:,jp_tem) = t_bkg(:,:,:) + t_bkginc(:,:,:)    
     760            tsn(:,:,:,jp_sal) = s_bkg(:,:,:) + s_bkginc(:,:,:)    
    762761 
    763762            ! Optional salinity fix 
     
    766765                  DO jj = 1, jpj 
    767766                     DO ji= 1, jpi 
    768                         sn(ji,jj,jk) = MAX( sn(ji,jj,jk), salfixmin ) 
     767                        tsn(ji,jj,jk,jp_sal) = MAX( tsn(ji,jj,jk,jp_sal), salfixmin ) 
    769768                     END DO 
    770769                  END DO 
     
    772771            ENDIF 
    773772 
    774             tb(:,:,:) = tn(:,:,:)                        ! Update before fields 
    775             sb(:,:,:) = sn(:,:,:) 
    776  
    777             CALL eos( tb, sb, rhd, rhop )                ! Before potential and in situ densities 
     773            tsb(:,:,:,:) = tsn(:,:,:,:)                        ! Update before fields 
     774 
     775            CALL eos( tsb, rhd, rhop )                ! Before potential and in situ densities 
    778776          
    779777#if ! defined key_c1d 
    780778            IF( ln_zps ) & 
    781                &  CALL zps_hde( nit000, tb, sb, rhd,  &  ! Partial steps: before horizontal derivative 
    782                &                gtu, gsu, gru,        &  ! of T, S, rd at the bottom ocean level 
    783                &                gtv, gsv, grv ) 
     779               &  CALL zps_hde( nit000, jpts, tsb,   &  ! Partial steps: before horizontal derivative 
     780               &                gtsu, gtsv, rhd,        &  ! of T, S, rd at the bottom ocean level 
     781               &                gru , grv ) 
    784782#endif 
    785783 
  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/ASM/asmtrj.F90

    r2128 r2239  
    224224#if defined key_ldfslp 
    225225            IF(lwp) WRITE(numout,*) ' Compute the slopes of neutral surface...' 
    226             CALL bn2( tb, sb, rn2 )  
     226            CALL bn2( tsb, rn2 )  
    227227            CALL ldf_slp( nit000, rhd, rn2 ) 
    228228#endif 
     
    280280 
    281281END MODULE asmtrj 
    282     
    283     
Note: See TracChangeset for help on using the changeset viewer.