[5246] | 1 | SUBROUTINE read_vent(debutphy, step, nbjour, u10m_ec, v10m_ec) |
---|
| 2 | USE dimphy |
---|
| 3 | USE mod_grid_phy_lmdz |
---|
| 4 | USE mod_phys_lmdz_para |
---|
[5270] | 5 | USE netcdf, ONLY: nf90_get_var, nf90_open, nf90_inq_varid, nf90_nowrite |
---|
[5246] | 6 | ! USE write_field_phy |
---|
[5271] | 7 | USE dimensions_mod, ONLY: iim, jjm, llm, ndm |
---|
[5272] | 8 | USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & |
---|
| 9 | ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm |
---|
[5271] | 10 | IMPLICIT NONE |
---|
| 11 | |
---|
[5246] | 12 | ! INCLUDE "dimphy.h" |
---|
[5272] | 13 | |
---|
[5246] | 14 | ! |
---|
| 15 | INTEGER :: step, nbjour |
---|
| 16 | LOGICAL :: debutphy |
---|
| 17 | real :: u10m_ec(klon), v10m_ec(klon) |
---|
| 18 | real :: u10m_ec_glo(klon_glo), v10m_ec_glo(klon_glo) |
---|
| 19 | ! |
---|
| 20 | ! real u10m_nc(iip1,jjp1) !, v10m_nc(iip1,jjm) ! dim 97x72 |
---|
| 21 | ! real v10m_nc(iip1,jjp1) ! dim 97x73 |
---|
| 22 | real :: u10m_nc_glo(nbp_lon+1,nbp_lat) !, v10m_nc(iip1,jjm) ! dim 97x72 |
---|
| 23 | real :: v10m_nc_glo(nbp_lon+1,nbp_lat) ! dim 97x73 |
---|
[5270] | 24 | integer :: ncidu1, varidu1, ncidv1, varidv1, rcode |
---|
[5246] | 25 | save ncidu1, varidu1, ncidv1, varidv1 |
---|
[2630] | 26 | !$OMP THREADPRIVATE(ncidu1, varidu1, ncidv1, varidv1) |
---|
[5246] | 27 | integer :: start(4),count(4), status |
---|
| 28 | integer :: i, j, ig |
---|
[2630] | 29 | |
---|
| 30 | |
---|
[5246] | 31 | ! |
---|
[2630] | 32 | !$OMP MASTER |
---|
[5246] | 33 | IF (is_mpi_root .AND. is_omp_root) THEN |
---|
| 34 | if (debutphy) then |
---|
| 35 | ! |
---|
[5270] | 36 | ncidu1=nf90_open('u10m.nc',nf90_nowrite,rcode) |
---|
| 37 | varidu1=nf90_inq_varid(ncidu1,'U10M',rcode) |
---|
| 38 | ncidv1=nf90_open('v10m.nc',nf90_nowrite,rcode) |
---|
| 39 | varidv1=nf90_inq_varid(ncidv1,'V10M',rcode) |
---|
[5246] | 40 | ! |
---|
| 41 | endif |
---|
| 42 | ! |
---|
| 43 | start(1)=1 |
---|
| 44 | start(2)=1 |
---|
| 45 | start(4)=0 |
---|
[2630] | 46 | |
---|
[5246] | 47 | ! count(1)=iip1 |
---|
| 48 | count(1)=nbp_lon+1 |
---|
| 49 | ! count(2)=jjp1 |
---|
| 50 | count(2)=nbp_lat |
---|
| 51 | count(3)=1 |
---|
| 52 | count(4)=0 |
---|
| 53 | ! |
---|
| 54 | start(3)=step |
---|
| 55 | ! |
---|
[5249] | 56 | status = nf90_get_var(ncidu1, varidu1, u10m_nc_glo, start, count) |
---|
| 57 | |
---|
[5246] | 58 | ! print *,status |
---|
| 59 | ! |
---|
[5249] | 60 | status = nf90_get_var(ncidv1, varidv1, v10m_nc_glo, start, count) |
---|
| 61 | |
---|
[5246] | 62 | ! |
---|
[2630] | 63 | |
---|
[5246] | 64 | ! print *,'beforebidcor u10m_nc', u10m_nc(1,jjp1) |
---|
| 65 | ! print *,'beforebidcor v10m_nc', v10m_nc(1,jjp1) |
---|
[2630] | 66 | |
---|
[5246] | 67 | ! print *,status |
---|
| 68 | ! call correctbid(iim,jjp1,u10m_nc) |
---|
| 69 | ! call correctbid(iim,jjp1,v10m_nc) |
---|
| 70 | call correctbid(nbp_lon,nbp_lat,u10m_nc_glo) |
---|
| 71 | call correctbid(nbp_lon,nbp_lat,v10m_nc_glo) |
---|
[2630] | 72 | |
---|
[5246] | 73 | ! print *,'afterbidcor u10m_nc', u10m_nc(1,jjp1) |
---|
| 74 | ! print *,'afterbidcor v10m_nc', v10m_nc(1,jjp1) |
---|
| 75 | ! |
---|
| 76 | !--upside down + physical grid |
---|
| 77 | ! |
---|
| 78 | ! u10m_ec(1)=u10m_nc(1,jjp1) |
---|
| 79 | ! v10m_ec(1)=v10m_nc(1,jjp1) |
---|
| 80 | u10m_ec_glo(1)=u10m_nc_glo(1,nbp_lat) |
---|
| 81 | v10m_ec_glo(1)=v10m_nc_glo(1,nbp_lat) |
---|
| 82 | ig=2 |
---|
| 83 | ! DO j=2,jjm |
---|
| 84 | ! DO i = 1, iim |
---|
| 85 | DO j=2,nbp_lat-1 |
---|
| 86 | DO i = 1, nbp_lon |
---|
| 87 | ! u10m_ec(ig)=u10m_nc(i,jjp1+1-j) |
---|
| 88 | ! v10m_ec(ig)=v10m_nc(i,jjp1+1-j) |
---|
| 89 | u10m_ec_glo(ig)=u10m_nc_glo(i,nbp_lat+1-j) |
---|
| 90 | v10m_ec_glo(ig)=v10m_nc_glo(i,nbp_lat+1-j) |
---|
| 91 | ig=ig+1 |
---|
| 92 | ! print *,u10m_ec(ig) ,v10m_ec(ig) |
---|
| 93 | ENDDO |
---|
| 94 | ENDDO |
---|
| 95 | u10m_ec_glo(ig)=u10m_nc_glo(1,1) |
---|
| 96 | v10m_ec_glo(ig)=v10m_nc_glo(1,1) |
---|
[2630] | 97 | |
---|
| 98 | |
---|
[5246] | 99 | ! end if master |
---|
| 100 | ENDIF |
---|
[2630] | 101 | !$OMP END MASTER |
---|
| 102 | !$OMP BARRIER |
---|
[5246] | 103 | CALL scatter(u10m_ec_glo,u10m_ec) |
---|
| 104 | CALL scatter(v10m_ec_glo,v10m_ec) |
---|
[2630] | 105 | |
---|
[5246] | 106 | ! print *,'JE tamagno viento ig= ', ig |
---|
| 107 | ! print *,'READ_VENT U = ',SUM(u10m_ec),MINVAL(u10m_ec), |
---|
| 108 | ! . MAXVAL(u10m_ec) |
---|
| 109 | ! print *,'READ_VENT V = ',SUM(v10m_ec),MINVAL(v10m_ec), |
---|
| 110 | ! . MAXVAL(v10m_ec) |
---|
| 111 | ! print *,'u v 1 ', u10m_ec(1),v10m_ec(1) |
---|
| 112 | ! print *,'u v klon ', u10m_ec(klon),v10m_ec(klon) |
---|
| 113 | RETURN |
---|
| 114 | END SUBROUTINE read_vent |
---|
[2630] | 115 | |
---|
[5246] | 116 | ! added by JE from the nh SPLA, dyn3d/read_reanalyse.F which is not available any more |
---|
| 117 | subroutine correctbid(iim,nl,x) |
---|
| 118 | integer :: iim,nl |
---|
| 119 | real :: x(iim+1,nl) |
---|
| 120 | integer :: i,l |
---|
| 121 | real :: zz |
---|
[2630] | 122 | |
---|
[5246] | 123 | do l=1,nl |
---|
| 124 | do i=2,iim-1 |
---|
| 125 | if(abs(x(i,l)).gt.1.e10) then |
---|
| 126 | zz=0.5*(x(i-1,l)+x(i+1,l)) |
---|
| 127 | ! print*,'correction ',i,l,x(i,l),zz |
---|
| 128 | x(i,l)=zz |
---|
| 129 | endif |
---|
| 130 | enddo |
---|
| 131 | enddo |
---|
[2630] | 132 | |
---|
[5246] | 133 | return |
---|
| 134 | end subroutine correctbid |
---|
[2630] | 135 | |
---|
| 136 | |
---|
| 137 | |
---|