----- 原始邮件 ----
发件人： XUNQIANG BI <bixq@ictp.it>
收件人： 子 夫 <pine_guest@yahoo.com.cn>
抄送： regcnet@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@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
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~