[RegCNET] 回复： a question in running RegCM3

[子 夫]

Sorry. I forgot to add the domain.param.




----- 原始邮件 ----
发件人： XUNQIANG BI <bixq at ictp.it>
收件人： 子 夫 <pine_guest at yahoo.com.cn>
抄送： regcnet at lists.ictp.it
已发送： 2008/6/13(周五), 下午4:24:26
主题： Re: [RegCNET] a question in running RegCM3


Hi, Zhao:

You'd better let us know which UNIX computer you used, as well as the
domain.param file ?

I guess that it's not a problem of 10 minute resolution datasets,
(of course, you can easily check this point by using 5 minute instead),
at least it's not a problem for LINUX users.

By the way, Ashraf said that he think the dust package (under RegCM/Main)
need to be improved. And he will tell the details later this month.

So you'd better wait for the updated version for dust simulation.

Best Regards,
Bi Xunqiang

On Thu, 12 Jun 2008, [gb2312] вс ╥Р wrote:

> Hi. 
б═б═б═б═ I am a new one to use RegCM3 model. When I ran it to simulate the dust transportation, I found some blocks prevented me from the fluent run. 
б═б═б═ The problem is associated with readingб═SOILCAT.10 in Preproc/Terrain/terrain.f. I have setб═ntypec=10 in advance. I printed the following sentences intentionally to know what reason caused the mistake. Here are the print sentence and signaled mistakes:
Print sentences in terrain.f:
б═б═б═б═б═б═б═б═б═б═б═ б═б═б═ print*,'AERTYP2',AERTYP
б═б═б═ б═б═б═ б═б═б═б═б═б═б═ print*,'б═б═ LANDUSE READ1IN', ilat1,ilat2,irec,'==',ch_tex(1,1)
б═б═б═ б═б═б═б═б═б═б═ б═б═б═ print*,'б═б═ LANDUSE READ2IN', nxmax,ntex,ibyte,ntypec ; !!pause
б═б═б═б═б═б═б═ б═б═б═ б═б═б═ if(AERTYP(7:7).eq.'1')
б═б═б═ б═б═б═ б═б═б═б═ &б═б═б═ б═б═б═ read(45,rec=irec) ((ch_tex(j,k),k=1,ntex),j=1,nxmax)
б═б═б═ б═б═б═б═б═б═б═ б═б═б═ if(AERTYP(7:7).eq.'1') then
б═б═б═б═б═б═б═б═б═б═б═б═б═б═б═  print*,'б═б═ LANDUSE READ3IN', ilat1,ilat2,irec,ch_tex(1,3)
б═б═б═б═б═б═б═ б═б═б═ б═б═б═ pause
б═б═б═ б═б═б═ б═б═б═б═б═б═б═ endif
б═б═б═б═б═б═б═ c........process slice and store in array stores
Prompted mistakes:
б═б═б═ б═б═б═ б═б═б═ ELEVATION STD DEV READ INб═б═б═б═б═б═б═б═ 217б═б═б═б═б═б═б═б═ 437б═б═б═б═б═б═б═б═ 271
б═AERTYP2AER00D1
б═б═б═ LANDUSE READ1INб═б═б═б═б═б═б═б═ 217б═б═б═б═б═б═б═б═ 437б═б═б═б═б═б═б═б═ 271 ==
б═б═б═ LANDUSE READ2INб═б═б═б═б═б═б═ 2160б═б═б═б═б═б═б═б═б═ 17б═б═б═б═б═б═б═б═б═б═ 4б═б═б═б═б═б═б═б═б═ 10
forrtl: severe (36): attempt to access non-existent record, unit 45, file /lun02/kys/MODEL/RegCMv3/PreProc/Terrain/../DATA/SURFACE/SOILCAT.10
Imageб═б═б═б═б═б═б═б═б═б═б═б═б═ PCб═б═б═б═б═б═б═б═б═б═б═б═б═б═б═ Routineб═б═б═б═б═б═б═б═б═б═б═ Lineб═б═б═б═б═б═б═ Sourceб═б═б═б═б═б═б═б═б═б═б═б═ 
terrainб═б═б═б═б═б═б═б═б═б═б═ 400000000021F2D0б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 400000000021CDE0б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 40000000001921D0б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 40000000001159D0б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 40000000001149A0б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 400000000013CB10б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 400000000013A510б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 400000000008DF70б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 400000000000C240б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 40000000000036D0б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
libc.so.6.1б═б═б═б═б═б═б═ 2000000000427C20б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
terrainб═б═б═б═б═б═б═б═б═б═б═ 4000000000003500б═ Unknownб═б═б═б═б═б═б═б═б═б═б═б═б═б═ Unknownб═ Unknown
б═б═б═б═б═б═б═ I think that the SOILCAT.10б═ may have some problems, such as not including global but regional area. Here isб═size of the data file given by commanding 'ls -lt SOILCAT.10'б═under UNIX:
-rw-r--r-- 1 kys usersб═ 39657600 2006-08-13 04:12 SOILCAT.10
б═б═б═б═б═б═б═ I have tyied to download another SOILCAT.10, but the same result. I donot know why. Can you do me a favor.
б═б═б═б═б═б═б═ Thanks
б═б═б═б═б═б═б═ Sincerely
б═б═б═б═б═б═б═б═Jianhua Zhao
б═б═б═б═б═б═б═б═ Jun 12 2008


       ___________________________________________________________
  И⌡┘Х≥▌И┌╝Г╝╠О╪▄Ф┌╗Г └Г╩┬Г■÷И┌╝Г╝╠О╪│ 
http://cn.mail.yahoo.com/ 

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   Dr. Xunqiang Bi         email:bixq at ictp.it
   Earth System Physics Group
   The Abdus Salam ICTP
   Strada Costiera, 11
   P.O. BOX 586, 34100 Trieste, ITALY
   Tel: +39-040-2240302  Fax: +39-040-2240449
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



      ___________________________________________________________ 
 雅虎邮箱，您的终生邮箱！ 
http://cn.mail.yahoo.com/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20080613/36d3cd32/attachment-0002.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: domain.param
Type: application/octet-stream
Size: 8213 bytes
Desc: not available
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20080613/36d3cd32/attachment-0002.obj>