Changeset 2239 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/ASM
- Timestamp:
- 2010-10-13T00:03:11+02:00 (14 years ago)
- 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 47 47 & ub, un, ua, & 48 48 & vb, vn, va, & 49 & tb, tn, ta, & 50 & sb, sn, sa, & 49 & tsb, tsn, tsa, & 51 50 & sshb, sshn, & 52 51 & rhd, rhop, & 53 52 & rotb, rotn, & 54 53 & hdivb, hdivn, & 55 & gt u, gsu, gru,&56 & gt v, gsv, grv54 & gtsu, gru, & 55 & gtsv, grv 57 56 USE divcur, ONLY : & ! Horizontal divergence and relative vorticity 58 57 & div_cur … … 724 723 ! Update the tracer tendencies 725 724 DO jk = 1, jpkm1 726 t a(:,:,jk) = ta(:,:,jk) + t_bkginc(:,:,jk) * zincwgt727 sa(:,:,jk) = sa(:,:,jk) + s_bkginc(:,:,jk) * zincwgt725 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 728 727 END DO 729 728 … … 733 732 DO jj = 1, jpj 734 733 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 ) 736 735 END DO 737 736 END DO … … 758 757 759 758 ! Initialize the now fields with the background + increment 760 t n(:,:,:) = 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(:,:,:) 762 761 763 762 ! Optional salinity fix … … 766 765 DO jj = 1, jpj 767 766 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 ) 769 768 END DO 770 769 END DO … … 772 771 ENDIF 773 772 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 778 776 779 777 #if ! defined key_c1d 780 778 IF( ln_zps ) & 781 & CALL zps_hde( nit000, tb, sb, rhd,& ! Partial steps: before horizontal derivative782 & gt u, gsu, gru, & ! of T, S, rd at the bottom ocean level783 & g tv, 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 ) 784 782 #endif 785 783 -
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/ASM/asmtrj.F90
r2128 r2239 224 224 #if defined key_ldfslp 225 225 IF(lwp) WRITE(numout,*) ' Compute the slopes of neutral surface...' 226 CALL bn2( t b,sb, rn2 )226 CALL bn2( tsb, rn2 ) 227 227 CALL ldf_slp( nit000, rhd, rn2 ) 228 228 #endif … … 280 280 281 281 END MODULE asmtrj 282 283
Note: See TracChangeset
for help on using the changeset viewer.