source: trunk/INTERP2/make_geogrid.pro @ 48

Last change on this file since 48 was 48, checked in by pinsard, 10 years ago

fix thanks to coding rules
File size: 16.8 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;
5; NAME: make_geogrid.pro
6;
7; PURPOSE: Creates for a given ORCA grid, ORCA_GEO which is
8; the closest geographical grid to this ORCA grid
9;
10; CATEGORY : Subroutine
11;
12; CALLING SEQUENCE : make_geogrid
13;
14; INPUTS :
15;         All the variables of a given ORCA meshmask file
16; OUTPUTS :
17;         A netcdf format meshmask file that has the same size as the
18;         one given in input and contains exactly the same variable,
19;         although they describe another grid.
20;
21; COMMON BLOCKS:
22;         None
23;
24; SIDE EFFECTS:
25;
26; RESTRICTIONS:
27;
28; EXAMPLE:
29;
30; MODIFICATION HISTORY: 08/2002 Robinson Hordoir
31;
32;------------------------------------------------------------
33;------------------------------------------------------------
34;------------------------------------------------------------
35
36
37pro make_geogrid
38
39@initorca
40@init_path2
41@naminterp2
42
43;------------------------------------------------------------
44; Creation of a regular grid based on a given ORCA grid
45;------------------------------------------------------------
46
47
48;******************************************
49; 0. DEFINE OUTPUT VALUES
50;******************************************
51
52
53; Output values are set to input values
54
55glamt_reg=glamt
56glamu_reg=glamu
57glamv_reg=glamv
58glamf_reg=glamf
59
60gphit_reg=gphit
61gphiu_reg=gphiu
62gphiv_reg=gphiv
63gphif_reg=gphif
64
65e1t_reg=e1t
66e1u_reg=e1u
67e1v_reg=e1v
68e1f_reg=e1f
69
70e2t_reg=e2t
71e2u_reg=e2u
72e2v_reg=e2v
73e2f_reg=e2f
74
75tmask_reg=tmask
76umask_reg=tmask
77vmask_reg=tmask
78fmask_reg=tmask
79
80
81;*****************************************************************
82; 1. FIRST STEP : DEFINE  THE LATITUDE, LONGITUDE & SCALE FACTORS
83;*****************************************************************
84
85; Let's get our reference latitude. Search for northernmost place
86; of the ORCA grid that is still geographical
87
88northernmost=where(gphit EQ max(gphit))
89
90jnorth=northernmost/jpi
91
92inorth=min(northernmost-jnorth*jpi)
93
94ref_latitude=gphit(inorth,*)
95
96; Let's get the reference longitude
97
98equador=where(abs(gphit) EQ min(abs(gphit)) )
99
100ref_longitude=glamt(equador)
101
102; Duplication for all T points
103
104glamt_reg=ref_longitude#replicate(1,n_elements(ref_latitude))
105
106
107gphit_reg=replicate(1,n_elements(ref_longitude))#ref_latitude
108
109; COMPUTATION OF LONGITUDES AND LATITUDES AT ALL POINTS
110
111; From this we can deduce glamu_reg and gphiu_reg
112
113glamu_reg=glamt_reg+grid_def/2.
114gphiu_reg=gphit_reg
115
116; As well as glamv_reg and gphiv_reg
117
118glamv_reg=glamt_reg
119gphiv_reg=gphit_reg+(gphit_reg-shift(gphit_reg,0,-1))/2.
120gphiv_reg(*,jpj-1)=gphit_reg(*,jpj-1)+(gphit_reg(*,jpj-1)-gphit_reg(*,jpj-2))/2.
121
122; Same thing for glamf_reg and gphif_reg
123
124glamf_reg=glamu_reg
125gphif_reg=gphiv_reg
126
127
128; COMPUTATION OF SCALE FACTORS
129
130
131
132; Scale factor at equator
133
134eq_j=where(gphit(inorth,*) EQ 0.)
135
136e1t_eq=e1t(*,eq_j)
137
138; E1t scale factor
139
140e1t_reg=( e1t_eq#replicate(1.,jpj) )*(replicate(1.,jpi)#cos(!PI*gphit(inorth,*)/180.))
141
142
143; E2t scale factor
144
145e2t_reg=replicate(1.,jpi)#e2t(inorth,*)
146
147
148; E1u scale factor
149; e1u_reg=e1t_reg because the grid is fully regular along x axis
150; Same remark for e2u
151
152e1u_reg=e1t_reg
153e2u_reg=e2t_reg
154
155;E1v scale factor
156
157e1v_reg=( e1t_eq#replicate(1.,jpj) )*(replicate(1.,jpi)#cos(!PI*gphiv(inorth,*)/180.))
158
159; E2v scale factor
160
161e2v_reg=replicate(1.,jpi)#e2v(inorth,*)
162
163;E1f scale factor
164
165e1f_reg=( e1t_eq#replicate(1.,jpj) )*(replicate(1.,jpi)#cos(!PI*gphit(inorth,*)/180.))
166
167;E2f scale factor
168
169e2f_reg=replicate(1.,jpi)#e2f(inorth,*)
170
171
172;************************************************
173;  2nd STEP : Define the mask of the outpût grid
174;************************************************
175
176
177; We now need to define the grid mask...
178
179
180mask_reg=tmask
181
182; Let's get entire coordinates of meshmask
183
184
185sph_coord=replicate(1.,3,jpi*jpj)
186
187; We put some crazy values of radius on half of the North Fold line in order
188; to avoid any problem
189
190sph_coord(2,jpj*jpi-1-jpi/2:jpj*jpi-1)=10.
191
192sph_coord(0,*)=reform(glamt,jpi*jpj)
193sph_coord(1,*)=reform(gphit,jpi*jpj)
194
195rec_grid=cv_coord(/DEGREES,FROM_SPHERE=sph_coord,/TO_RECT)
196
197xx_grid=rec_grid(0,*)
198yy_grid=rec_grid(1,*)
199zz_grid=rec_grid(2,*)
200
201
202not_found=1
203for jj=0, jpj-1 do begin
204
205lat_max_ORCA=max( gphit(*,jj)     )
206lat_max_reg =max( gphit_reg(*,jj) )
207
208
209if ((lat_max_ORCA NE lat_max_reg) AND (not_found EQ 1)) then begin
210   jpj_start=jj-2
211   not_found=0
212endif
213
214end
215
216
217
218; Loop on all elements of our geo grid
219
220for jj=jpj/2,jpj-1 do begin
221    for ji=0,jpi-1 do begin
222
223
224; Coordinates of our point
225
226    xx=glamt_reg(ji,jj)
227    yy=gphit_reg(ji,jj)
228
229    sph_coord=replicate(1.,3)
230    sph_coord(0)=xx
231    sph_coord(1)=yy
232
233; Convert them into cartesian
234
235    rec_coord=cv_coord(/DEGREES,FROM_SPHERE=sph_coord,/TO_RECT)
236
237    xx_bis=replicate(rec_coord(0),jpi*jpj)
238    yy_bis=replicate(rec_coord(1),jpi*jpj)
239    zz_bis=replicate(rec_coord(2),jpi*jpj)
240
241    modules=sqrt((xx_grid-xx_bis)*(xx_grid-xx_bis)  $
242                +(yy_grid-yy_bis)*(yy_grid-yy_bis)  $
243                +(zz_grid-zz_bis)*(zz_grid-zz_bis))
244
245
246    sort_mod=where(modules EQ min(modules))
247
248    answj=sort_mod(0)/jpi
249
250    answi=sort_mod(0)-answj*jpi
251
252
253    mask_reg(ji,jj,*)=tmask(answi,answj,*)
254
255
256    endfor
257endfor
258
259;*****************************************
260; 3rd STEP : PUT VALUES INTO ATTRIBUTES
261;*****************************************
262
263;********** T MASK*********************
264
265; Create Global Mask file for T point
266
267tmask_reg_glo=replicate(0.,jpiglo,jpjglo,jpk)
268
269tmask_reg_glo(1:jpiglo-2,0:jpjglo-2, *)=mask_reg
270
271
272; Data Duplication
273
274; North
275
276tmask_reg_glo(*,jpjglo-1,*)=tmask_reg_glo(*,jpjglo-2,*)
277;E/W
278
279tmask_reg_glo(0,*,*)=tmask_reg_glo(jpiglo-2,*,*)
280tmask_reg_glo(jpiglo-1,*,*)=tmask_reg_glo(1,*,*)
281
282; Global T mask in created
283
284; Creates U mask, V mask and F mask
285
286
287umask_reg=tmask_reg_glo(1:jpiglo-2,0:jpjglo-2,*)*tmask_reg_glo(2:jpiglo-1,0:jpjglo-2,*)
288
289vmask_reg=tmask_reg_glo(1:jpiglo-2,0:jpjglo-2,*)*tmask_reg_glo(1:jpiglo-2,1:jpjglo-1,*)
290
291fmask_reg=tmask_reg_glo(1:jpiglo-2,0:jpjglo-2,*)*tmask_reg_glo(2:jpiglo-1,0:jpjglo-2,*) $
292*tmask_reg_glo(1:jpiglo-2,1:jpjglo-1,*)*tmask_reg_glo(2:jpiglo-1,1:jpjglo-1,*)
293
294; Create Global Mask file for U,V & F points
295
296umask_reg_glo=replicate(0.,jpiglo,jpjglo,jpk)
297umask_reg_glo(1:jpiglo-2,0:jpjglo-2, *)=umask_reg
298
299vmask_reg_glo=replicate(0.,jpiglo,jpjglo,jpk)
300vmask_reg_glo(1:jpiglo-2,0:jpjglo-2, *)=vmask_reg
301
302fmask_reg_glo=replicate(0.,jpiglo,jpjglo,jpk)
303fmask_reg_glo(1:jpiglo-2,0:jpjglo-2, *)=fmask_reg
304
305; Data Duplication
306
307; North
308
309umask_reg_glo(*,jpjglo-1,*)=umask_reg_glo(*,jpjglo-2,*)
310vmask_reg_glo(*,jpjglo-1,*)=vmask_reg_glo(*,jpjglo-2,*)
311fmask_reg_glo(*,jpjglo-1,*)=fmask_reg_glo(*,jpjglo-2,*)
312
313;E/W
314
315umask_reg_glo(0,*,*)=umask_reg_glo(jpiglo-2,*,*)
316umask_reg_glo(jpiglo-1,*,*)=umask_reg_glo(1,*,*)
317
318vmask_reg_glo(0,*,*)=vmask_reg_glo(jpiglo-2,*,*)
319vmask_reg_glo(jpiglo-1,*,*)=vmask_reg_glo(1,*,*)
320
321fmask_reg_glo(0,*,*)=fmask_reg_glo(jpiglo-2,*,*)
322fmask_reg_glo(jpiglo-1,*,*)=fmask_reg_glo(1,*,*)
323
324; ***** NAV_LON & NAV_LAT **********
325
326
327; Create Global longitude and latitude
328
329glamt_reg_glo=replicate(0.,jpiglo,jpjglo)
330
331glamt_reg_glo(1:jpiglo-2,0:jpjglo-2)=glamt_reg
332
333gphit_reg_glo=replicate(0.,jpiglo,jpjglo)
334
335gphit_reg_glo(1:jpiglo-2,0:jpjglo-2)=gphit_reg
336
337nav_lon_reg_glo=replicate(0.,jpiglo,jpjglo)
338
339nav_lat_reg_glo=replicate(0.,jpiglo,jpjglo)
340
341nav_lon_reg_glo=glamt_reg_glo
342
343nav_lat_reg_glo=gphit_reg_glo
344
345; Data Duplication
346
347; E/W Transition
348
349nav_lon_reg_glo(0,*)=nav_lon_reg_glo(jpiglo-2,*)
350nav_lon_reg_glo(jpiglo-1,*)=nav_lon_reg_glo(1,*)
351
352nav_lat_reg_glo(0,*)=nav_lat_reg_glo(jpiglo-2,*)
353nav_lat_reg_glo(jpiglo-1,*)=nav_lat_reg_glo(1,*)
354
355; North
356
357nav_lon_reg_glo(*,jpjglo-1)= nav_lon_reg_glo(*,jpjglo-2)
358nav_lat_reg_glo(*,jpjglo-1)= nav_lat_reg_glo(*,jpjglo-2)
359
360; ***** GLAMT *****
361; We don't do it to save memory :  same as nav_lon
362
363; ***** GLAMU *****
364
365glamu_reg_glo=replicate(0.,jpiglo,jpjglo)
366glamu_reg_glo(1:jpiglo-2,0:jpjglo-2)=glamu_reg
367
368; North
369
370glamu_reg_glo(*,jpjglo-1)=glamu_reg_glo(*,jpjglo-2)
371
372; E/W
373
374glamu_reg_glo(0,*)=glamu_reg_glo(jpiglo-2,*)
375glamu_reg_glo(jpiglo-1,*)=glamu_reg_glo(1,*)
376
377; ***** GLAMV *****
378
379glamv_reg_glo=replicate(0.,jpiglo,jpjglo)
380glamv_reg_glo(1:jpiglo-2,0:jpjglo-2)=glamv_reg
381
382; North
383
384glamv_reg_glo(*,jpjglo-1)=glamv_reg_glo(*,jpjglo-2)
385
386; E/W
387
388glamv_reg_glo(0,*)=glamv_reg_glo(jpiglo-2,*)
389glamv_reg_glo(jpiglo-1,*)=glamv_reg_glo(1,*)
390
391
392; ***** GLAMF *****
393
394glamf_reg_glo=replicate(0.,jpiglo,jpjglo)
395glamf_reg_glo(1:jpiglo-2,0:jpjglo-2)=glamf_reg
396
397; North
398
399glamf_reg_glo(*,jpjglo-1)=glamf_reg_glo(*,jpjglo-2)
400
401; E/W
402
403glamf_reg_glo(0,*)=glamf_reg_glo(jpiglo-2,*)
404glamf_reg_glo(jpiglo-1,*)=glamf_reg_glo(1,*)
405
406
407; ***** GPHIT *****
408; We don't do it to save memory :  same as nav_lat
409
410; ***** GPHIU *****
411
412gphiu_reg_glo=replicate(0.,jpiglo,jpjglo)
413gphiu_reg_glo(1:jpiglo-2,0:jpjglo-2)=gphiu_reg
414
415; North
416
417gphiu_reg_glo(*,jpjglo-1)=gphiu_reg_glo(*,jpjglo-2)
418
419; E/W
420
421gphiu_reg_glo(0,*)=gphiu_reg_glo(jpiglo-2,*)
422gphiu_reg_glo(jpiglo-1,*)=gphiu_reg_glo(1,*)
423
424; ***** GPHIV *****
425
426gphiv_reg_glo=replicate(0.,jpiglo,jpjglo)
427gphiv_reg_glo(1:jpiglo-2,0:jpjglo-2)=gphiv_reg
428
429; North
430
431gphiv_reg_glo(*,jpjglo-1)=gphiv_reg_glo(*,jpjglo-2)
432
433; E/W
434
435gphiv_reg_glo(0,*)=gphiv_reg_glo(jpiglo-2,*)
436gphiv_reg_glo(jpiglo-1,*)=gphiv_reg_glo(1,*)
437
438
439; ***** GPHIF *****
440
441gphif_reg_glo=replicate(0.,jpiglo,jpjglo)
442gphif_reg_glo(1:jpiglo-2,0:jpjglo-2)=gphif_reg
443
444; North
445
446gphif_reg_glo(*,jpjglo-1)=gphif_reg_glo(*,jpjglo-2)
447
448; E/W
449
450gphif_reg_glo(0,*)=gphif_reg_glo(jpiglo-2,*)
451gphif_reg_glo(jpiglo-1,*)=gphif_reg_glo(1,*)
452
453; ***** FF *******
454
455omega = 2.0 * !PI / 86164.0
456
457ff_reg_glo = 2.*omega*sin(gphif_reg_glo*!PI/180.)
458
459
460; *********** E1T***************
461
462e1t_reg_glo=replicate(0.,jpiglo,jpjglo)
463e1t_reg_glo(1:jpiglo-2,0:jpjglo-2)=e1t_reg
464
465; North
466
467e1t_reg_glo(*,jpjglo-1)=e1t_reg_glo(*,jpjglo-2)
468
469; E/W
470
471e1t_reg_glo(0,*)=e1t_reg_glo(jpiglo-2,*)
472e1t_reg_glo(jpiglo-1,*)=e1t_reg_glo(1,*)
473
474
475
476; *********** E2T***************
477
478e2t_reg_glo=replicate(0.,jpiglo,jpjglo)
479e2t_reg_glo(1:jpiglo-2,0:jpjglo-2)=e2t_reg
480
481; North
482
483e2t_reg_glo(*,jpjglo-1)=e2t_reg_glo(*,jpjglo-2)
484
485; E/W
486
487e2t_reg_glo(0,*)=e2t_reg_glo(jpiglo-2,*)
488e2t_reg_glo(jpiglo-1,*)=e2t_reg_glo(1,*)
489
490; *********** E1U***************
491
492e1u_reg_glo=replicate(0.,jpiglo,jpjglo)
493e1u_reg_glo(1:jpiglo-2,0:jpjglo-2)=e1u_reg
494
495; North
496
497e1u_reg_glo(*,jpjglo-1)=e1u_reg_glo(*,jpjglo-2)
498
499; E/W
500
501e1u_reg_glo(0,*)=e1u_reg_glo(jpiglo-2,*)
502e1u_reg_glo(jpiglo-1,*)=e1u_reg_glo(1,*)
503
504
505; *********** E2U***************
506
507e2u_reg_glo=replicate(0.,jpiglo,jpjglo)
508e2u_reg_glo(1:jpiglo-2,0:jpjglo-2)=e2u_reg
509
510; North
511
512e2u_reg_glo(*,jpjglo-1)=e2u_reg_glo(*,jpjglo-2)
513
514; E/W
515
516e2u_reg_glo(0,*)=e2u_reg_glo(jpiglo-2,*)
517e2u_reg_glo(jpiglo-1,*)=e2u_reg_glo(1,*)
518
519; *********** E1V***************
520
521e1v_reg_glo=replicate(0.,jpiglo,jpjglo)
522e1v_reg_glo(1:jpiglo-2,0:jpjglo-2)=e1v_reg
523
524; North
525
526e1v_reg_glo(*,jpjglo-1)=e1v_reg_glo(*,jpjglo-2)
527
528; E/W
529
530e1v_reg_glo(0,*)=e1v_reg_glo(jpiglo-2,*)
531e1v_reg_glo(jpiglo-1,*)=e1v_reg_glo(1,*)
532
533; *********** E2V***************
534
535e2v_reg_glo=replicate(0.,jpiglo,jpjglo)
536e2v_reg_glo(1:jpiglo-2,0:jpjglo-2)=e2v_reg
537
538; North
539
540e2v_reg_glo(*,jpjglo-1)=e2v_reg_glo(*,jpjglo-2)
541
542; E/W
543
544e2v_reg_glo(0,*)=e2v_reg_glo(jpiglo-2,*)
545e2v_reg_glo(jpiglo-1,*)=e2v_reg_glo(1,*)
546
547; *********** E1F **************
548
549e1f_reg_glo=replicate(0.,jpiglo,jpjglo)
550e1f_reg_glo(1:jpiglo-2,0:jpjglo-2)=e1f_reg
551
552; North
553
554e1f_reg_glo(*,jpjglo-1)=e1f_reg_glo(*,jpjglo-2)
555
556; E/W
557
558e1f_reg_glo(0,*)=e1f_reg_glo(jpiglo-2,*)
559e1f_reg_glo(jpiglo-1,*)=e1f_reg_glo(1,*)
560
561; *********** E2F***************
562
563e2f_reg_glo=replicate(0.,jpiglo,jpjglo)
564e2f_reg_glo(1:jpiglo-2,0:jpjglo-2)=e2f_reg
565
566; North
567
568e2f_reg_glo(*,jpjglo-1)=e2f_reg_glo(*,jpjglo-2)
569
570; E/W
571
572e2f_reg_glo(0,*)=e2f_reg_glo(jpiglo-2,*)
573e2f_reg_glo(jpiglo-1,*)=e2f_reg_glo(1,*)
574
575
576
577;****************END OF ATTRIBUTES****************
578
579
580;*************************************************
581; 4th STEP : CREATES NCDF MESHMASK FILE
582;*************************************************
583
584; Creates Netcdf File to save this
585
586
587vargrid = 'T'
588
589; Name
590idout = NCDF_CREATE(outputgrid,/clobber)
591print, 'Creation du fichier Netcdf'
592NCDF_CONTROL, idout, /nofill
593; Dimension
594
595xidout = NCDF_DIMDEF(idout, 'x', jpiglo)
596yidout = NCDF_DIMDEF(idout, 'y', jpjglo)
597didout = NCDF_DIMDEF(idout, 'z', jpk)
598x_a    = NCDF_DIMDEF(idout, 'x_a',1)
599y_a    = NCDF_DIMDEF(idout, 'y_a',1)
600z_a    = NCDF_DIMDEF(idout, 'z_a',1)
601tidout = NCDF_DIMDEF(idout, 'time_counter', /unlimited)
602
603; Attributes
604
605id0  = NCDF_VARDEF(idout, 'nav_lon'     , [xidout, yidout                ],/FLOAT )
606id1  = NCDF_VARDEF(idout, 'nav_lat'     , [xidout, yidout                ],/FLOAT )
607id2  = NCDF_VARDEF(idout, 'nav_lev'     , [                didout        ],/FLOAT )
608id3  = NCDF_VARDEF(idout, 'time'        , [                        tidout],/FLOAT )
609id4  = NCDF_VARDEF(idout, 'time_steps'  , [                        tidout],/LONG  )
610id5  = NCDF_VARDEF(idout, 'glamt'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
611id6  = NCDF_VARDEF(idout, 'glamu'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
612id7  = NCDF_VARDEF(idout, 'glamv'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
613id8  = NCDF_VARDEF(idout, 'glamf'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
614id9  = NCDF_VARDEF(idout, 'gphit'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
615id10 = NCDF_VARDEF(idout, 'gphiu'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
616id11 = NCDF_VARDEF(idout, 'gphiv'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
617id12 = NCDF_VARDEF(idout, 'gphif'       , [xidout, yidout, z_a   , tidout],/DOUBLE)
618id13 = NCDF_VARDEF(idout, 'e1t'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
619id14 = NCDF_VARDEF(idout, 'e1u'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
620id15 = NCDF_VARDEF(idout, 'e1v'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
621id16 = NCDF_VARDEF(idout, 'e1f'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
622id17 = NCDF_VARDEF(idout, 'e2t'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
623id18 = NCDF_VARDEF(idout, 'e2u'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
624id19 = NCDF_VARDEF(idout, 'e2v'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
625id20 = NCDF_VARDEF(idout, 'e2f'         , [xidout, yidout, z_a   , tidout],/DOUBLE)
626id21 = NCDF_VARDEF(idout, 'tmask'       , [xidout, yidout, didout, tidout],/DOUBLE)
627id22 = NCDF_VARDEF(idout, 'umask'       , [xidout, yidout, didout, tidout],/DOUBLE)
628id23 = NCDF_VARDEF(idout, 'vmask'       , [xidout, yidout, didout, tidout],/DOUBLE)
629id24 = NCDF_VARDEF(idout, 'fmask'       , [xidout, yidout, didout, tidout],/DOUBLE)
630id25 = NCDF_VARDEF(idout, 'ff'          , [xidout, yidout, z_a   , tidout],/DOUBLE)
631id26 = NCDF_VARDEF(idout, 'gdept'       , [x_a   , y_a   , didout, tidout],/DOUBLE)
632id27 = NCDF_VARDEF(idout, 'gdepw'       , [x_a   , y_a   , didout, tidout],/DOUBLE)
633id28 = NCDF_VARDEF(idout, 'e3t'         , [x_a   , y_a   , didout, tidout],/DOUBLE)
634id29 = NCDF_VARDEF(idout, 'e3w'         , [x_a   , y_a   , didout, tidout],/DOUBLE)
635
636; Variable 0
637NCDF_ATTPUT, idout, id0, 'units', 'degrees_east'
638NCDF_ATTPUT, idout, id0, 'long_name', 'Longitude'
639NCDF_ATTPUT, idout, id0, 'nav_model', 'Default grid'
640; Variable 1
641NCDF_ATTPUT, idout, id1, 'units', 'degrees_north'
642NCDF_ATTPUT, idout, id1, 'long_name', 'Latitude'
643NCDF_ATTPUT, idout, id1, 'nav_model', 'Default grid'
644; Variable 2
645NCDF_ATTPUT, idout, id2, 'units','meters'
646NCDF_ATTPUT, idout, id2, 'long_name','Depth'
647NCDF_ATTPUT, idout, id2, 'nav_model','Default grid'
648; Variable 3
649NCDF_ATTPUT, idout, id3, 'units', 'seconds since 0001-01-15 12:00:00 '
650NCDF_ATTPUT, idout, id3, 'calendar','noleap'
651NCDF_ATTPUT, idout, id3, 'title', 'Time'
652NCDF_ATTPUT, idout, id3, 'long_name', 'Time axis'
653NCDF_ATTPUT, idout, id3, 'time_origin','0000-DEC-15 00:00:00'
654
655
656NCDF_CONTROL, idout, /ENDEF
657
658
659
660; Writting
661
662
663NCDF_VARPUT, idout, id0 , nav_lon_reg_glo
664NCDF_VARPUT, idout, id1 , nav_lat_reg_glo
665NCDF_VARPUT, idout, id2 , gdept
666NCDF_VARPUT, idout, id3 , 1.
667NCDF_VARPUT, idout, id4 , 1.
668NCDF_VARPUT, idout, id5 , nav_lon_reg_glo
669
670NCDF_VARPUT, idout, id6 , glamu_reg_glo
671NCDF_VARPUT, idout, id7 , glamv_reg_glo
672NCDF_VARPUT, idout, id8 , glamf_reg_glo
673NCDF_VARPUT, idout, id9 , nav_lat_reg_glo
674NCDF_VARPUT, idout, id10, gphiu_reg_glo
675NCDF_VARPUT, idout, id11, gphiv_reg_glo
676NCDF_VARPUT, idout, id12, gphif_reg_glo
677NCDF_VARPUT, idout, id13, e1t_reg_glo
678NCDF_VARPUT, idout, id14, e1u_reg_glo
679NCDF_VARPUT, idout, id15, e1v_reg_glo
680NCDF_VARPUT, idout, id16, e1f_reg_glo
681NCDF_VARPUT, idout, id17, e2t_reg_glo
682NCDF_VARPUT, idout, id18, e2u_reg_glo
683NCDF_VARPUT, idout, id19, e2v_reg_glo
684NCDF_VARPUT, idout, id20, e2f_reg_glo
685NCDF_VARPUT, idout, id21, tmask_reg_glo
686NCDF_VARPUT, idout, id22, umask_reg_glo
687NCDF_VARPUT, idout, id23, vmask_reg_glo
688NCDF_VARPUT, idout, id24, fmask_reg_glo
689NCDF_VARPUT, idout, id25, ff_reg_glo
690
691dummy=replicate(0.,1,1,jpk,1)
692dummy(0,0,*,0)=gdept
693
694NCDF_VARPUT, idout, id26, dummy
695
696 dummy=replicate(0.,1,1,jpk,1)
697 dummy(0,0,*,0)=gdepw
698
699NCDF_VARPUT, idout, id27, dummy
700
701 dummy=replicate(0.,1,1,jpk,1)
702 dummy(0,0,*,0)=e3t
703
704NCDF_VARPUT, idout, id28, dummy
705
706 dummy=replicate(0.,1,1,jpk,1)
707 dummy(0,0,*,0)=e3w
708
709NCDF_VARPUT, idout, id29, dummy
710
711NCDF_CLOSE, idout
712
713
714
715end
716
Note: See TracBrowser for help on using the repository browser.