[MITgcm-support] OB*file interpolation
Riema Rachmayani
imoth_22 at yahoo.com
Fri Mar 16 05:04:58 EDT 2007
dear MITgcm,
sorry i send message again in the same subject...there was an error in my yahoo
Martin, i'm trying your suggestion about OB file in your last message :
Martin :
exp4 gives an example of how to specify CONSTANT (in time) open
boundary values from files:
useOBCSprescribe = .TRUE., (enables reading from files in general)
OBNuFile = 'OBmeridU.bin', (example of a file)
[...]
these files were generated with exp4/input/gendata.m. The script tell
you that a filed for the northern (southern) boundary should have the
dimension (nx,nr) and for eastern (western), (ny,nr).
in data.obcs of exp4 :
# OBNuFile = 'OBmeridU.bin',
# OBSuFile = 'OBmeridU.bin',
# OBWuFile = 'OBzonalU.bin',
# OBEuFile = 'OBzonalU.bin',
# OBWsFile = 'OBzonalS.bin',
i checked dimension of OBmeridU ==> (80,8) for 2 iteration time independent ??
i checked dimension of OBzonalU ==> (42,8) for 2 iteration time independent ??
Martin :
If you want to specify time dependent boundary values you can do that
too. For example, you have hourly data you have create a field
(nx,nr,nt), where nt is the number of "hours" you have. and then you
use the parameters
periodicExternalForcing=.TRUE.,
externForcingPeriod=3600.,
in the file "data", in namelist PARM03. If your data is the same
after say 12 hour (ie., that you want to use the fields repeatedly)
then add this parameter in data
externForcingCycle=43200.,
in data file my case :
# Time stepping parameters
&PARM03
niter0=0,
nTimeSteps=1641600,
deltaT=1.5,
abEps=0.1,
pChkptFreq=0.0,
chkptFreq=0.0,
dumpFreq=11160.,
monitorFreq=2500.,
# for time dependent open boundary conditions comment
# in the following three lines
periodicExternalForcing=.TRUE.,
externForcingPeriod=3600.,
externForcingCycle=1641600.,
&
is that correct???
i have value of velocity at one point for 19 days (attachment file =U_AOK.IN), that file just for an example, i have other file for other point. I mean i pick out a point to other point, for example at north open boundary, i picked out point A (tip point at west)==>(Nx,Ny,Nr)==>(1,1,1), and then i picked out point B (tip point at east))==>(Nx,Ny,Nr)==>(111,1,1).
i still confused, in my attachment file U_AOK.IN, i only have 1 value for 1 point every 1 hour for 19 days like i told you before, mean while in exp4 case, OBmeridU was already in (nx,nr,nt)==>(80,8,2)...so in data.obcs wrote OBNuFile = 'OBmeridU.bin, what about my case? i only have 1 point for 19 days
what i want to do is to interpolate value in U_AOK.IN in spatial (inner of point A and point B) and in time just like you said...(Martin : Forgot to say: the model does the interpolation in time to each timestep, but that's "only" a linear interpolation)
so what i have to do with my U_AOK.IN file to get what i want?? interpolate in time first (how the code configuration in data file or data.obcs or other routine) and then interpolate in spatial (how the code configuration in data file or data.obcs or other routine) or vice versa???
thank you martin, i hope my questions is clear enough for you...
i'll wait for your reply...i really need this...
thank you...
best regards,
rima
Send instant messages to your online friends http://uk.messenger.yahoo.com
-3.0460
-3.7236
-3.6199
-2.7506
-1.3155
0.3511
1.8643
2.8838
3.1977
2.7720
1.7554
0.4382
-0.8247
-1.7053
-1.9864
-1.6151
-0.7142
0.4523
1.5491
2.2586
2.3634
1.8039
0.6941
-0.7062
-2.0577
-3.0264
-3.3686
-2.9920
-1.9782
-0.5621
0.9274
2.1490
2.8341
2.8504
2.2329
1.1700
-0.0465
-1.0984
-1.7221
-1.7747
-1.2680
-0.3625
0.6767
1.5509
2.0056
1.8973
1.2306
0.1570
-1.0632
-2.1279
-2.7689
-2.8200
-2.2584
-1.2101
0.0827
1.3219
2.2273
2.6054
2.3948
&nbs p; 1.6783
0.6581
-0.3973
-1.2222
-1.6210
-1.5155
-0.9617
-0.1323
0.7308
1.3801
1.6274
1.3919
0.7212
-0.2224
-1.2027
-1.9707
-2.3273
-2.1733
-1.5344
-0.5536
0.5443
1.5088
2.1252
2.2661
1.9205
1.1934
0.2769
-0.5995
-1.2271
-1.4686
-1.2905
-0.7680
-0.0625
0.6219
1.0922
1.2156
0.9534
0.3697
-0.3862
-1.1182
-1.6336
-1.7920
-1.5402
-0.9257
-0.0839
0.7945
1.5114
1.9081
1.9035
1.5116
0.8358
0.0420
-0.6836
-1.1796
-1.3476
-1.1740
-0.7299
-0.1506
0.4010
0.7754
0.8748
&n bsp; 0.6765
0.2370
-0.3234
-0.8517
-1.2016
-1.2702
-1.0237
-0.5064
0.1706
0.8558
1.3939
1.6628
1.6022
1.2269
0.6226
-0.0755
-0.7153
-1.1632
-1.3348
-1.2131
-0.8499
-0.3509
0.1514
0.5308
0.6981
0.6221
0.3357
-0.0733
-0.4853
-0.7791
-0.8627
-0.6965
-0.3040
0.2328
0.7928
1.2442
1.4757
1.4233
1.0862
0.5275
-0.1396
-0.7774
-1.2555
-1.4797
-1.4138
-1.0866
-0.5843
-0.0295
0.4485
0.7446
0.8022
0.6266
0.2825
-0.1227
-0.4664
-0.6418
-0.5856
-0.2953
0.1685
0.6954
1.1514
1.4113
1.3904
1.0669
0.4896
-0.2311
-0.9470
-1.5065
-1.7897
-1.7369
-1.3632
-0.7559
-0.0538
0.5858
1.0245
1.1758
1.0265
0.6386
0.1338
-0.3377
-0.6360
-0.6684
-0.4138
0.0690
0.6550
1.1831
1.4965
1.4840
1.1106
0.4298
-0.4257
-1.2733
-1.9235
-2.2256
-2.1055
-1.5849
-0.7779
0.1354
0.9520
1.4951
1.6571
1.4262
0.8887
0.2064
-0.4245
-0.8249
-0.8794
-0.5679
0.0275
0.7417
1.3687
1.7142
1.6478
1.1398
0.2734
-0.7721
-1.7625
-2.4655
-2.7083
-2.4220-1.6609
-0.5909
0.5488
1.5042
2.0686
2.1335
1.7153
0.9500
0.0588
-0.7099
-1.1449
-1.1274
-0.6638
0.1151
0.9862
1.6921
2.0082
1.8035
1.0789
-0.0286
-1.2758
-2.3720
-3.0503
-3.1357
-2.5906
-1.5268
-0.1803
1.1457
2.1548
2.6308
2.4911
1.8064
0.7820
-0.2956
-1.1330
-1.5061
-1.3191
-0.6300
0.3628
1.3719
2.0956
2.2978
1.8725
0.8758
-0.4840
-1.8948
-3.0177
-3.5717
-3.4062
-2.5405
-1.1620
0.4182
1.8415
2.7902
3.0663
2.6399
1.6558
0.3950
-0.7980
-1.6046
-1.8148
-1.3818
-0.4350
0.7534
1.8387
2.4945
2.4994
1.7976
0.5182
-1.0534
-2.5453
-3.5913
-3.9240
-3.4434
-2.2446
-0.5964
1.1242
2.5249
3.2955
3.2842
2.5346
1.2726
-0.1524
-1.3584
-2.0294
-1.9967
-1.2831
-0.0962
1.2248
2.2981
2.7995
2.5487
1.5595
0.0401
-1.6566
-3.1200
-3.9854
-4.0268
-3.2151
-1.7278
0.0940
1.8296
3.0830
3.5808
3.2396
2.1849
0.7168
-0.7679
-1.8782
-2.3294
-2.0160
-1.0384
0.3240
1.6857
2.6556
2.9398
2.4198
; 1.1857
-0.4838
-2.1895
-3.5124
-4.1189
-3.8460
-2.7437
-1.0636
0.8031
2.4227
3.4264
3.6017
2.9436
1.6530
0.0841
-1.3483
-2.2743
-2.4638
-1.8845
-0.7081
0.7342
2.0392
2.8361
2.8837
2.1336
0.7450
-0.9541
-2.5492
-3.6436
-3.9585
-3.4042
-2.1028
-0.3585
1.4185
2.8135
3.5108
3.3703
2.4607
1.0389
-0.5167
-1.8042
-2.4996
-2.4390
-1.6582
-0.3815
1.0376
2.2086
2.8048
2.6482
1.7552
0.3333
-1.2708
-2.6592
-3.4827
-3.5295
-2.7780
-1.4021
0.2720
1.8475
2.9559
3. 3477
2.9524
1.8936
0.4556
-0.9910
-2.0826
-2.5545
-2.3060
-1.4228
-0.1529
1.1585
2.1589
2.5782
2.2977
1.3791
0.0491
-1.3595
-2.4912
-3.0584
-2.9126
-2.0818
-0.7625
0.7298
2.0394
2.8582
3.0015
2.4512
1.3581
0.0016
-1.2809
-2.1792
-2.4861
-2.1468
-1.2697
-0.0958
1.0641
1.9089
2.2229
1.9288
Send instant messages to your online friends http://uk.messenger.yahoo.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20070316/90f9fbce/attachment.htm>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: U_AOK.IN
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20070316/90f9fbce/attachment.el>
More information about the MITgcm-support
mailing list