Changing the terrain hight in idealized WRF model (em-les)

[skateboy366]

Hi !

Is there any ways to modify the terrain height in the em_les mode ?
Actually, I want to run an idealized WRF-LES case for a real terrain type.
I have found some approaches hadling the terrain slope in idealzed cases, but it looks not enough for a real terrain type .

Thanks.
Zihan Zhao

 

[Ming Chen]

Please go to the code "dyn_em/module_initialize_ideal.F", and starting from around line 567, you can find the following piece of code:
567 CASE (les)
568
569 CALL nl_set_mminlu(1, ' ')
570 CALL nl_set_iswater(1,0)
571 CALL nl_set_cen_lat(1,40.)
572 CALL nl_set_cen_lon(1,-105.)
573 CALL nl_set_truelat1(1,0.)
574 CALL nl_set_truelat2(1,0.)
575 CALL nl_set_moad_cen_lat (1,0.)
576 CALL nl_set_stand_lon (1,0.)
577 CALL nl_set_pole_lon (1,0.)
578 CALL nl_set_pole_lat (1,90.)
579 CALL nl_set_map_proj(1,0)
580
581
582 ! here we initialize data that currently is not initialized
583 ! in the input data
584
585 DO j = jts, jte
586 DO i = its, ite
587 grid%msftx(i,j) = 1.
588 grid%msfty(i,j) = 1.
589 grid%msfux(i,j) = 1.
590 grid%msfuy(i,j) = 1.
591 grid%msfvx(i,j) = 1.
592 grid%msfvx_inv(i,j)= 1.
593 grid%msfvy(i,j) = 1.
594 grid%sina(i,j) = 0.
595 grid%cosa(i,j) = 1.

.....

This is where you can add terrain for LES case.

 

[skateboy366]

OK！ Thanks a lot for your reply.

 

[Sauljipe]

Hi Zihan Zhao!

Did you solve the problem?
I am trying to do the same with no successful results yet.

Hope you can answer!

My best regards