Dear Dr. Fabien

     Thank you very much for your reply.  with setting the ichebdy=0 it works.
What is the differences between ichdustemd = 1 and ichdustemd = 2 ? If I'd like to generate dust boundary
conditions, how I can do it and where the global dust emission data?

thank you
Rahul




From: Solmon Fabien <fsolmon@ictp.it>
To: R Climate <rclimate@yahoo.com>
Cc: "regcnet@lists.ictp.it" <regcnet@lists.ictp.it>
Sent: Saturday, January 12, 2013 10:35:16 AM
Subject: Re: [RegCNET] Dust Run



Hello !

I think you need to rerun terrain with the option ltexture = true. After
that you can check that you have the texture fields ( needed by dust
emission) in your DOMAIN000.nc file.
Also you need to perform the chbc step to generate dust boundary
conditions if you want them, or set ichebdy=0 if you want to relax toward
zero concentration at the boundary ( in which case the chbc step is not
necessary).

Finally the rdtsemfac might also be adjusted.
Usually I use rdstemfac=1 (default)  when I use the new option
ichdustemd = 2 ( which I think is better at least over the sahara). With
the option ichdustemd =1 , the emission more intense and I remember using
a rdstemfac = 0.7 to get more realistic emissions. But again all of this
domain dependant ..

Fabien


> Hello Dr. Giulian and Regcnet,
>      I'm running a dust simulation with the RegCM4.3.5.1 and I put the
> following selection in the name-list:
>
>  ltexture = .true.,
>
>  chemsimtype = 'DUST    ',
>  ichsolver = 1,
>  ichsursrc = 1,
>  ichdrdepo = 1,
>  ichebdy = 1,
>  ichcumtra = 1,
>  ichremlsc = 1,
>  ichremcvc = 1,
>  ichdustemd = 1,
>  ichdiag = 0,
>  idirect   = 1,
>  rdstemfac = 1.0,
>
>
> Then, when I run the simulation I got the following:
>
>  Constant hor. diff. coef. =  0.360000E+05 m^2 s-1
>   Maximumt hor. diff. coef. =  0.375000E+06 m^2 s-1
>  Calculate solar declination angle at    2009010100
>  JDay         0.00 solar declination angle = -23.05862917 degrees
>                    solar TSI irradiance    =    1367.0000 W/m^2
>  READY IC DATA for 2009-01-01 00:00:00 UTC        
>  SEARCH BC data for    2009010106
>  READY  BC from        2009010100  to    2009010106
>   # Ozone profiles #
>     0.496E-07  0.496E-07  0.496E-07  0.496E-07  0.496E-07  0.519E-07 
> 0.527E-07
>     0.574E-07  0.623E-07  0.710E-07  0.846E-07  0.106E-06  0.148E-06 
> 0.223E-06
>     0.223E-06  0.409E-06  0.793E-06  0.141E-05  0.317E-05   0.00   
>  Calling inidust
>  No such file or
> directory                                                      
> --------------------------------------------------------------------------
> MPI_ABORT was invoked on rank 7 in communicator MPI_COMM_WORLD
> with errorcode 1.
>
> NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
> You may or may not see output from other processes, depending on
> exactly when Open MPI kills them.
>
>
> Any help is welcome
>
>
> Thanks
> Rahul
> Manger of Climate Division_______________________________________________
> RegCNET mailing list
> RegCNET@lists.ictp.it
> https://lists.ictp.it/mailman/listinfo.cgi/regcnet


--
+39 040 2240206 (office) | Earth System Physics
+39 040 2240449 (fax)    | ICTP (UNESCO-IAEA)
esp.ictp.it              | Strada Costiera 11
                        | 34014 Trieste, Italy