Changeset 13781
- Timestamp:
- 2020-11-11T18:40:04+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice
- Files:
-
- 1 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/OCE/IOM/iom.F90
r13557 r13781 1914 1914 IF( iom_use(cdname) ) THEN 1915 1915 #if defined key_iomput 1916 IF( SIZE(pfield2d, dim=1) == jpi .AND. SIZE(pfield2d, dim=2) == jpj ) THEN 1917 CALL xios_send_field( cdname, pfield2d(Nis0:Nie0, Njs0:Nje0) ) ! this extraction will create a copy of pfield2d 1918 ELSE 1919 CALL xios_send_field( cdname, pfield2d ) 1920 ENDIF 1916 CALL xios_send_field( cdname, pfield2d ) 1921 1917 #else 1922 1918 WRITE(numout,*) pfield2d ! iom_use(cdname) = .F. -> useless test to avoid compilation warnings … … 1930 1926 IF( iom_use(cdname) ) THEN 1931 1927 #if defined key_iomput 1932 IF( SIZE(pfield2d, dim=1) == jpi .AND. SIZE(pfield2d, dim=2) == jpj ) THEN 1933 CALL xios_send_field( cdname, pfield2d(Nis0:Nie0, Njs0:Nje0) ) ! this extraction will create a copy of pfield2d 1934 ELSE 1935 CALL xios_send_field( cdname, pfield2d ) 1936 ENDIF 1928 CALL xios_send_field( cdname, pfield2d ) 1937 1929 #else 1938 1930 WRITE(numout,*) pfield2d ! iom_use(cdname) = .F. -> useless test to avoid compilation warnings … … 1946 1938 IF( iom_use(cdname) ) THEN 1947 1939 #if defined key_iomput 1948 IF( SIZE(pfield3d, dim=1) == jpi .AND. SIZE(pfield3d, dim=2) == jpj ) THEN 1949 CALL xios_send_field( cdname, pfield3d(Nis0:Nie0, Njs0:Nje0,:) ) ! this extraction will create a copy of pfield3d 1950 ELSE 1951 CALL xios_send_field( cdname, pfield3d ) 1952 ENDIF 1940 CALL xios_send_field( cdname, pfield3d ) 1953 1941 #else 1954 1942 WRITE(numout,*) pfield3d ! iom_use(cdname) = .F. -> useless test to avoid compilation warnings … … 1962 1950 IF( iom_use(cdname) ) THEN 1963 1951 #if defined key_iomput 1964 IF( SIZE(pfield3d, dim=1) == jpi .AND. SIZE(pfield3d, dim=2) == jpj ) THEN 1965 CALL xios_send_field( cdname, pfield3d(Nis0:Nie0, Njs0:Nje0,:) ) ! this extraction will create a copy of pfield3d 1966 ELSE 1967 CALL xios_send_field( cdname, pfield3d ) 1968 ENDIF 1952 CALL xios_send_field( cdname, pfield3d ) 1969 1953 #else 1970 1954 WRITE(numout,*) pfield3d ! iom_use(cdname) = .F. -> useless test to avoid compilation warnings … … 1978 1962 IF( iom_use(cdname) ) THEN 1979 1963 #if defined key_iomput 1980 IF( SIZE(pfield4d, dim=1) == jpi .AND. SIZE(pfield4d, dim=2) == jpj ) THEN 1981 CALL xios_send_field( cdname, pfield4d(Nis0:Nie0, Njs0:Nje0,:,:) ) ! this extraction will create a copy of pfield4d 1982 ELSE 1983 CALL xios_send_field (cdname, pfield4d ) 1984 ENDIF 1964 CALL xios_send_field (cdname, pfield4d ) 1985 1965 #else 1986 1966 WRITE(numout,*) pfield4d ! iom_use(cdname) = .F. -> useless test to avoid compilation warnings … … 1994 1974 IF( iom_use(cdname) ) THEN 1995 1975 #if defined key_iomput 1996 IF( SIZE(pfield4d, dim=1) == jpi .AND. SIZE(pfield4d, dim=2) == jpj ) THEN 1997 CALL xios_send_field( cdname, pfield4d(Nis0:Nie0, Njs0:Nje0,:,:) ) ! this extraction will create a copy of pfield4d 1998 ELSE 1999 CALL xios_send_field (cdname, pfield4d ) 2000 ENDIF 1976 CALL xios_send_field (cdname, pfield4d ) 2001 1977 #else 2002 1978 WRITE(numout,*) pfield4d ! iom_use(cdname) = .F. -> useless test to avoid compilation warnings … … 2204 2180 ! 2205 2181 CALL iom_set_domain_attr("grid_"//cdgrd, ni_glo=Ni0glo,nj_glo=Nj0glo,ibegin=mig0(Nis0)-1,jbegin=mjg0(Njs0)-1,ni=Ni_0,nj=Nj_0) 2206 CALL iom_set_domain_attr("grid_"//cdgrd, data_dim=2, data_ibegin = 0, data_ni = Ni_0, data_jbegin = 0, data_nj = Nj_0)2182 CALL iom_set_domain_attr("grid_"//cdgrd, data_dim=2, data_ibegin = -nn_hls, data_ni = jpi, data_jbegin = -nn_hls, data_nj = jpj) 2207 2183 !don't define lon and lat for restart reading context. 2208 2184 IF ( .NOT.ldrxios ) & … … 2303 2279 CALL dom_ngb( 180.0_wp, 90.0_wp, ix, iy, 'T' ) ! i-line that passes near the North Pole : Reference latitude (used in plots) 2304 2280 CALL iom_set_domain_attr("gznl", ni_glo=Ni0glo, nj_glo=Nj0glo, ibegin=mig0(Nis0)-1, jbegin=mjg0(Njs0)-1, ni=Ni_0, nj=Nj_0) 2305 CALL iom_set_domain_attr("gznl", data_dim=2, data_ibegin = 0, data_ni = Ni_0, data_jbegin = 0, data_nj = Nj_0)2281 CALL iom_set_domain_attr("gznl", data_dim=2, data_ibegin = -nn_hls, data_ni = jpi, data_jbegin = -nn_hls, data_nj = jpj) 2306 2282 CALL iom_set_domain_attr("gznl", lonvalue = real(zlon, dp), & 2307 2283 & latvalue = real(RESHAPE(plat(Nis0:Nie0, Njs0:Nje0),(/ Ni_0*Nj_0 /)),dp)) 2308 CALL iom_set_zoom_domain_attr("ptr", ibegin=ix-1, jbegin=0, ni=1, nj=Nj _0)2284 CALL iom_set_zoom_domain_attr("ptr", ibegin=ix-1, jbegin=0, ni=1, nj=Nj0glo) 2309 2285 ! 2310 2286 CALL iom_update_file_name('ptr') -
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf.sh
r13733 r13781 18 18 # => i_sea_ice: whether to compute fluxes over sea-ice as well 19 19 # => SFORC: string sufficient to copy relevant files as in "*${SFORC}*.nc" 20 FORCING="PAPA" ; i_sea_ice=0 ; SFORC="Station_PAPA_50N-145W"21 #FORCING="ERA5_arctic" ; i_sea_ice=1 ; SFORC="ERA5_arctic_surface_81N_36p75E_1h" ; # "ERA5_arctic" WITH ice/air flux computation20 #FORCING="PAPA" ; i_sea_ice=0 ; SFORC="Station_PAPA_50N-145W" 21 FORCING="ERA5_arctic" ; i_sea_ice=1 ; SFORC="ERA5_arctic_surface_81N_36p75E_1h" ; # "ERA5_arctic" WITH ice/air flux computation 22 22 #FORCING="ERA5_arctic" ; i_sea_ice=0 ; SFORC="ERA5_arctic_surface_81N_36p75E_1h" ; # "ERA5_arctic" WITHOUT ice/air flux computation 23 23 -
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_ICE.py
r13723 r13781 36 36 ############################################ 37 37 crealm = 'sea-ice' 38 L_VNEM = [ 'Cd_ice', 'C e_ice', 'qla_ice' , 'qsb_ice' , 'qt_ice' , 'qlw_ice' , 'qsr_ice' , 'taum_ai' ]39 L_VARO = [ 'Cd' , 'C e' , 'Qlat' , 'Qsen' , 'Qnet' , 'Qlw' , 'Qsw' , 'Tau' ]40 L_VARL = [ r'$C_{D}$', r'$C_{ E}$', r'$Q_{lat}$', r'$Q_{sens}$' , r'$Q_{net}$' , r'$Q_{lw}$' , r'$Q_{sw}$' , r'$|\tau|$' ]38 L_VNEM = [ 'Cd_ice', 'Ch_ice', 'qla_ice' , 'qsb_ice' , 'qt_ice' , 'qlw_ice' , 'qsr_ice' , 'taum_ai' ] 39 L_VARO = [ 'Cd' , 'Ch' , 'Qlat' , 'Qsen' , 'Qnet' , 'Qlw' , 'Qsw' , 'Tau' ] 40 L_VARL = [ r'$C_{D}$', r'$C_{H}$', r'$Q_{lat}$', r'$Q_{sens}$' , r'$Q_{net}$' , r'$Q_{lw}$' , r'$Q_{sw}$' , r'$|\tau|$' ] 41 41 L_VUNT = [ '' , '' , r'$W/m^2$' , r'$W/m^2$' , r'$W/m^2$' , r'$W/m^2$' , r'$W/m^2$' , r'$N/m^2$' ] 42 42 L_BASE = [ 0.0005 , 0.0005 , 5. , 5. , 5. , 5. , 5. , 0.05 ]
Note: See TracChangeset
for help on using the changeset viewer.