Changeset 48 for trunk/INTERP2/newinterp.pro
- Timestamp:
- 03/16/14 20:38:39 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/INTERP2/newinterp.pro
r2 r48 12 12 ; 13 13 ; INPUTS : 14 ; ORCA_GEO grid file + ORCA_GEO weight file + Input field 14 ; ORCA_GEO grid file + ORCA_GEO weight file + Input field 15 15 ; on ORCA grid 16 16 ; OUTPUTS : … … 108 108 109 109 for jj=0,jpj-1 do begin 110 111 110 111 112 112 if mask_reg(ji,jj,0) NE 0 then begin 113 113 … … 117 117 + weight(ji,jj,2)*input_field(pointer_i(ji,jj,2)+1,pointer_j(ji,jj,2),*,*) $ 118 118 + weight(ji,jj,3)*input_field(pointer_i(ji,jj,3)+1,pointer_j(ji,jj,3),*,*) 119 119 120 120 endif 121 121 … … 131 131 ; Data of the input grid is duplicated for the northern boundary condition 132 132 133 zdata_inv(0:jpi/2 -1 ,*,*,*) = zdata_north(jpi/2:jpi-1,*,*,*) 134 zdata_inv(jpi/2:jpi-1 ,*,*,*) = zdata_north(0:jpi/2-1,*,*,*) 133 zdata_inv(0:jpi/2 -1 ,*,*,*) = zdata_north(jpi/2:jpi-1,*,*,*) 134 zdata_inv(jpi/2:jpi-1 ,*,*,*) = zdata_north(0:jpi/2-1,*,*,*) 135 135 zdata_inv = reverse(zdata_inv,2) 136 136 zdata_north(*,jpj:jpj+jplus-1,*,*) = zdata_inv(*,jpj:jpj+jplus-1,*,*) … … 178 178 179 179 ; Attributes 180 180 181 181 id0 = NCDF_VARDEF(idout, 'nav_lon' , [xidout, yidout ],/FLOAT ) 182 182 id1 = NCDF_VARDEF(idout, 'nav_lat' , [xidout, yidout ],/FLOAT ) … … 207 207 208 208 209 209 210 210 NCDF_CONTROL, idout, /ENDEF 211 211 … … 222 222 NCDF_CLOSE, idout 223 223 224 225 224 225 226 226 227 227 end 228 229 230 231 232 233 234 235
Note: See TracChangeset
for help on using the changeset viewer.